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Abstract 



Superconductivity, discovered by Kamerlingh Onnes in 1911, continues to be a fas- 
cinating subject of condensed matter physics today Much interest has been devoted 
to the study of the superconductivity induced in a metal which by itself is not su- 
perconducting but is in electrical contact with a superconductor. As the carriers 
of superconductivity, the Cooper pairs, diffuse across the contact into the metal 
they remain correlated, although the pairing mechanism is lifted; we call this the 
proximity effect. The observation of these superconducting correlations has come 
within the reach of experiments in the last decade. With state-of-the-art fabrica- 
tion techniques mesoscopic samples have been produced which are small and clean 
enough for the quantum mechanical coherence of the electrons to be preserved over 
the sample size. This thesis focuses on the variety of signatures of single-particle 
physics that appear in the electrical transport and the magnetic screening properties 
of these systems. 

The main contribution of this thesis to the field of mesoscopic superconductivity 
consists in the study of the thermodynamic properties in the presence of a magnetic 
induction. We map out the phase diagram of a normal-metal-superconductor (NS) 
hybrid system in the space spanned by the temperature (T) and the magnetic field 
(H). The region of the H — T diagram of interest lies well below the transition tem- 
perature and the critical field line, i.e. in the Meissner phase, of the superconductor. 
The involved energy scales are determined by the geometric system size d and the 
mean free path I = vpT of the normal metal and are given by the single-particle 
level spacing ~ hvp/d or the Thouless energy ~ hvpl/d 2 , respectively, in ballistic or 
diffusive systems. 

We distinguish two phases in the H — T diagram: a diamagnetic state (at low 
temperature and field) where the applied field is effectively screened, and a nor- 
mal state where the magnetic field penetrates into the normal metal. The magnetic 
breakdown of the Meissner field expulsion occurring at the frontier of the two regimes 
is interpreted as the critical field of the proximity layer. We give a complete charac- 
terization of this first order transition for a ballistic NS structure, determining the 
spinodals, the magnetization jump, as well as the latent heat, and find a critical 
point, beyond which the two phases are connected by a continuous cross-over. 
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ABSTRACT 



We study the magnetic screening in the diamagnetic phase for a NS hybrid layer. 
We find that the current-field relation is typically non-local in the proximity effect 
producing a high sensitivity of the diamagnetic susceptibility to the mean free path. 
We classify the behavior into several regimes between the clean (/ — ► oo) and the 
dirty (7 — > 0) limit. The reduction of the screening density and the non-locality 
range with decreasing purity affect the screening in an opposite way, competing in 
the diamagnetic susceptibility observed in the experiments. 

The screening properties can be traced back to the quasi-particle excitation spec- 
trum in the absence of the field. We thus establish a connection to the transport 
properties of NS junctions: We explain the resonance structure in the conductance 
and the shot noise at finite voltage, which appears as a signature of the quantum 
behavior of the discrete electrons. 

Back to the H — T phase diagram: in the field penetration phase, where the 
screening is negligible, we analyze the geometry of a ballistic SNS Josephson link. 
We show that the combination of the non-locality and the finite size of the junction 
produces an anomalous critical-current-flux dependence. As a function of increas- 
ing field H we find a cross-over from the usual single to a double flux quantum 
periodicity. 

In the low-temperature low- field corner of the H — T phase diagram we discover a 
new paramagnetic phase by considering the effect of a repulsive electron-electron in- 
teraction in the proximity metal. Quasi-particles are found to be trapped at the NS 
interface, accumulating to a density of states peak at the Fermi energy. This peak 
induces a paramagnetic instability towards spontaneous magnetic moments. On 
top of the Meissner phase we find a first order transition as the moments polarize 
in the applied field. The hysteresis and the dissipation accompanying a paramag- 
netic reentrance in the susceptibility provide the characteristic features observed in 
fascinating experiments on NS cylinders. 

Finally, we consider the transport in Josephson junctions with unconventional 
superconductors, where the current-phase relation exhibits a double frequency. We 
propose to exploit this 7r-periodicity for the implementation of a quantum computer. 



Zusammenfassung 



Supraleitung, in 1911 durch Kamerlingh Onnes entdeckt, is heute noch ein faszinieren- 
des Gebiet der Festkorperphysik. Die Untersuchung der Supraleitung, welche in 
einem Metall induziert wird, das selbst nicht supraleitend ist, jedoch im elektrischen 
Kontakt mit einem Supraleiter steht, hat ein grosses Interesse erfahren. Wenn die 
Trager der Supraleitung, die Cooper-Paare, durch den Kontakt in das Metall dif- 
fundieren, behalten sie ihre Korrelationen, obwohl der Paarungsmechanismus ver- 
schwindet. In der letzten Dekade ist die Beobachtung der supraleitenden Korrelatio- 
nen zwischen einzelnen Quasi- Teilchen in die Reichweite der Experimente geriickt. 
Die modernsten Fabrikationstechnologien haben die Herstellung von mesoskopischen 
Systemen ermoglicht, welche klein und sauber genug sind, um die quantenmechani- 
sche Phasenkoharenz der Elektronen iiber die Probengrosse zu gewahrleisten. Diese 
Dissertation widmet sich der Vielfalt der Einteilchen-Physik, welche im elektrischen 
Transport und der magnetischen Abschirmung in diesen Systemen auftritt. 

Der Hauptbeitrag dieser Dissertation zum Feld der mesoskopischen Supraleitung 
besteht im Studium der thermodynamischen Eigenschaften unter dem Einfluss einer 
magnetischen Induktion. Wir bestimmen fur eine Metall-Supraleiter (NS) Hybrid- 
struktur das Verhalten im Phasendiagramm, welches durch die Temperatur (T) und 
das Magnetfeld (if) aufgespannt wird. Der untersuchte Bereich im. H — T Dia- 
gramm liegt deutlich unter der Ubergangstemperatur des Supraleiters sowie der 
kritischen-Feld-Grenze. Die geometrischen Masse d und die freie Weglange I = vpr 
bestimmen die typische Energieskala im Problem: die Einteilchen-Energie ~ hvp/d 
beziehungsweise die Thouless-Energie ~ hvpl/d 2 in ballistischen oder diffusiven Sys- 
temen. 

Wir unterscheiden zwei Phasen im H — T Diagramm: einen diamagnetischen 
Zustand (bei kleinen Temperaturen und Feldern), der sich durch die Abschirmung 
der angelegten Felder auszeichnet, und einen normalleitenden Zustand, wo das Mag- 
netfeld ungehindert in das Metall eindringt. Wir bestimmen das kritische Feld, das 
die beiden Phasen trennt, und charakterisieren diesen Ubergang erster Ordnung in 
einem ballistischen NS Schichtsystem durch die Spinodalen, den Sprung in der Mag- 
netisierung sowie die latente Warme. Wir finden, dass die Phasengrenze oberhalb 
eines kritischen Punktes verschwindet. 
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ZUSAMMENFASSUNG 



In der Beschreibung der diamagnetischen Phase finden wir eine fur die induzierte 
Supraleitung typische nicht-lokale Strom-Feld Relation, welche die diamagnetische 
Suszeptibilitat zu einem feinen Indikator der freien Weglange macht. Wir unter- 
scheiden mehrere Verhalten zwischen dem reinen (/ — > oo) und dem schmutzigen 
(7 — > 0) Grenzfall. Mit abnehmender Reinheit haben die Reduktion der abschir- 
menden Dichte sowie der Reichweite der Strom-Feld Relation eine entgegengesetzte 
Wirkung auf die Abschirmung und stehen daher im Wettbewerb in der experimentell 
gemessenen diamagnetischen Suszeptibilitat. 

Das Abschirmverhalten kann auch auf das Quasiteilchen-Spektrum zuruckgefuhrt 
werden. Dies ergibt eine Verbindung zu den Transport-Eigenschaften in NS Struk- 
turen: Wir untersuchen die Resonanzen im Leitwert und finden sie wieder in den 
quantenmechanischen Stromfluktuationen aufgrund der Diskretheit der Elektronen- 
ladung. 

In der normalleitenden Phase, wo die Abschirmung vernachlassigt werden kann, 
untersuchen wir die Anordnung eines ballistischen SNS Josephson Kontaktes. Wir 
zeigen, dass die Verbindung von Nichtlokalitat und der endlichen Ausdehnung des 
Metalls zu einer anomalen Abhangigkeit des kritischen Stromes vom magnetischen 
Fluss fiihrt. Mit zunehmendem Feld finden wir einen Ubergang von der einfachen 
Periode eines Fluss-Quantums zur doppelten. 

Wir entdecken die Existenz einer neuen paramagnetischen Phase, welche durch 
eine repulsive Elektron-Elektron Wechselwirkung verursacht wird, in der Nahe des 
Ursprungs des H — T Phasendiagramms. Quasiteilchen werden an der NS Grenz- 
flache festgehalten und tragen zu einer Spitze in der Zustandsdichte bei der Fermi- 
Energie bei. Diese verursacht eine paramagnetische Instabilitat, welche spontane 
magnetische Momente hervorbringt. Dem diamagnetischen Verhalten uberlagert fin- 
den wir das die Polarisation der magnetischen Momente zu einem Phaseniibergang 
erster Ordnung bei Feld Null fiihrt. Die Hysterese und Dissipation, welche das para- 
magnetische Signal in der Suszeptibilitat begleiten, sind in qualitativer Ubereinstim- 
mung mit bisher noch unerklarten Experimenten auf mesoskopischen Zylindern. 

Schliesslich betrachten wir den Transport in Josephson Kontakten mit unkonven- 
tionellen Supraleitern, wo eine doppelte Josephson-Frequenz auftritt. Wir schlagen 
vor, diese 7r-Periodizitat in der Phase zur Implementation eines Quanten-Computers 
zu benutzen. 



Chapter 1 
Introduction 



In this thesis we study the behavior of electrons in a metal which is in electrical 
contact with a superconductor. The superconductor (S) can be viewed as a reservoir 
where due to an effective attraction the electrons condense into Cooper pairs, while 
the normal metal (N) is idealized as a gas of free electrons. Through the diffusion 
across the NS contact the metallic electrons acquire superconducting correlations; 
we call this the proximity effect. 

The Cooper pairs are essentially described by one wavefunction (the order pa- 
rameter of the superconductor) producing quantum behavior on a macroscopic scale. 
The quantum interference effects between the free electrons instead rely on the co- 
herence of the single-particle wavefunctions and are much more fragile. The observ- 
able quantum phenomena are limited to the mesoscopic scale given by the phase 
coherence length here. The proximity effect offers an experimental access to the co- 
herent single-particle physics by combining it with the macroscopic superconducting 
correlations. 

Here we study the electrical transport and the orbital magnetism of electrons in 
normal-metal - superconductor hybrid structures. Our theoretical work is motivated 
by an intense experimental activity in this area of mesoscopic physics triggered by 
recent advances in nanostructure fabrication technology. For an introduction to 
the subject we refer the reader to the review on "Mesoscopic Electron Transport" 



Kouwenhoven 1997 



The microscopic process underlying the physics of proximity induced supercon- 
ductivity is the Andreev reflection (AR) at a NS interface: An electron incident 
from the normal metal into the superconductor is retro-reflected into a hole, trans- 
mitting the double electron charge into the superconductor [[Andreev 1964j| . In the 



following, we give an introduction to the basic properties of the AR and relate the 
works of this thesis to this fundamental process. 
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CHAPTER 1. INTRODUCTION 



The quasi-particle transfer at a NS interface is best described by the Bogoliubov- 



de-Gennes (BdG) equations ||de Gennes 1966a] , which determine the two-component 



wavefunctions $ of electron and hole close to the Fermi level, 

h o~V A \f u \_,f u \ $=f M )e ik ^ x (1.1 ; 



A* — h* Q + /i J \ v J \ v J ' ' \ v 

We use the linearized Hamiltonian ho — \i — — iftvpV, where vp = hkp/m is the 
Fermi velocity, A provides the off-diagonal pair-potential and e is the quasi-particle 
energy with respect to the chemical potential. Let us consider a NS junction as 
described by the pair potential 



A e^ x > (S) 
While in the normal metal the electron and hole states decouple , 

$+ = ^ J J e ifc+x , = ^ ° J e ifc " x (1.3) 

(fc ± = k x ±e/h~v x , hk x /m = v x = t^cost?), the quasi-particles in the superconductor 
are coupled evanescent modes, 

(7 = exp[— i arccos(e/A)], k e> h = k x ± «a/Aq — e 2 /Hv x , e < A). By using the ansatz 
<fr + +rA*&- for x < and $ e for x > we easily find the Andreev reflection amplitude 
from the continuity of the wavefunction, 

r A = = exp(— z arccos ^) e - * 1 ^. (1.5) 

For e < A, the incident quasi-particle excitation is reflected with probability \r^\ = 
1. Similarly the Andreev amplitude describing the conversion of an incident hole to a 
reflected electron is given by jexp(iip). Andreev thus found that the incident quasi- 
particle excitation is prevented from entering the superconductor, and consequently 
heat (or energy) transport across the NS interface is suppressed | |Andreev 1965| . 
On the other hand, the double charge 2e (of incident electron and reflected hole) 
is transferred across the NS interface. In the superconductor the decaying quasi- 
particle current is converted into supercurrent ||Blonder 1982|. 



The Andreev reflection (AR) coherently couples electron and hole states and 
thus induces the superconducting correlations F(x) in the normal metal, 

F{x) = (^(x)* T (s)) = Y, u a (x)v* a (x)f(e a ), (1.6) 



3 



where the index a denotes the eigenstates and /(e) is the Fermi occupation number, 
see also appendix |A]. Instead of a Cooper pair condensate, the proximity metal 
exhibits correlated electron- hole pairs, which are mutually independent. Here lies 
the fragility of the proximity effect: The thermally excited Andreev pairs at energy 
e ~ T are dephased over the thermal length £at(T) = hvp/^irT and the correlation 
function (Of) is suppressed as ~ exp[— x/£n(T)] (see Chp. |3T2]). 



Let us consider in the following a normal metal slab of thickness d on top of a bulk 
superconductor. Due to the AR process, the quasi-particles at subgap energies e < A 



are bound to the normal layer, see Fig. 1.1. The discrete energies of the Andreev 



bound states follow from a constructive interference condition for an electron and a 
hole traveling back and forth across the normal layer and being converted into each 
other at the NS interface (n = 0, 1, 



2k + d — 2k_d — 2 arccos ^ = 2ri7r, (1-7) 
e n (<) = —(n7r + arccos—) w ~l^f( n + ( L8 ) 

The bound states are quantized at energies ~ hv x /d depending only on the compo- 
nent of the Fermi velocity in forward direction (v x — Vp cos$). The energy levels are 
repelled from the chemical potential, as a finite energy mismatch of the wave vector 
k± of electron and hole is needed to compensate for the phase arccos(e/A) w n/2 
of the AR. The discrete Andreev levels contribute to the finite density of states at 
e < A filling the superconducting gap. 

The transport through the ballistic NS structure coupled by means of a tunneling 
barrier (I) to a NINS junction, see Fig. |1.1| , reveals a resonance structure at the 
energies of the Andreev (quasi)-bound states e ~ hvp/d [ [Rowell 1966|| . A similar 



behavior is expected in a diffusive system, where the average path traveled d is to 
be replaced by d 2 /l (I mean free path). According to the constructive interference 
condition ed 2 /FlvfI ~ (2n + 1)tt the resonances are found at the Thouless energy 
E c ~ hvpl/d 2 . In chapter |]we analyze the spectral conductance G(e) and the shot 
noise S(e) in NS junctions with an arbitrary scattering region in the normal lead, 
which typically includes a barrier at the NS interface and additional disorder in the 
normal lead. We describe the single-particle resonance structure that follows from 
the interplay between the normal scattering and the AR processes and enhances the 
subgap transport. 

The linear transport coefficients G(0) and S(0) in a NS junction are suppressed 
by the transmission T 2 of the N conductor, reflecting the two-particle nature of 
the AR process [[Blonder 1982| , |Beenakker 1992a| |de Jong 1994|| . This explains why 



the first observation of an enhanced subgap conductance could be misinterpreted 
as being due to supercurrents [[Kastalsky 1991||. Actually, the iterative scattering 
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CHAPTER 1. INTRODUCTION 



AR 



AR 




Figure 1.1: Trajectories of electron- hole Andreev levels in a proximity normal-metal 
layer: electron and hole travel back and forth and are converted into each other by 
the Andreev reflection (AR). Top: Coupling of the NS slab by a tunnel contact to 
a NINS junction. 



processes between the disorder and the interface barrier increase the possibility 
for an AR as compared to a simple tunneling junction [[Melsen 1994|| . If, on the 
other hand, the interface barrier is weak, the resonance structure of the Andreev 
states, which are el-hole quasi-bound states between the disorder and the super- 
conductor, determine the transport coefficients. The conductance G(e) as well 
as the shot noise S(e) exhibit peaks at the Andreev or the Thouless energy, re- 
spectively, for a ballistic or a diffusive conductor. Characteristically these reso- 
nances are pinned at a finite energy above the Fermi level, which we understand in 
terms of the ir/2 phase shift in the AR process, and survive in the multichan- 
nel limit. The reentrance effect observed in NS junctions with a diffusive lead 
[Pctrashov 1995| , |Charlat 1996| , [Poirier 1997a|| reflect the non-monotonicity of the 
spectral conductance due to the underlying Andreev resonances. While theoreti- 
cal works using the Green's function technique have explained these experiments 
accurately [|Yip 1995| , |Nazarov 1996| , |Stoof 199q , |Colubov 19971 , IWilhelm 19981 , the 
Landauer-Biittiker approach offers an intuitive understanding here, based on a bal- 
listic viewpoint. The diffusive system is described by averaging over an ensemble 
of impurity configurations, see | |Marmokos 1993| , |Argaman 1997| [Lambert 1998|| . In 
Chp. we present new results for the shot noise, which are of interest with respect 



to ongoing noise experiments in NS junctions |Jehl 1999, Kozhevnikov 1999 



The quantization of the Andreev levels at discrete energies away from the Fermi 
level suppresses the density of states (DOS) of the proximity layer as compared to 



5 



the bulk metal. From the low energy spectrum (|1.8|) we obtain the Andreev DOS 



de Gennes 1963 



N(E) = y *h^ S [E - eo{v w )] (1.9) 
~ N -^- (1.10) 



(v x = h^J kp — ky — kg/m, N = mk F /h 2 Tx 2 ). The excitation gap A is softened in 

the normal metal to a linearly vanishing DOS below Hvp/d. We note that similar 
results for the DOS are observed in diffusive NS hybrid structures ||Gueron 199(3 . 
Belzig 19961 . 



The DOS suppression close to the Fermi level affects the response properties 
of the proximity layer. In chapter ||| we consider the magnetic response, i.e., the 
screening currents induced by an applied magnetic field. Generically, the linear 
current response j = jdia + jpara is given by the sum of the diamagnetic current 
describing the response of the rigid wavefunction and the paramagnetic current 
induced by the perturbation of the wavefunction in the applied field. While the 
diamagnetic current is trivially proportional to the bulk electron density, the para- 
magnetic current depends on the DOS at the Fermi surface |Schricffcr 1988] . The 



solution of the screening problem requires the full knowledge of the dispersive re- 
lation j(q) = — (c/X 2 )K(q)A(q) which has to be solved together with the Maxwell 
equation j(q) = q 2 A(q). (A = (Anne 2 /mc 2 ) -1 / 2 is the London penetration depth). 
The variation of the field on wavelength 1/q typically induces quasi-particle excita- 
tions at energy E ~ hqvp giving a paramagnetic current. Depending on the avail- 
ability of the finite energy excitations, the paramagnetic current either vanishes or 
compensates for the diamagnetic current. 

In a BCS superconductor with excitation gap A, the paramagnetic current is 
gapped (the DOS vanishes) for long wavelength and the response is purely diamag- 
netic, K(q) ~ 1 for hqvF < A. Above the gap the paramagnetic current is found 
to cancel the diamagnetic current leaving a small response K(q) ~ A/hqvp. In real 
space this implies a range of the current-field relations of the order of £o ~ hvp/A- 
For a London length A > £o we find a local response, while for A < £o the current- 
field dependence is nonlocal according to a Pippard-type relation [ [Fippard 1953|| . In 



a bulk normal metal, the continuous density of states at the Fermi level implies a 
cancellation of dia- and paramagnetic currents, leaving no net screening response. 

In a clean proximity layer, the suppression of the DOS produces a net diamag- 
netic current, but the lack of a quasi-particle gap translates to a diverging nonlocality 
range. It turns out that for finite q the paramagnetic current is able to compensate 
the diamagnetic current, producing the kernel K(q) ~ S(q), and thus the nonlocal 
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CHAPTER 1. INTRODUCTION 



current functional llZaikin 1982 



3{x) 



\ 2 d 



dx'A(x') 



Xll) 



We study this screening problem in more detail in Clip. ^, where the magnetic 
response of a NS slab with arbitrary impurity concentration is discussed. We intro- 
duce the powerful quasi-classical Green's function technique to derive the general 
linear current functional. We find that a finite mean free path provides a cutoff for 
the nonlocality range, but that the current-field relations in the proximity effect 
nevertheless remain typically nonlocal. We distinguish several regimes for the mag- 
netic response depending on the geometric thickness d, the thermal length £,n{T) 
and the mean free path I. The reduction of the screening density and the nonlocality 
range with decreasing purity affect the susceptibility in different ways, producing 
a non-monotonic dependence in the diamagnetic susceptibility with changing mean 
free path. This understanding provides the basis for a quantitative fit with the ex- 
periments ||Muller-Allinger 1999| . 



The nonlocal field dependence is also found in the spectrum of the Andreev 
states. Replacing the gradient in ( |1 . 1| ) by V + ieA(x)/hc we can account for the 
vector potential by adding an Aharonov-Bohm-type phase to the constructive inter- 
ference condition (|1.7|), yielding the bound state energies 



hv x 
~Ad 



2n f 

(2n + 1W <h Adx 

% J 



;i.i2) 



The circular integral gives the flux $ = § Arfx enclosed between the quasi-classical 
trajectory and the NS interface, see Fig. |TT] (hatched region). The flux enters 
modulo an integer number of flux quanta $o = hc/2e. This leads to a random 
distribution of the Andreev levels at large field $ ^> $o? suppressing the energy 
gap and inducing the magnetic breakdown in the proximity layer. The magnetic 
breakdown denotes the proximity layer's own critical field H^, at which the Meissner 
field expulsion ceases to be effective and the applied field penetrates into the normal 
layer by a first order transition. While in a small applied magnetic field the Andreev 
pairs produce diamagnetic screening currents, at large fields the superconducting 
correlations are destroyed by the random Aharonov-Bohm phases in ( |1.12| ) and the 
net screening current vanishes. In Chp. f| we determine the H — T phase diagram 
for this behavior, finding the first order transition line separating the diamagnetic 
phase of effective field expulsion and the normal state of field penetration. 

The magnetic breakdown has been observed recently in increasingly clean sam- 
ples [ pda 19801 , |Mota 1982] , Pergmann 1987| , |Visani 1990b|| and we find good agree- 
ment of both temperature and thickness dependence of the breakdown field with ex- 
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Figure 1.2: Trajectories of electron-hole Andreev levels in a normal-metal - su- 
perconductor - normal metal (SNS) junction. Note the existence of trajectories 
reflecting at the metal boundary. 



perimental data on quasi-ballistic Ag-Nb samples ||Mota 1989[ |Miiller-Allinger 1999 . 



We study a different aspect of nonlocality in Chp. [|, where we consider a ballistic 
normal metal layer connected by two superconducting reservoirs to a SNS junction. 
The two NS interfaces act as a confining potential for the electrons and holes in the 
normal interlayer. The junction exhibits current- carrying bound states consisting of 
an electron traveling forward and a hole traveling backwards which are converted 
into each other by the AR |[Kulik 19~T0 



see Fig. 1.2. The constructive interference 



condition, 



(k+ - k^)d =f A<p = {2n + 1)tt, 



;i.l3) 



now features the phase difference Aip of the two superconducting reservoirs (the sign 
of Af depending on the orientation of the quasi-particle current) and implies the 
spectrum of bound states 



(mr H — ± Aip). 

d K 2 yi 



;i.i4) 



These bound states are known to transmit the supercurrent between the two reser- 
voirs: While at phase difference Aip = the left- and right-going states are degener- 
ate and their respective currents cancel, by changing A(p an asymmetry is introduced 
producing the supercurrents in the junction. 

Including a magnetic induction as above, the bound state energies 



~d 



7T 



±[Acp 



2tt 

$ 



A rfx 
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Figure 1.3: Trajectories of representatives from two Andreev level populations 
shifted by a phase difference tt. 



depend on the gauge invariant phase difference (7) between the two endpoints 1 
and 2 of the quasiparticle trajectories connecting the superconducting reservoirs, 
see Fig. [L2|. The currents j oc dE/dj depend non-locally on the field A(x). The 
Josephson current through the junction is given by the sum of the currents carried 
along all the quasi-particle trajectories across the junction. In Chp. ^| we find that 
in a junction of finite width the trajectories with multiple reflections at the metallic 
boundary alter the field dependence of the Josephson current, as they experience 
a different gauge invariant phase difference, see Fig. |1.2| . The current-flux relation 
/ = /($) crosses over from the usual $0- to the anomalous 2$ -periodicity as a 
function of increasing field. The cross-over takes place when the Josephson vortex 
distance a = Q /Hd becomes smaller than the nonlocality range of the junction. 
Our results explain the experimental data recently obtain on S-2DEG-S junctions 
Heida 199811 , where the ballistic link was provided by a two-dimensional electron 



gas (2DEG). 



Let us consider the Josephson relation / = I(A(p) in a junction between a con- 
ventional, S-wave superconductor (S) and a High-T c superconductor with D-wave 
symmetry (D). The order parameter in the D-wave superconductor is direction- 
dependent as sin 2$ exhibiting two lobes with positive sign and two lobes with 



negative sign, see Fig. |1.3| . The phase difference Aip across the SND junction thus 
depends on the orientation of the quasi-particle trajectory, varying by 7r for tra- 
jectories connecting with the negative rather than the positive lobes of the D-wave 
material. For the alignment of the D-wave superconductor with a gap A oc sin 2$ 
with respect to the x-axis, the two quasiparticle populations with d > and $ < 0, 
see Fig. |1.3j , correspond to two conventional Josephson junctions shifted by tt with 
respect to each other. The spectrum of the bound states is given by Eq. ( |1.13| ), 
where Aip differs by 7r for the two quasi-particle populations. Therefore, the 2tt 
periodicity of the critical current of each junction is reduced to a 7r periodicity in 
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the SND junction. 

In chapter [7| we propose to exploit the doubly periodic Josephson relation for the 
implementation of a quantum bit (qubit). The 7r-periodic junction can be viewed as 
a macroscopic two-level system exhibiting two degenerate ground states which carry 
no current and are thus optimally isolated from the environment. The macroscopic 
quantum tunneling between these two states opens the possibility for a quantum 
time evolution of the junction, which we propose to use for quantum computation. 



In chapter [5] we return to the orbital magnetism, synthesizing various elements 
of proximity induced superconductivity discussed above. The motivation for this 
work has been provided by the puzzling low- temperature anomalies observed in 



normal-metal - superconductor cylinders by Mota and co-workers [Visani 1990a 



Mota 1994 . In the low-temperature low-field corner of the H-T phase diagram, the 
Ag-Nb and Cu-Nb samples develop a paramagnetic signal on top of the saturated 
Meissner response, which produces a reentrance in the magnetic susceptibility (for a 



detailed discussion of the experiments see also | |Visani 1990c| , |Frassanito 1996| ). We 



note that the observed behavior is not believed to be related to the low-temperature 
reentrance in the conductance of a disordered NS junction discussed above, which 
is due to the finite energy properties entering the non-equilibrium transport. The 
magnetic response instead is a thermodynamic property of the system. 

In chapter |5| we consider the consequence of a repulsive electron-electron interac- 
tion in the normal metal proximity layer. The repulsive coupling > induces a 
finite gap A(x) = V^F{x) in the normal layer which is proportional to the induced 
superconducting correlations ( |1.6|) . The order parameter changes sign across the 
NS interface as the coupling constant changes from attractive to repulsive. The NS 
structure behaves like a Josephson junction with a phase difference of tc across the 
NS interface. We know from Eq. ( |1.13| ) that a junction with Atp = tc traps density 
of states at the Fermi energy e = 0, according to 

(k + -k_)L = 2nn, (1.16) 

which remains true in the limit of a vanishing junction length L — > 0. The phase 
difference 7r across the junction cancels the phase shifts 7r/2 of the two AR. As all 
transverse wavefunction (k y , k z ) exhibit a zero energy bound state, the DOS shows 
a peak a zero energy with macroscopic weight 

N(E) ~ N k 2 F 5(E), (1.17) 

per unit surface. 

As we show in Chp. [|, this peak implies a paramagnetic instability in the current 
response and leads to a spontaneous onset of currents along the interface. Sponta- 
neous moments form at the NS interface producing a low-temperature reentrance 
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in the magnetic susceptibility. The qualitative agreement of our results with the 
observed low-temperature anomaly is indicative of a repulsive electron-electron in- 
teraction in these systems. 



Chapter 2 



Nonlinearity in normal-metal — 
superconductor transport: 
Conductance and shot noise 



2.1 Scattering matrix approach to normal- metal 
— superconductor transport 

2.1.1 Introduction 

The study of electronic transport in normal-metal-superconductor (NS) or semicon- 
ductor-superconductor (SmS) sandwiches has attracted a considerable amount of 
interest in the past years. At sufficiently low temperature and in high quality 
mesoscopic samples, the phase-breaking length of the electrons is larger than the 
typical system size, resulting in directly observable quantum coherence effects. Of 
special interest is the effect of the electronic phase-coherence in a normal metal- 
superconductor system. In the standard theory of the proximity effect, the influence 
of the superconductor on the normal metal can be understood in terms of the coher- 
ent coupling of electrons and holes in the metal as described by the Bogoliubov-de 
Gennes (BdG) equations ||de Gennes 1966a[| . The correlation between electrons and 



holes is produced by the process of Andreev reflection [[Andreev 1964j , [Blonder 1982 ] 
at the NS interface, by which an incident electron is back- reflected into a hole, con- 
verting quasiparticle current to supercurrent. This microscopic picture of electron- 
hole correlation is equivalent to that provided by a pair correlation function in the 
normal metal, which is induced by the superconductor through the contact at the NS 
interface. The scattering matrix approach makes use of the microscopic single par- 
ticle picture of coupled electron and hole channels providing a straightforward tool 



within a formalism of the Landauer-Biittiker type [[Landauer 1970| , [Biittiker 1985| . 
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Beenakker 1992aj 



The quality of the interface as well as the phase breaking processes determine the 
strength of the proximity effect and naturally have their impact on the current-vol- 
tage characteristics (CVC). A few fascinating transport experiments ||Pothier 1994 , 



Petrashov 1995|, |(Jourtois 199t}| have been carried out recently, investigating temper- 



ature and voltage dependence, as well as the flux modulation, of both NS and SmS 
junctions. Interestingly, the relative strength of the interface barrier and the elastic 
scattering in the normal region is crucial for the features of the CVC. The ratio of 
the two determines whether subgap conductance peaks arise at zero or finite voltage 
Volkov 1993|, [Marmokos 1993|, |Yip 1995] . The investigation of these so-called zero 



and finite bias anomalies in the subgap conductance have been the object of recent 
experiments |[Kastalsky 199Tj |JNguyen 199^ , [INitta 1994| , |Bakker 1994 |Magnee 1991] , 
Poirier 1997a|| . The present work draws much of its motivation from the ongoing 
discussions and experiments in this area, see also [[Aleiner 1996|1 - In our case study 



of a double barrier NINIS junction we observe zero and finite bias anomalies and 
shed light on the underlying mechanism. 

Recently in [[Nazarov 1996| , |Stoof 1996| as well as fColubov 1997|| , the carrier 
transport in disordered NS junctions has been described in terms of an energy- 
dependent diffusion constant, successfully explaining the recent experiments on 
reentrance in the conductivity at low temperatures [ Petrashov 1995 , Charlat 1996 ] . 
Their work uses the quasi-classical Green's functions technique [[Rammer 1986| , 
which allows to describe transport both close to equilibrium [[Hekking 1994|| and 
far away from equilibrium | ]Volkov 1993| , | Volkov 1994| , ^hou 19915] and facilitates the 
averaging over disorder in diffusive conductors. This approach is quite indispensable 
if phase-breaking processes are to be included. An appealing alternative approach 
is the scattering matrix technique, which relies on the quasiparticle wavefunctions 
described by the BdG equations. While being valid in a general context, it de- 
scribes the transport in mesocopic systems from the intuitive ballistic point of view. 
The transport through normal or superconducting leads is expressed through the 
properties of a multichannel scattering matrix accounting for all elastic scattering 
processes, whether they be due to a geometric constriction, single impurities, or dis- 
order (inelastic processes are excluded from such a description). By these means, the 
transport problem is reduced to solving a ballistic problem at the interfaces of the 
normal and superconducting leads. The current and the conductance of the system 
are determined analytically in terms of the transmission and reflection amplitudes of 
the scattering matrix. Adhering to this formalism rather than the Green's function 
technique helps us to improve our understanding of zero and finite bias anomalies. 

The study of normal-superconducting junctions goes back to the works of Kulik 
Kulik 1970|j on SNS junctions and of Blonder et al. fBlonder 1982|| on NIS junctions, 



who studied nonlinear transport within the framework of the BdG equations consid- 
ering quasi one- dimensional models. On the other hand, the scattering matrix tech- 
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nique was developed by Landauer and co-workers ||Landauer 1970| , [Biittiker 1985 ] 
in the linear response regime, resulting in the well known conductance formula for 
a normal metal. Lambert [ [Lambert 1991| and Takane and Ebisawa flakane 1991 



Takane 1992|| extended the approach to include superconducting segments, on the 
basis of which Beenakker [ [Beenakker 1992a| derived a zero temperature, linear re- 
sponse conductance formula for the transport through NS junctions. A few stud- 
ies have been carried out recently [[Marmokos 1993| , [Brouwer 1 9 9 5| , |Claughton 1996 



combining the scattering matrix approach with the finite voltage transport model 
of Ref . [Blonder 1982|| . They have limited themselves to the energy dependence of 
the scattering states to extract finite voltage properties of the CVC in the subgap 
regime. Here, we extend these works to voltages above the gap and additionally take 
into account the full voltage dependence of the transport problem. This is impor- 
tant within the context of the sign reversal symmetry of the differential conductance 
which we discuss below. At the same time, we provide a common framework for the 
above studies, tracing them back to a single general formula. 

In the present chapter, we derive the general expression for the current through 
a NS junction in the scattering matrix approach, valid for multiple channels, finite 
voltage, and nonzero temperature. We review the derivation of the current-voltage 
relation and express it in terms of a spectral conductance formula, thereby account- 
ing for the full voltage dependence of the transport problem. The reflection at the 
NS interface is made explicit using the Andreev approximation and a conductance 
formula is obtained expressing the result in terms of the normal scattering matrix. 
We illustrate this formula in the new regime of voltages above the gap and expose 
its connection to previously obtained limits. We describe the existence of resonances 
due to quasi-bound Andreev states, and show that they produce sharp conductance 
peaks in both the single and multichannel junction. We explain the generic mech- 
anism underlying the appearance of zero and finite bias anomalies in the ballistic 
two barrier system. Finally, we interpret our results in connection to the recent 
experiments on zero and finite bias anomalies in disordered NS junctions. 

Chapter ^TT] closely follows our paper published in ||Lesovik 1997| |. We would 
like to point out further works, both theoretical and experimental, which have been 
carried out since then [tiartog 1997, Lambert 1998, Wilhelm 1998| , Lesovik 1998 . 



We consider a normal-superconducting junction with quasi one-dimensional, bal- 

The pair potential 



listic normal and superconducting leads, as shown in Figure 2.1 



vanishes in the normal part, due to the absence of electron-electron interactions. 
Between the normal lead and the NS interface, the electrons traverse a disordered 
region, the transmission through and the reflection thereof are described by a scat- 
tering matrix. Between the scattering region and the NS interface, a small ballistic 
normal region serves to separate the scattering in the normal part, which mixes all 
electron channels at a given energy, from the scattering at the NS interface, where 
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Figure 2.1: Schematic structure of a disordered NS junction. Ballistic normal (Ni, 
N 2 ) and superconducting leads (S) are coupled to reservoirs at chemical potential 
H — eV and /x, respectively. Scattering is limited to the hatched region between leads 
Ni and N 2 . 



electron and hole channels are mixed in the reflection process. The coherent scatter- 
ing in the disorder region of the normal metal is described by the 4N x 4N scattering 
matrix, 



( ot\ 

\O k 2 ) 



( r n (e) t l2 (e) 
r* u (-e) 
hi (e) r 22 (e) 





(- 




t*u H) 



*Si H) 







22 



H)/ 



\ Ih 2 J 



ni(e) £i 2 (e 
r22(e 



(2.1) 



which we denote by Sat. The matrix connects the incoming N electron (hole) chan- 
nels If [If) on each side to the equal energy outgoing channels Of (Of) according 
to Figure [2.1| (i = 1,2). The channels represent the different transverse states, 



rn and are N x N reflection and transmission matrices for electron channels, 
fjj and ijj the comprehensive 2N x 2iV matrices including the complex conjugated 
reflection and transmission amplitudes for holes []. Following usual convention, we 
include the propagation in the ballistic region N 2 in the scattering matrix. For states 
normalized to carry unit probability current [ [Blonder 1982| ], the continuity equation 
implies the unitarity of the scattering matrix. By allowing for a dependence of the 
scattering matrices on voltage we can account for the full voltage dependence of the 
scattering problem. We can thus describe the deformation of the states due to the 



1 We denote the complex conjugate of a matrix m by m* , the transposed matrix by m T and the 
adjoint by rrv. 
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space dependent potential including, e.g., the voltage dependent Schottky barrier at 
a semiconductor-superconductor (SmS) interface. 

We define an analogous unitary scattering matrix for the NS interface S/ (see 
Fig. O) by 



1 2 

Th 

1 2 

o% 



ri(e) ifj(e 
ii(e) f f j(e) 



( °S \ 

Is 



(2.2) 



The 2N x 2N matrices 77, f'j, ij, i\ give the reflection and transmission amplitudes of 
the states normalized to unit probability current in the normal and superconducting 
leadQ. In the evaluation of the current, we will consider the matrix ( |2.1|) to be 
specified by an arbitrary model for the disorder, while the matrix ( |2.2| ) will be made 
explicit below using the BdG equations. The total effect of all scattering processes 
in the disorder region and at the NS interface can be described by a unitary, global 
AN x AN scattering matrix. We restrict ourselves to its sub-matrix f describing the 
reflection into the normal region, 



Of 

o\ 



r(e,V) 



1 i 



r ee (e,V) r eh (e,V) 
rhe(t,V) r hh (e,V) 



1 i 



(2.3) 



r eh, r he, and r hh are again N x N reflection matrices, f (e, V) can be expressed 



(2.4) 



through the given scattering matrices Q2.1|) and (|2 



f (e, V) = hi (e) + tia (e) [1 - fj (e) f 22 (e) ] X f 7 (e) t 21 (e) 



representing the sum over all scattering paths of an incident electron or hole exci- 
tation, multiply scattered between the disorder region and the NS interface. Apart 
from the direct reflection at the disorder region, the simplest process contributing 
consists of an excitation, which is first transmitted (£21) through the disorder region, 
reflected (77) at the NS interface, and transmitted (£12) back to the normal lead. All 
further paths result from iterative scattering processes between the disorder region 
and the NS interface. 

We now derive the current in the normal lead, expressing it through the global 
reflection matrix f(e, V). Applying a voltage V (denoting the voltage in a two 
point measurement) on the normal side has two consequences. First, the voltage 
induces an electrostatic potential drop over the disorder region in the NS junction, 
resulting in a deformation of scattering states. The coupling of incident and outgoing 
channels is thus voltage dependent in general, as described by f (e, V). A stationary 

2 Note that the connection the unit charge currents on either side of the interface will not produce 
unitarity, since quasiparticle charge current is not a conserved quantity in the superconductor, see 
also appendix [A] 
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state incident from the normal lead (of energy e in channel v) consists of the incident 
electron and the reflected electron- and hole- states and carries the current (e = |e|), 



L{e, V) de 



h 




(2.5) 



Second, the applied voltage shifts the chemical potential of the reservoir attached to 
the normal lead by — eV with respect to the reservoir on the superconducting side. 
The deformation of the states by itself produces no net current^. The net current flow 
results exclusively from the difference in occupation of the (finite voltage) scattering 
states incident from the left and right reservoirs. Writing the sum over channels as 
a trace, we obtain the current-voltage relation 



de - 

e 



f(e)-f(e + eV) G s (e, V) 



with the spectral conductance 

1 - rl (e, V) r ee (e, V) + r{ e (e, V) r he (e, V) 



2e 2 

G s (e, V) = -^Tr 



(2.6) 



(2.7) 



A factor two accounts for the spin degeneracy. The defined spectral conductance 
G s (e, V) describes the current contribution of the incident scattering states at energy 
e, at a given voltage V (by convention, the energy is measured with respect to 
the chemical potential in the superconductor). Formulas fl2.6| ) and (|2.7| ) imply the 
differential conductance 



dl 

dV 



J def'(e + eV) G s (e,V) 



(2.8) 



+ 



de - 

e 



f(e)-f(e + eV) 



dG s (6, V) 
dV 



with the expansion dl/dV \y= G s (— eV, 0) + 2VdyG s (e, V) \ e= _ eV y =0 +... at zero 
temperature. This differs from the differential conductance of Ref. fBlonder 1982|| , 
dl/dV = G s (—eV, 0), by accounting for the change in the conductance of the open 
channels with increasing voltage. 

The expressions ( |2.4j ) , ( |2.6| ) , and ( |2.7| ) determine the general form of the current- 



voltage relation of a disordered NS junction, for an arbitrary scattering and pair 
potential at the NS interface. 

We close this section with a discussion of the symmetry of the CVC with respect 
to sign reversal of the applied voltage. In the subgap regime e\V\ < A, the incoming 



3 This can be shown using the unitarity of the global scattering matrix including the Andreev 
process by an extension of the argument given for a normal junction. 
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quasiparticle excitations may not enter the superconductor. The probability current 
of the states |e| < A is totally reflected and thus the global reflection matrix f (e, V) 



of 



' e' ee 



is unitary. The unitarity produces the relations r\ e r ee + r\ e rh. 
J 



1 and 



+ r eh^ e h = 1- The symmetry of electron- and hole- type excitations in the 



BdG equations guarantees r e h (e, V) 
conductance takes the form 



he 



-e, V). As a consequence, the subgap 



4e 2 r ♦ 

G s (e,V) = —T T \ri(e,V)r he (e,V) 



4e V 



rjfc (e, V) r eh (e, V) 
rl(-e,^)r te (-e,y) = G.(-e,V). 



(2.9) 



A subtle issue is that this symmetry does not yet imply a symmetry in the CVC 
under reversal of voltage [[Leadbeater 1996|1 . The latter requires that G s (e, V) = 
G s (— e, — V), which amounts to G s (e,V) being independent of voltage. Then we 
have G s (e) \ e =- e v = dI/dV\v and the differential conductance is invariant un- 
der sign reversal of the voltage. Indeed, in recent experiments on SmS junctions 



Magnee 1994] , Po frier 1997a|| , an asymmetry in the CVC was found in the subgap 



regime, which can be understood on the basis of the above discussion taking into 
account the voltage dependent Schottky barrier at the SmS interface. An explicit 
account of the voltage dependence of G s requires the scattering matrix to be 
determined in the applied electrostatic potential. In principle, this task demands 
the self-consistent solution of the scattering problem and the Poisson equation, see 
Christen 1996||. In the next section, we evaluate the spectral conductance (|2.7|) by 



using the Andreev approximation for the scattering at the NS interface. 



2.1.2 Spectral conductance in the Andreev approximation 

The stationary states in the ballistic leads are solutions of the BdG equations 



de Gennes 1966a| and are of the plane wave type. A step function model for the 
pair potential A (x) = A e tx 9 (x) is assumed, which neglects the suppression of the 
pair potential in the superconductor on the scale of a coherence length. The NS 
interface connects electrons and holes of the same channel with a reflection ampli- 
tude depending on the reduced chemical potential \i v = \i — h 2 k 2 L /2m. In the limit 
e, A p, v , the BdG equations are simplified by linearizing the dispersion relation 
around the effective Fermi wave number = ^J2m^ v /h. In this approximation 
made by Andreev, incoming electrons are purely reflected into holes and vice versa. 
The reflection matrix at the NS interface is given by 
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depending on the channel independent, scalar Andreev amplitudes 



e — sign (e) \/ e 2 — A 2 



7T 

iVA 2 - e 2 



A 

2|e|' 



exp (- 



-i arccos 



I) 



|e| > A, 
lei < A. 



(2-11) 



Inserting in ( |2.4|) and (|2.7p , we obtain the multichannel spectral conductance formula 



G.(e,V) = ^(l+|r(c) | 2 )Tr^4(e) 



T*(e)r. 



22 



-e r 



22 



[i-|r 



|2 „T 



r 22 {-e)r 



22 



;-c)] [l-r(e) 2 r 22 (e)r 2 * 2 (- e )]- 1 t 21 (6)}. (2.12) 



Combined with equation (|2.6|) , Eq. ( |2.12| ) provides the finite voltage, finite tem- 
perature CVC of a disordered normal-superconducting junction in the Andreev 
approximation. The spectral conductance depends on the scattering matrices of 
the electrons at energies ±e as a signature of the presence of Andreev reflection. 
The dependence of this formula on the phases of the reflection and transmission 
amplitudes proves crucial in determining the resonance peaks in the conductance. 
The elementary process, which contributes to these phases is the propagation of an 
electron and a hole between the disorder region and the NS interface. 

If no inter-channel mixing takes place, i. e., the matrices Uj and ra are diagonal, 
the conductance G s (e, V) reduces to the quasi one-dimensional form 



2e 2 



N 

E 



(1+ | T (e) | 2 ) T v (e, V) [1- | T (e) | 2 R v (-e, V)} 



1 1+ | T (e) | 4 R v (e, V) R u (-e, V) - 2Re T (e) z r v (e, V) r* (-e, V) 



(2.13) 



r v = ir 22) vv is the reflection amplitude for a state coming in from the right side of 
the scattering region (Fig. |2~T| ) and R v = \r u \ 2 and T v = 1 — R v denote the reflection 
and transmission probabilities of the t^-th channel. 

For voltages well above the gap, |e|, V ^> A (still assuming |e| /i), the Andreev 
reflection is strongly suppressed and drops according to T (e, V) ~ A/2|e| 0. The 
spectral conductance ( [2.12|) asymptotically approaches the expression for a normal 
junction, 



G s (e,y) = ^Tr[4 (e, V) t 2 \ {e,V) 
h L 

At voltages below the gap, we make use of | Y (e) 
obtain the spectral conductance G s (e, V) in the form 



1 and T v = 1 



4e 2 



N 



T u (e,V)T u (-e,V) 



h ^ 1 + R v ( e , V) R v (-e, V) - 2Re [T (e) 2 r u (e, V) r* (-e, V)] 



(2.14) 
R v and 

(2.15) 
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The reflection and transmission coefficients at ±e are symmetrically involved in this 
formula, which results in the symmetry of the CVC discussed in the section above. 
In contrast, the spectral conductance ( |2.13| ) at voltages above the gap becomes 
increasingly asymmetric as it asymptotically approaches the Landauer expression. 
The sensitivity of ( |2.15| ) to the phase of the reflection amplitudes allows for the 
distinction of zero and finite bias peaks in the double barrier NS junction discussed 
below. 

The linear response limit (e, V — > 0) of ( |2.15| ) can be determined using T (0) 2 = 1 
and R v = 1 — T v and takes the form [[Beenakker 1992aH 



G(0) 



4^ 
~h 



E 



T v (Of 



[2-2; 



(2.16) 



2.1.3 Double barrier NINIS junction 

In this section we apply the above results to a double barrier NINIS junction, which 
allows to study the interplay between normal- and Andreev- levels in a Fabry-Perot 
type I1NI2 interferometer. First, we discuss the structure in the conductance of a 
single channel NI1NI2S junction, which we trace back to the presence of Andreev 
resonances. Second, we present numerical results for a multichannel junction, show- 
ing that the typical resonance structure of a single channel survives the summation 
over the channels. This stability is a peculiarity of the superconducting system, 
absent in a normal NI1NI2N double barrier junction. 

Since the channels separate in the double barrier problem, we can make use of 
the result Q2.13| ) for the conductance G s . G s depends on the phases ip (±e) of the 
reflection amplitudes r (±e) as well as on the complex amplitude T (e) of the Andreev 
reflection. We use the notation T (±e) = T±, R (±e) = R±, and r (±e) = y/R^e^^ 
for the reflection amplitude, the phase factors being determined by the potential 
barriers I\ and I2 and the propagation between them. We rewrite the amplitude 
of the Andreev reflection as T (e) = |r|e _i1?< - e ' ) with the phase $ (e) = arccos (e/A) 
below the gap and vanishing above. The conductance G s (e) simplifies to 
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1 + 
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*r + r_ -2|r|v^M*- 


cosb (e)-p(-e)- 20(e)] 



(2.17) 



which is always less or equal to the universal value 4e 2 /h. Note that the Andreev 
reflection is suppressed above the gap (|T| < 1), while the phase i? (e) vanishes. 
Below the gap the conductance exhibits resonances, which is maximal if R± are 
equal and the scattering phase difference of the electrons and holes compensates for 
the phase of the Andreev reflection, according to the resonance condition 



cos [<p (e) -(p(-e)-2ti (e)] = 1. 



(2.18) 
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Let us assume that the barrier at the NS interface vanishes, I2 = 0, and denote 
the length of the interferometer by d. In the limit of a high potential barrier Ii, 
R + and R- are energy independent and equal, and (p (±e) = it + 2k±d. Using the 
wave number k± = mv u ±e/v u [v v is the Fermi velocity of channel u), the resonance 
condition (J2.18[) yields the spectrum of Andreev levels, 



e n ,u = 7n (nir + arccos ^) , n = 0, 1, ... (2.19) 

which predicts resonances in the conductance of a typical width proportional to the 
transmission T of the barrier. The phase $ (e) varies between —n/2 and from 
e = to e = A and guarantees the existence of at least one Andreev resonance for 
arbitrarily small d. In the — ^ limit, the weight of this resonance lines up with 
the gap voltage and we recover the NIS junction, featuring a suppressed subgap 
conductivity and a peak in the differential conductance at the gap voltage. This 
peak can be understood in terms of the Andreev resonance which moves to the gap 
energy for d — > 0. 

We now turn to double barrier NI1NI2S junction with arbitrary barriers. Here 
ip (±e) represent the phases for the reflection of electrons entering the double barrier 
scattering region from the right (superconducting side). The corresponding reflection 
amplitudes are given by 

r ( ±£ ) =r2+ l-V 1 r 2 e^ ' (2 - 20) 

where r^t{ are the amplitudes of the left (i = 1) and the right (i = 2) barrier. The 
phase of this reflection amplitude plays the major role in determining the structure 
of the conductance, as it controls the existence of resonances according to (|2.18 ). 



Let us fix the barrier L and increase I2 slowly, keeping their strengths L > I 2 . In 
this situation, the INI interferometer develops pronounced Andreev resonances. For 
r i r 2 , the phase <p (±e) of the reflection amplitude r (±e) ~ t 2 l r 1 e 2lk±d produces a 
linear energy dependence of the phase ip (e), which changes by 2ir on the scale vpjd 
and results in (nearly) equidistant resonances, in accordance with (|2.19|) . As the 



strength of I2 is increased, the resonances pair up as is illustrated in Figure |2.2| . The 
phase function <p (e), as displayed in Figure |2]^ (solid line), can be used to determine 
the location of the resonances by finding those combinations of energies ±e, which 
have a phase difference^ Atp (e) = (p (e) — (p (— e) = 71 + 2mr. The period of Aip (e) 
reduced by half with respect to the period of cp (e) accounts for the pairing of the 
resonances. 

When the strengths of the barriers become of the same order, L ~ I 2 , the spectral 
weight of the INI interferometer is shared by Andreev quasi-bound states of a mixed 



4 The energy dependence of the Andreev reflection amplitude produces slight deviations close 
to the gap energy. 
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Figure 2.2: Andreev resonances and resonance condition for the phase (p. Bottom: 
phase (p (e) of the reflection amplitude versus energy. The solid line represents the 
I1NI2 interferometer with barriers strengths R\ = 0.8 and R2 = 0.5, the dashed 
line stands for the inverse barrier sequence (Ri = 0.5, R2 = 0.8). The Andreev 
resonance condition for the phase is met for a pair of energies ±e ra with phase 
difference Atp (e n ) = ti + 2mr. This phase condition can be fulfilled only by the 
first barrier sequence (Ri > R2, solid line) at those energies indicated in the graph. 
The dotted line shows the transmission probability T as determined by the normal 
resonances. Top: conductance (arbitrary units) of the double barrier NS junction 
versus energy, the solid line again representing the barriers Ri = 0.8, R 2 = 0.5, 
and the dashed line the barriers Ri = 0.5, R2 = 0.8. Note the symmetry of the 
resonances with respect to e = 0. 



electron-hole character and normal electron quasi-bound states. Due to the large 
gradient of the phase close to the normal resonances, see Fig. [2.2| , the Andreev 
resonances tend to be pinned to normal resonances at either +e or — e. While the 
Andreev bound states contribute to the current, normal bound states do not couple 
to the superconductor and thus do not participate in the charge transport. This is 
reflected by the subgap symmetry of the G s (e) under reversal of voltage which is is 
observed for all barrier strengths, see Fig. [2.2| . 

As the barrier strength is increased further to I2 > Ii, the Andreev resonances 
are weakened and eventually disappear. Although the normal resonances dominate 
in the INI interferometer in this regime, only the weak Andreev resonances show up 
in the conductance, which thus exhibits only a weak subgap structure. The phase 
function ip (e) in Figure |2.2| (dashed line) becomes nearly constant for 7*1 <C r 2 , see 
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( |2.20D , and the phase condition for resonance ( |2.18| ) cannot be met. 



A comparison of the double barrier systems NI1NI2S and NI2NI1S, i.e., with 



inverse sequences of the barriers Ii and I2, is given if Fig. |2.2| . The transmission 
T (e) (dotted line) is identical for both cases. Let us assume that I x I 2 . In 
this first barrier sequence, we have a strong energy dependence of the phase (p (e) 
(solid line in Fig. |2.2j ), which implies the existence of Andreev type resonances at 
finite bias. The electrons entering the INI interferometer from the normal lead, are 
given enough time to build up an Andreev resonance and preferably leave into the 
superconductor. For the inverse barrier sequence, the barrier I2 at the NS interface 
dominates. The weak energy dependence of the phase (p (e) of reflection allows 
no sharp resonances to build up (dashed line in Fig. |2.2| ). This reflects the fact 
that the electrons which enter the INI region leave through I x into the normal lead 
before an Andreev resonance can build up. In summary, the spectral density in 
the INI interferometer changes radically with the coupling strengths of the normal 
and superconducting leads. Normal resonances dominate when the interferometer 
is coupled more strongly to the normal lead, whereas Andreev resonances take over 
in weight when the coupling to the superconductor is stronger. At any instant, 
however, only the Andreev states participate in the charge transport. 

We turn to the numerical analysis of a multichannel NINIS junction. With Eq. 
( [2.13|) we extend the linear response study of Ref. ||Melsen 1994j| to finite voltage. We 



investigate an NI1NI2S junction with two 5-function barriers of typical strength H = 
J V(x)dx ~ hvp and corresponding reflection probability R = H 2 / (H 2 + h 2 v 2 F ) as- 
suming values between R = 0.2 and R — 1. We vary the relative barrier strengths 
to cover the range between the two limits Ii > I 2 and I x < I 2 discussed above. The 
distance between the barriers is chosen to be of the order of or larger than the coher- 
ence length of the superconductor, providing the forward channel with one to a few 
Andreev resonances. The number of resonances increases with the incidence angle of 
the channels. We chose leads with a cross section area of (lOO/fc^) 2 , which amounts 
to about 800 transverse channels. The ratio of the energy gap to the Fermi energy 
is assumed to be A/ep = 0.002. Each channel features the typical conductance 
structure of paired Andreev resonances exposed above. Their positions and widths 
depend on the ratio of the barrier strengths Ii and I 2 as well as the longitudinal 
kinetic energy of the single channels. Remarkably, the overall conductance, which is 
obtained through summation of single channel conductances, exhibits a characteris- 
tic subgap structure signaling the presence of Andreev resonances. In contrast, the 
overall conductance of the corresponding N^N^N junction is practically constant, 
the normal resonances of the INI region having been averaged out. 



The numerical study of a 3-dimensional NINS junction has shown that both 
the positions and the number of resonances in the overall conductance correspond 
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Figure 2.3: Differential conductance for a multichannel NINIS junction with in- 
terbarrier distance d = 2vp/A = 2ttC,. The average conductance per channel is 
plotted versus voltage at temperature T = 0. The barrier strengths are chosen to be 
Ri = 0.2 (dotted line), 0.5 (dashed), 0.7 (dot-dashed), and 0.8 (solid) in the normal 
region, while R2 = 0.5 at the NS interface. The normal state conductance, which is 
roughly independent over this voltage range, is indicated on the left. With increas- 
ing barrier strength Ri the zero bias anomaly develops into a finite bias anomaly as 
the Andreev resonance is formed for R\ > R%. 



to those in the forward channel [jChaudhuri 1995]! 5 !. In the NINIS junction, we do 
not find a direct correspondence of the resonances of the total conductance with 
the forward channel nor with any other specific channel, although a clear resonance 
structure still survives the summation over the channels. 

Let us concentrate on the conductance formula ( j2.15| ), valid at subgap voltages, 
and on the properties of the CVC close to zero voltage. For I x > I 2 , the denominator 
of ( ^.15| ) changes rapidly with the strong energy dependence of the phase of r v (e) 
which is responsible for the appearance of conductance peaks at finite voltage. The 
pronounced structure in the conductance survives the summation over the channels 
as displayed in Figures |2.3| and [2.4| (solid lines). The repulsion of the Andreev levels 



3 This is connected to the decrease of the transmission probability with increasing incidence 



angle, and the non-isotropic distribution of the channels over the incidence angles [Chaudhuri 1995 1 
which is oc sin (2z9). 
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around zero voltage produces a minimum in dl/dV at zero voltage. For Ii < I2, 
the phase of the reflection amplitudes r v (e) has negligible energy dependence and 
the numerator of ( |2.15| ) dictates the features of the conductance. The expansion of 
the product T v (e) T v (— e) = T 2 — e 2 T^ 2 about zero energy shows the existence of a 
zero bias maximum. The zero bias anomaly shows up as a characteristic property 
of the overall conductance, see Figures [2.3| and |2.4j (dotted lines). The zero bias 
maximum coincides with the maximum of the conductance product G (e) G (— e) of 
the corresponding NININ junction at zero energy. Figs. |0| and |2]^ illustrate the 
crossover from zero to finite bias anomalies for two different interbarrier distances 
d as the strength of barrier Ii is increased and I2 is kept fixed. For an interbarrier 
distance larger than the coherence length of the superconductor, several Andreev 
resonances show up (see Fig. [2.4| ). Note that the inversion of the barrier sequence 
transforms the zero voltage conductance from a local minimum to a local maximum, 
while keeping the same zero voltage conductance. 
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Figure 2.4: Differential conductance for a multichannel NINIS junction of width 
d = 4v p/ A = 47r£. The average conductance per channel is plotted versus voltage for 
temperature T = 0. The barrier strengths are Ri = 0.04 (dotted line), 0.2 (dashed), 
0.4 (dot-dashed), and 0.54 (solid) in the normal region, and fixed at R2 = 0.2 at the 
NS interface. The normal state conductivity is indicated on the left. As the strength 
Ri is increased, the zero bias anomaly turn into an finite bias anomaly and several 
Andreev resonances appear. For the last choice of barriers we have interchanged the 
barrier sequence: the conductance at zero voltage remains the same while changing 
from a zero bias minimum to a zero bias maximum (short solid line). 
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The interference of multiple scattering processes between the scattering region 
and the NS interface thus produces an interesting structure in the differential conduc- 
tance. Recently, experiments in disordered NS and SmS junctions [[Kastalsky 1991 



Nguyen 1992] , |Nitta 1994] , Pakker 1994] , |Magnee 19941 have concentrated on the ob- 
servation of zero and finite bias anomalies. It has been understood theoretically 



Volkov 1993| , |Marmokos 1993| , |Yip 1995] , [van Wees 1992|| that these features are due 



to the interplay between the barrier at the NS interface and disorder in the normal 
lead. At small disorder, the differential conductance exhibits a zero bias maximum, 
while at large enough disorder, a finite bias peak is expected |[Yip 1995| at a voltage 
of the order of the Thouless energy E c of the normal lead [[Edwards 1972 ]. This 
has recently been confirmed in an experiment [Poirier 19974 Poirier 1997b|| . Here, 
we have found the existence of the analog features in a ballistic double barrier NS 
junction (in addition, the double barrier junction shows higher harmonics in the res- 
onances of the conductance). The ballistic point of view applied to the disordered 
NS junction thus explains the finite bias anomaly as a superposition of resonances 
due to quasi-bound Andreev states between the superconductor and the disorder. 
This interpretation applies in the same way to recent experiments showing a low- 
temperature reentrance in the conductance IPetrashov 1995| , |Charlat 1996|| . The 
non-monotonicity is again due to the resonance structure in the spectral conduc- 
tance G(e), which now appears as a function of temperature. 



2.1.4 Discussion 

We have described the current- volt age characteristics of NS junctions by the scat- 
tering matrix approach and expressed it through the spectral conductance, which 
takes into account the explicit voltage dependence of the scattering problem. The 
spectral conductance is symmetric under sign reversal of voltage in the in the subgap 
regime. We have derived Eq. ( |2.13| ) giving the spectral conductance in terms of the 
normal scattering amplitudes. This result has enabled us to carry out a study of a 
double barrier NINIS junction at finite voltage. As the ratio of the barrier strengths 
of the INI interferometer is changed, and we go over from normal to Andreev reso- 
nances, the conductance shows a cross-over from a zero bias to a finite bias peak. 
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2.2 Finite voltage shot noise in normal- metal — 
superconductor junctions 

2.2.1 Introduction 

Shot noise arises from the current fluctuations in transport as a consequence of 
the discrete nature of the carriers, first predicted by Schottky for the vacuum tube 
Schottky 1918| . In a coherent conductor the quantum fluctuations, which follow 



from the probabilistic nature of the backscattering restricted by the Pauli princi- 
ple, are reduced in comparison to the Schottky result [[Lesovik 1989| , Piittiker 1992 ] . 
In a disordered normal-metal-superconductor (NS) junction, the shot noise is pro- 
duced by the normal scattering processes as well as the imperfect Andreev reflection 
Khlus 1987| , |lVluzykantskii 1994] , |de Jong 1994] , |Hessling 199lj . The Andreev reflec- 



tion introduces fluctuations proportional to the double electron charge, which may 
be iteratively increased to give fluctuations of several charge quanta in biased SNS 
junctions [[Averin 1996|| . For a review, we refer the reader to ||de Jong 1997|| . In this 
chapter we study the shot noise at finite voltage in a disordered NS junction as well 
as its structure due to the iterative scattering processes between the disordered nor- 
mal lead and the NS interface. We derive a formula which expresses the differential 
shot noise of the dirty NS junction in terms of the scattering matrix of the normal 
lead and the (scalar) amplitude of Andreev reflection. In the systems of a single 
barrier normal-metal-insulator-normal-metal-superconductor (NINS) junction and 
a double barrier NINIS junction, we explain the existence of a non-trivial resonance 
structure. The results of chapter have been published in [ [b'auchere 1998| |. 

We consider a disordered NXS junction, shown in the inset of Fig. 2~ 5, with an 
arbitrary elastic scattering region X in the normal lead, whose effect on the noise 
is to be determined. We restrict ourselves to voltages below the superconducting 
gap eV <C A. The low-frequency power spectrum of the current fluctuations is 
determined by the irreducible current-current correlator 



P 



dte tuJt ((I(t)I(0))), u^O 



(2.21) 



(((I (t) I (0))) = ((I(t) - {I}) (1(0) - (I))) is the second cumulant in time). The 
time-dependence of the current operator is defined by 



/ (t) = exp [i (H - fiN) t] I exp [-i (H - fiN) t] 



(2.22) 



In the mean field approximation for the Hamiltonian H we do not account for the 
fluctuations of the order parameter in the S region. The effective Hamiltonian is 
diagonalized by a Bogoliubov transformation fiNifio 198"4]]. After integration over the 
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Figure 2.5: Differential conductance (top) and shot noise (bottom) in a multichan- 
nel NINIS junction (see inset, symmetrical barriers of strength J dxV (x) = 3hvp, 
T\ = T2 ~ 0.05, interbarrier distance d = 20vp/A, 8 x 10 4 channels). The average 
shot noise per channel [Eq. (|2.28|) 1 exhibits maxima at the resonances of the con- 
ductance (solid lines). Note the enhanced structure of the noise with respect to the 
conductance [Eq. (|2.26|) 1. The dashed lines show the corresponding results for the 
normal state NININ junction. Inset: schematic NXS junction, e.g., X = INI. 



cross-section, the net current operator can be expressed through the positive energy 
eigenfunctions, 



I (*) = —? Y] / dydz(u* m d x Urala n - v* m d x v n ^ n ^e l{er> 
mi J 



[2.23) 



em,£n>0 



which may be evaluated in the normal lead. The operators 7 m belonging to the 
wavefunctions (u m , v m ) annihilate the scattering states of the disordered NS junc- 
tion (note ud x v = ud x v — vd x u). We note that the symmetry of the BdG eigen- 
functions (u n ,v n ), e n with respect to sign reversal of energy to (— t>*,-u*), — e n is the 
consequence of spin degeneracy in Nambu space, allowing the the negative energy 
states to be eliminated in favor of the positive energy states, see Appendix At 
voltages below the superconducting gap (eV < A), the quasi-particles are injected 
from the normal reservoir only, and the wavefunctions depend on the N x N global 
reflection matrices r ee (e), The (e), r e h (e), and Thh (e) of the disordered NS junction, 
see Eq. ( |2.3| ). The scattering states indexed by m = fi, e) consist of an incident 
e(lectron) or h(o\e) like quasi-particle in channel fi with energy e superposed with 
the reflected electron and hole states. The occupation numbers are given by the 
Fermi-Dirac distribution f e h = f (e ± eV). 

The noise power ( j2.21j) is determined by the transitions induced by the cur- 
rent operator / from initial states |n) = \n n = l,n m = 0) to intermediate states 
\ m ) — \ n n — 0, n m = 1) and back, which differ only by their occupation with respect 
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to two single particle states with indices n and m. E.g., the transition between an 
incident electron n = [e, u, e] and an incident hole state m = [h, fi, e] is produced 
by the matrix element between the reflected electrons and holes of the two scat- 
tering states with respect to the current operator, (m\I\n) oc f e (1 — fh) {r\ h r ee — 



r hh r he 



' [IV 



Summed over the channels, the transitions contribute to the fluctua- 



eh r ee 



tions with the weight Y i ^\( m \ I \ n )\ 2 ^ fe (1 - fh) Tr{ {r\ e r eh - r[ e r hh )(r 

r hh r he)} = / e (l - fh)Tr{rl e rhe(l - r{e r he)}- Following the quantum mechani- 
cal formalism outlined, we obtain the low-frequency limit of the power spectrum 



Anantram 19961, IMartin 19961 valid for T, eV < A 



o 



[/e(l-A) + A(l-/e)]TV 



+ \f e (l-f e )+f h (l-f h )]Tr 



r{ e r he ( 1 - r\r h 



ie 1 



(2.24) 



The first term describes the transitions between states of opposite current sign, while 
the second represents the transitions between states of equal current sign. At zero 
temperature, the second term is suppressed by the Pauli principle and the first term 
produces the shot noise [[Lesovik 1989| , |Lesovik 1994| , [Levitov 1996|| . The thermal 
fluctuations describing the Johnson- Nyquist noise ||Johnson 1927| , [Nyquist 1928|| are 
due to both terms. Eq. Q2.24 ) is manifestly invariant under sign reversal of voltage, 
see Chp. |1[ 



2.2.2 Interplay between normal scattering and Andreev re- 
flection 



We express the electron-hole reflection matrix r/j e (e) in Eq. (|2.24|) in terms of the 



normal reflection- and transmission matrices ru and of the scattering (X) region 
and the amplitude of Andreev reflection T at the XS interface, 

r he (e) = t\ 2 (- e ) [1 - T 2 ( e ) r 22 (e) r* 22 (-e)] ~' T (e) t 21 (-e) , (2.25) 

see Eqs. ( (2.4Q and ( |2.10| ), recalling T (e) = e 1 ®^, (e) = arccos (e/A)). Inserting 
The (e) in Eq. ( [2.24|) provides us with a general multichannel expression for the shot 
noise in a disordered NXS junction, explicit in the normal scattering matrix and the 
Andreev reflection amplitude. 

Next we restrict ourselves to the situation of a junction with uniform transverse 
cross-section, which permits the separation of channels and will be used in the single 
and double barrier systems below. We recall that the differential conductances of 
the channels v = 1...N are given by 

h ;Gu(e)=(rlr he ) , (2.26) 



4e 2 



in.' 
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see Eq. ( |2.15| ). At zero temperature, the shot noise fl2.24| ) exhibits the noise power 



P = - deVS„(e), (2.27) 
Jo 



e 



with the differential (low-frequency) shot noise of channel u, 

—S v (e) = 4 (rl e r he (l - r\ ie r hf ^j ^ 

AT U (e) T v (-e) (_K (e) + (-e) - 2Re [r (e) 2 (e) (-e)] ) 
(T„ (e) T„ (-e) + ^ (e) + ^ (-e) - 2Re [r (e) 2 r v (e) r* (-e)] )' 



(2.28) 



r u (e) is the reflection amplitude of channel v in the matrix r22 (e) describing the 
reflection from the right-hand side of the X region, R u (e) = \r„(e) | 2 , and T u (e) = 
1 — R u (e). The two energy dependencies of the reflection probabilities R u (e) and 
(the phases of) the reflection amplitudes r v (e) translate into the voltage dependence 
of the shot noise. The term sensitive to the phases of the reflection amplitudes is 
due to the states multiply scattered between the disorder and the NS interface. Note 
the similarity of Eq. ( |2.28| ) to the conductance formula ( |2.13 ). In the limit e — ► 



— > — i) the phase dependencies drop out of Eq. (|2.28|) and we recover the linear 
response result ||de Jong 1997|| , 



In order to describe the weak coupling of a normal lead to an NS proximity sand- 
wich, we consider a normal-metal-insulator-normal-metal-superconducting (NINS) 
junction with a potential barrier of low transmission (T 1) placed at a distance d 
away from the perfect NS interface. The NINS junction serves as a model system for 
a tunneling experiment from a metallic tip to a thin film NS layer structure, which 
has permitted the observation of the Rowell-McMillan oscillations [[Rowell 1966|| . In 
this elementary model, the reflection amplitudes r v (e) = \fRe llfv ^ have a roughly 
constant modulus R and an energy dependent phase ip u (e) w 2 (k u + e/v v ) d + <p , 
accumulated during the propagation between the potential barrier and the NS inter- 
face (k u , v v denote the Fermi wave number and velocity of channel v). Consequently 
the shot noise ( [2. 28] ) for the channel u, 



h . 8T 2 R[1 - cosa u (e)] 
2e3 uK ! ~ (T 2 + 2i?[l - cosq;^ (e)]) 



Su (e) = ,S, ""/Li , (2-30) 



depends only on the phase difference a v (e) = <p v (e) — f u (— e) — 2$ (e) of electron 
and hole propagation through the X (=IN) region. The resonance structure of S v is 
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intimately connected to the voltage (energy e = eV) dependence of the differential 
conductance ( |2.15|) , 



h T 2 

4^° u (C) = T 2 + 2R[l-cosa u (e)] ' (2 ' 31) 

Figures |2.6| (a) and |2l](b) show their generic dependence on a v . The minima of the 
denominator [coso^ = 1 =>- e v<n = v u /2d{nn + arccose^ n /A)] correspond to the 
resonances of the Andreev quasi-bound states. These conductance resonances are 
repelled from zero voltage, since the phase of the Andreev reflection d (e — > 0) = tt/2 
has to be compensated by the phase difference of electron and hole propagation 
<f v (e) — ip u (— e). The shot noise (|2.30|) vanishes at these resonances (S u (e v ,n) = 0). 



The noise maxima are found from cos a u = (2R — T 2 ) /2R at energies doubly peaked 
close to the resonances. The peak separation 5e ~ Tvp/d coincides with the width 
of the resonances. 

Interestingly, this non-trivial structure survives in the multichannel NINS junc- 
tion, whose results are shown in Figs. |2.6| (d) and |2.6| (e). The stability of the reso- 
nance structure is due to its pinning to the Fermi energy. In contrast, the resonance 
structure of a potential well in a normal conductor is washed out in the multichannel 
limit. As Figures [2.6| (b,e) show, the narrow double peak structure in the noise is 
smeared out in the multichannel limit, and the noise S takes a maximum rather 
than a minimum at the resonance, retracing the shape of the conductance. 

At large voltages (eV ^> vp/d), the shot noise approaches a constant as a con- 
sequence of the dephasing between the channels. The channel average may be eval- 
uated by averaging over the phase a u in Eqs. fl2.30|) and (|2.31|) and we obtain the 
result (a„ (l/iV£„ -> l/2vr ff* da), 



1 r 2w 2e 3 - 2e 2 

- daS u (a) = — T, G = —T. (2.32) 
In J h h 



-27T 



It follows that both the shot noise and the conductance approach the normal state 
values at voltages vpjd eV A. This limiting behavior demonstrates that the 
NIN junction is effectively decoupled from the NS interface in the large voltage limit 
and dominates both noise and conductance due to its low transparency. 

Let us concentrate on the overall ratio of noise power to current P/el [P = 
J^ V ^deS(e), I = Jq V de G(e)]. This ratio has already been used successfully in 
determining the unit of the charge carriers in the Fractional Quantum Hall Effect 
Glattli 1997| , Ide-Picciotto 19971| . While at small voltages (P/el = 2) the noise car- 



ries the signature of the Cooper pairs created by the Andreev reflection, at large 
voltages {yp/d <C eV A) it decays to the normal state value (P/el = 1). In- 
terestingly, this decay is immediate at the first Andreev resonance, as is seen from 
Figs. 2(c) and (f). In fact, in a single channel, the energy integration is equivalent to 
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Figure 2.6: Zero temperature conductance (a,d), differential shot noise (b,e), and 
finite voltage noise power to current ratio (c,f ) for a NINS junction with a barrier of 
strength f dxV(x) = 3Hvf, average transmission T = 0.05, e F = 500A, d = 20t>i?/A 
(vp/d < A < e F ). (a,b,c) follow from Eqs. (|2.30|) and (|2.31|) for one channel. (d,e,f) 
from Eqs. ( [2.15 ) and (2.28) are averaged over 8 x 10 4 channels. Note the pronounced 
resonance structure in both conductance (d) and noise (e), which approaches the 
values for the corresponding NIN junction at large voltages (dashed lines). The 
noise to current ratio (c,f) decays to the classical Schottky value P/el = 1 above 
the Andreev minigap. 



the integration over the phase a u , which at the midpoint between two resonances a u 
has covered the period 2n, thus producing P/el = 1 at a v = ir, see Fig. 2(c). In the 
transport through a normal lead weakly coupled to a NS sandwich, we thus find two 
characteristic regimes: At voltages below eV < vpjd the one-particle excitations in 
the NS structure are effectively gapped and we may only couple to the supercon- 
ducting condensate which introduces the charge 2e in the noise power. At voltages 
above the first Andreev resonance, we may couple directly to a one-particle density 
of states in the normal layer and we loose the signature of the superconductor at 
once. 



2.2.3 Double barrier junction 

In the study of normal versus superconducting junctions, it is of interest to investi- 
gate the impact the quantum coherence has on the macroscopic limit of transport 
(limit of an infinite number of channels) and, in particular, to compare the behavior 
of a specific scattering region connected to either a normal or a superconducting 
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reservoir. We consider an insulator-normal-metal-insulator (X=INI) interferometer 
which we connect to normal leads in a NININ junction and to a superconducting 
lead in a NINIS junction. Besides the intrinsic interest in the double barrier system 
as a paradigm, now to be combined with a superconductor, it may also serve as 
a qualitative model for a disordered conductor, due to its strong similarity in the 
transmission distribution fDorokhov 1982j , [Melsen 1994|| . As we consider voltages 
below the superconducting gap, we typically chose the thickness of the interlayer to 
be d » hv F /A. 

The resonance structure of the INI interferometer results in a bimodal distribu- 
tion of its (overall) transmission probabilities T. For symmetric barriers it takes the 
form (infinite channel limit, 7\ = T2 -C 1) ||de Jong 199711 , 



P(T) 
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vr 2 ^T 3 (1 - T) ' 



T e 



Tl 



71- 



(2.33) 



which has its analog in a disordered conductor [Porokhov 1982| , [Melsen 1994|| . We re- 
call that, in a NININ junction, the rich structure of this distribution has no impact on 
the macroscopic transport properties. The (linear) conductance and shot noise may 
be evaluated from the bimodal distribution (|2.33|) , (h/2e 2 ) G = J dT p (T) T = Xi/2, 
(h/2e 3 ) S = J dTp(T)T(l — T) = Xi/4, and show no characteristics of coher- 
ent transport [|Chen 1991|| . The conductance follows from the series resistance 
of the two barriers, and an incoherent model of a double barrier junction gives 
the same suppressed noise S/eG = 1/2 as a consequence of charge conservation 
Bcenakker 1992 b| . When considering the INI interferometer in a NINIS junction 



instead, we find a non-trivial answer, both at small and large voltages. The linear 
response follows from Eqs. (|2~T5| ), and (|2~33l) , 
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The ratio of shot noise to conductance is ||de Jong 1997|| S(0)/eG(0) = 3/4. A 
comparison of the shot noise to conductance ratios in a single barrier junction (X=I 
or X=IN) and a junction with disorder (X=D) is instructive, see Table 1. In the 
single barrier junction the distribution of transmissions is peaked at T C 1. The 
noise is due to the Schottky type fluctuations in nearly closed channels, and con- 
sequently the noise ratios Sjy/eG^ = 1 in the NIN and Ss/eGs = 2 in the NIS 
junction differ only by the charge quanta involved. In a junction with disorder, the 
noise-to-conductance ratio is also doubled for the NDS case, see Table 1. However, 
since in the bimodal distribution of the transmissions the noise is produced by the 
channels with intermediate transmission < T < 1 and the current by the open 
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channels with T — > 1, the doubling is a non-generic and thus yet unexplained fea- 
ture. This is demonstrated by the noise-to-conductance ratio of 3/4 in the NINIS 
as compared to 1/2 in the NININ junction. 



Table 2.1: Shot noise to conductance ratio S/eG of a NXN compared to a NXS 
junction; single barrier X = I, double barrier X = INI, disorder X = D. The results 
are valid for small transmissions T< 1, and many channels (N -> oo) ffKhlus 19871 , 
de Jong 1997] , |Chen 1991] , |Beenakker 1992g]. 





single 


double barrier 


disorder 




barrier T 1 


Ti = T 2 < 1 






1 


1/2 


1/3 




2 


3/4 


2/3 



At finite voltage, the shot noise is described by Eq. (|2.28 ) with the reflection 



amplitudes r u (±e) = r 2 + i^ne 2 ^ / (l ~~ n^e 21 *^* 6 ^) of the double barrier 
interferometer. We have evaluated this expression for a symmetric double barrier 
junction (T x = T 2 = 0.05) and display the results in Fig. |2"75| . At voltages of the 
order of the Andreev levels eV ~ vp/d, we find a resonance structure independent 
of the number of channels. We observe again that the differential noise follows the 
resonance peaks of the conductance. At large voltages (eV 3> vp/d), the resonances 
disappear as the electrons and holes dephase, and we approach the regime where the 
conductance and shot noise become indistinguishable from the incoherent addition 
of the NIN and NIS junctions. 



2.2.4 Discussion 

In conclusion, we have expressed the differential shot noise in a disordered NXS 
junction in Eq. ( |2.28[ ) in terms of the normal and Andreev scattering amplitudes. 
We have described the resonance structure in the shot noise found at finite voltage in 
coherent transport. The robustness of the resonance structure in the multichannel 
limit is owed to its pinning to the Fermi energy. We have pointed out the possibility 
of a non-trivial normal versus superconducting noise ratio in the NINIS junction as 
as consequence of the bimodal distribution. Finally, we have found a rapid decay of 
double shot noise in a NINS junction above the Andreev minigap. 
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Chapter 3 

Magnetic response of a 
normal-metal — superconductor 
structure: Nonlocality 



3.1 Introduction 

In this chapter we turn to the magnetic response of a normal metal in proximity to 
a superconductor. While the transport experiments discussed in Chapter ^| probe 
the Andreev quasi-particles in an out-of-equilibrium situation, which is kept up by 
the voltage difference across the leads, the magnetic screening is a thermodynamic 
property of the condensate wavefunction. Induced superconductivity, in fact, does 
not necessarily imply supercurrents. Though the Andreev reflection may enhance 
the conductance in a normal- metal - superconductor junction, the transport cur- 
rents remain nevertheless dissipative. In the response to an applied magnetic field 
considered here, however, the Andreev quasi-particles produce diamagnetic super- 
currents. The magnetism is a property of the H — T phase diagram of the proximity 
system. 

Here we study the linear current response to an applied magnetic field. Com- 
bined with the self-consistent solution of the Maxwell equation this allows to de- 
scribe the diamagnetic phase of the H — T phase diagram. To be specific, we 
consider an normal-metal slab of thickness d in contact with a bulk superconductor, 



as shown in Fig. 3.1. The pioneering works in the magnetic response of hybrid NS 



structures go back to the Orsay Group ||de Gennes 1964] , |de Gennes 1966bj| . Using 



the Ginsburg-Landau (GL) equations for superconductivity, see ||de Gennes 1966b| . 
Geshkenbein 1988| , they gained a good understanding of the rich phenomenology 
of these mostly dirty thin film systems. The normal metal in proximity showed 
an increasingly diamagnetic response with decreasing temperature, eventually ap- 
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proaching a perfect Meissner-Ochsenfeld effect [[Meissner f933|| . Although they ob- 
tained the correct qualitative picture, the GL equations were pushed beyond their 
strict domain of validity. The progress in the fabrication technology opened the 
way for interesting experiments on increasingly clean NS proximity systems in the 
last two decades |pda Tggg , |Mota Bgg, |Bergmann Tgg7] , |Mota Egg |Visani 1990a 



Mota 1994] , |Onoe l995f . Their quantitative understanding requires a more accurate 



microscopic description. Most importantly for this thesis was an astonishing finding 
made by Mota and co-workers ||Visani 1990a|| in the low-temperature behavior of 
normal-metal coated superconducting cylinders. These Ag-Nb cylinders were shown 
to exhibit a paramagnetic signal on top of the diamagnetic susceptibility. Before a 
closer analysis of this fascinating experiment became possible — we will return to it 
in Chapter |5| — a quantitative analysis of their experimental data at higher temper- 
ature was needed for a characterization of the sample with respect to the proximity 
contact and the impurity concentration. In fact, the theory by ||de Gennes 1966b|| 
could not give an quantitative description of the usual diamagnetic behavior. The 
work described in this chapter has laid the basis for this understanding, tracing the 
screening behavior back to the non-locality of the current-field relations shown to 
be typical for the proximity effect. 

A mile-stone in the theoretical development was the work by Zaikin [ [Zaikin 1982 
on the magnetic response of a clean normal-metal proximity layer. Applying the 
quasi-classical Green's function technique to this problem, he demonstrated the 
power of this approach, which has been predominantly used since | Narikiyo 1989 , 
Higashitani T995j Belzig 1996aj 



Section £372] is devoted to the quasi-classical tech- 
nique in the proximity effect. Zaikin found a non-local current-field dependence in 
this system, which exhibits a similarity to the Pippard equation fPippard 1953|| . In 
section |3.3| we derive the general linear response kernel for a NS sandwich with an 
arbitrary impurity concentration. We find that the clean limit represents just one 
facet of the typical nonlocal constitutive relation in the proximity effect. The rich 
phenomenology is due to the self-consistency with the Maxwell equations required 
by the screening problem. In section |3.4| , we study the sensitivity of the nonlocality 
to the presence of impurities. The superfluid density and the range of the current- 
field relations, both diminishing with decreasing mean free path, are shown to affect 
the screening ability in a contrary way and compete in the magnetic response. 



3.2 Quasi-classical equations and proximity effect 
3.2.1 Eilenberger equation 

The basic set of equations appropriate for describing spatially inhomogeneous su- 
perconductors was developed by Eilenberger [|Eilenberger 1968"| and by Larkin and 
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Ovchinnikov [[Larkin 1968|| . Starting from the widely used field theoretic method for 



superconductivity, see ||Abrikosov 1975| , |Rickayzen 1980| , |Mahan 1990|| , by averaging 



over the oscillations on the scale of the Fermi wavelength, a quasi-classical Green's 
function is obtained, which retains the full information on the length scales of inter- 
est: the coherence length and the field penetration depth. The 2x2 matrix Green's 
function g = <L n (x, v F ), 



9 f 
P '9 



(3.1) 



satisfies the Eilenberger equation^] 

- (vf • V) g Wn (x, v F ) = [(u n + iev F ■ A (x)) f 3 + a Wn (x), g Wn (x, v F )] , (3.2) 

with the normalization condition g 2 = 1. Here g = g Wn {^-, Vp) is the normal- and / = 
/aj n (x, Vp), p = /] (x,Vf) are the anomalous Green's functions, depending on the 
Matsubara frequency u n = (2n + 1)ttT (n = 0, ±1, ..), the center-of-mass coordinate 
x, and the direction of the Fermi velocity vp (t, denote the Pauli matrices, [-, ■] is the 
commutator). The Eilenberger equation allows to describe spatially inhomogeneous 
superconductors in equilibrium. The self-energy a = + <jW + ... contains the 
off-diagonal pair potential 

- ( A,x)- V ) ' <«) 

as well as other contributions from impurity scattering or higher order electron- 
electron interactions. The quasi-classical technique represents a self-contained for- 
malism, i.e, the self-energies can be expressed through the quasi-classical Green's 
functions themselves. Below we use the impurity scattering self-energy in the Born 
limit, 

^S(x) = ^(^(x,v F )>, (3.4) 

where r the scattering time and (...) denotes the angular average J dQ/4n(...) over 
v>. The order parameter is self-consistently determined by (N = mkp/h 2 TT 2 ) 

A(x) = —VNqtxT (/-n(x, v F )>, (3.5) 

where V < describes the effective coupling constant for the attractive electron- 
electron interaction. The systems of Eqs. ( |3.2| ) and ( |3.5| ) are completed by the 
quasi-classical current expression, 

j(x) = 2ieN nT (v F 9. n (x, v F )>. (3.6) 

w n >0 



1 e = |e|, h = Ub = c = 1 throughout this thesis. 
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In the Gorkov formalism [[Abrikosov 1975|| the current is separated into a diamag- 
netic and a paramagnetic contribution. Eq. (|3.6| ) contains both contributions, even 
if it is formally similar to the paramagnetic current in terms of the Gorkov-type 
Green's function. This is due to the renormalization procedure used in the quasi- 
classical theory and is explained in appendix |TJ. 

The Eilenberger equation is valid on all length scales ( with kpC ^ 1 ; an d 
for magnetic fields such that the Larmor radius = hkpc/eH exceeds the involved 
length scale tl ^> ( {tl being the cyclotron radius of a particle with Fermi velocity). 
The quasi-classical Green's functions need to be supplemented by a set of boundary 
conditions | Zaitsev 1984 , Kicsclmann 1987 , Hara 1993|j , because they are not con- 
tinuous at the interfaces of different materials. For more details on the derivation, 
the symmetries and the validity of the Eilenberger equations, we refer the reader 
to appendix || and several reviews on the subject ||Serene 1983| , |Shelankov 1985| . 



Rammer 1986| , |Kopnin 1997 



Let us illustrate the use of the quasi-classical equations for the simple case of a 
homogeneous superconductor. The Green's function takes the form 



1 



u n A 

A — 0J r . 



(3.7) 



where f2 n = \fu^ + A 2 (oj n > 0). Eq. (|3.5|) yields the usual gap equation for the 



superconductor. Within a field A(x) the solution of ( |3.2| ) is found by substituting 
u n — ■> uj n + ievp ■ A(x) in ( |3.7| ), assuming a slow variation of the field on the 
scale A £o larger than the superconducting coherence length £ = hvp/2A. The 
current-field relation that follows from Eq. (|37 



j(x) 



mc 



A(x), 



is the local London relation with the superconducting density 
n = -^vrT V — 

o;„>0 v "■' 



(3. 



(3.9) 



approaching the full electron density n = kp/37T 2 for T — > 0. In the limit of vanishing 
superconductivity A — >• we recover the trivial Green's function for the normal 
metal g = t%, and the current expression vanishes. In the following section we apply 
the quasi-classical formalism to the proximity effect in a normal metal slab. 



3.2.2 Clean NS slab 



In this chapter, we consider a normal-metal slab in contact with a bulk super- 
conductor, see Fig. |3.1| . The normal layer (0 < x < d) is described within the 
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x 



Figure 3.1: Normal- metal layer (thickness d) in proximity to a bulk superconductor. 
A magnetic field H parallel to the surface is applied, driving screening current j 
along the surface. Current and fields are inhomogeneous in x-direction. 



free electron gas approximation. We assume a perfect contact between the metal 
and the superconductor as well as a specularly reflecting metal boundary. A mag- 
netic field (0, 0, H) is applied parallel to the metal surface, driving screening cur- 
rents (0, j y (x), 0) along the surface. The magnetic induction is described by the 
vector potential (0, A y (x), 0). Taking a thickness d ^> £o> we can neglect the self- 
consistency of the pair potential, and approximate the order parameter by a step 
function A(x) = AQ(—x). 

The relevant length scales in the system are given by the geometric thickness d, 
the thermal coherence length £n{T) = Hvf/^ttT ', and the elastic mean free path I = 
vpT. The thermal length describes the distance over which quasi-particle excitations 
in a typical energy interval given by the temperature T are mutually dephased. The 
above length scales correspond to the respective energy scales T4 = hvp/d of the 
Andreev levels, the temperature T and the impurity self-energy h/r, all present in 
the Eilenberger equation ( |3.2| ). In the following we first consider a clean NS slab, 
i.e. the electrons propagate ballistically between the NS interface and the metal- 
vacuum boundary (r — > 00). The Green's functions in the normal layer, in the low 
temperature limit T4 A of interest, take the form ||Zaikin 1982| 



g = tanhx(d), 

n _ cxp[sgn(v x )x(d-x)] [6.L0) 
•> coshx(rf) 



where (sgn(u a .) = v x /\v x \) 

x ( x ) = ^ + ^/ dx'Ay(x'), (3.11) 
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(uj n > 0). The normal Green's function is constant as a signature of the extent of the 
quasi-particle over the whole normal layer. Carrying out the analytical continuation 
u n — ► —iE + and determining the poles of g(E) yields the energies of the Andreev 
states bound to the normal layer, 



d 



E « = ^H n +o)- e -7 / dx'A y (x'), n = 0,l,.... (3.12) 
2d 2 a J 

We recover the energies of Eq. (|2.19|) of the previous chapter (in the limit T4 < A), 
which are split in direction ±v y and shifted by the vector potential. Note the non- 
local field dependence of the bound state energies. 

From the anomalous Green's function we can derive the pair correlation function 
F(x) = i(^|(x)^ T (x)) which satisfies A(x) = -VF(x), 

F(x) = N rrT ^ (U n (x, v F )). (3.13) 

LO n >0 

We consider the case in the absence of fields A y = here. Neglecting the effect of 
the boundary at x = d we find 

Fix) « N — e- x/iN{T \ (3.14) 

x 

At high temperatures T ^> Ta, only the first Matsubara frequency u = nT con- 
tributes to the sum ( |3.13|) , F ~ exp —x/^n(T), while at zero temperature, the decay 
shows no length scale, F ~ 1/x. 



3.3 Magnetic screening 
3.3.1 Linear response kernel 

When studying the magnetic response, the task is to determine the current-field 
relation in the normal layer in functional dependence of the vector potential j[A(x)] 
and solve for the screening currents self-consistently with the Maxwell equation 

-d 2 x A y (x)=47ij y (x). (3.15) 

In this chapter we focus on the magnetic response linear in the field. While the 
ground-state wavefunction experiences a weak perturbation that can be treated 
linearly, the screening of the magnetic induction is strong and requires the self- 
consistent solution of the Maxwell equation. Therefore the discussion of screening 
relies on the full knowledge of the dispersive response function j(q) = K(q)A(q). 
While in a London superconductor the response K(q) — K — —c/AnXl is local, 
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implying a penetration depth A cff = A L = (mc 2 /47rn s e 2 ) 1 / 2 , in a Pippard supercon- 
ductor with £ > A [[Tinkham 1996|| with K(q) ~ Kq/(1 + q^o) the nonlocal response 



changes the penetration depth to A eff = yAf£o > A L . Here we derive the current 
response in functional dependence of the field for the NS slab geometry. We find 
that the current response in the proximity effect is nonlocal in general. 
The current functional linear in A(x) takes the generic form 

j y (x) = — J K(x,x')A y (x')dx' (3.16) 

of a convolution with the response kernel K(x,x'). In the following we derive 
K(x,x') for the NS sandwich shown in Fig. |3.1| , expressing it through the Green's 
function in absence of fields go- Inserting the appropriate Green's function go, we 
are able to describe the current in a normal metal layer with arbitrary impurity 
concentration, ranging from the clean to the dirty limit. To calculate the linear 
diamagnetic response, we separate the Green's function into a (real) zeroth order 
part and an (imaginary) first order part in A{x) (we drop the Matsubara index u n 
in the following), 

#(x, v F ) = go{x, v x ) + igi(x, v x , v y ). (3.17) 
In the absence of external fields Eq. ( |3.2| ) reduces to 

-v x d x g (x,v x ) = A(x) f (x,v x ) - f£(x,v x ) 

-v x d x f {x,v x ) = 2Q{x)fo{x,v x ) - 2A(x)g {x,v x ) (3.18) 
v x d x fo(x,v x ) = 2u(x)fl(x, v x ) - 2A(x)g (x,v x ) . 

We have introduced the effective frequency u(x) = u + (g (x)) /2t and pair potential 
A (re) = A(x) + (fo(x))/2r as diagonal and off-diagonal potentials, respectively. 
Equations ( |3.18| ) imply that fo(x,v x ) = fo(x, —v x ), as A has been chosen real. We 
consider the solution of Eqs. ( |3.18j ) as given in the following. The first-order parts 
obey the equations 

- v x d x fi(x, v x , v y ) = 2£u(x)f 1 (x, v x , v y ) - 2A(x)g 1 (x, v x , v y ) 

+2ev y Ay(x)fo(x,v x ), (3.19) 
v x d x fl(x,v x ,v y ) = 2uj(x)fl(x,v x ,v y ) - 2A(x)g 1 (x,v x ,v y ) 

+2ev y Ay{x)fl{ X, v x ), (3.20) 

where g\ = (fofl + fifo)/^9o follows from the normalization g 2 = 1. Eqs. ( |3.19| ) 
and (|3.20|) can be mapped to two uncoupled Riccati equations and integrated out 
formally in terms of go> /o ; an d fo- Determining the current (|3.(j|) we have extracted 
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the kernel K(x,x') appearing in ( |3.16| ). The derivation of the results given here is 
deferred to Appendix [C| see also [[Bclzig 1998|| . Introducing a 'propagator' 



m(v x ,x,x') 




A fx") 



-dx" 



f Q (x" , v x 

with the properties m(v x ,x,x') = m{v X) x' ,x)~ l and m(v x , x, x")m(v x , x", x') 
the linear response kernel takes the form 

l-u 2 /v 2 F 

■ uU 

U„>0 Q 



(3.21) 



Vp 

K(x,x') = e 2 v F N 7rT [ d u- 



U 



[l+g (x,u)}[l-g (x',u)} (3.22) 



Q(x — x')m(u, x, x') + Q(x' — x)m(—u, x, x') + m(—u, x, d)m(u, d, x') 



The propagator (|3.21|) determines the ratio of the current response at x to the 
field at x', thus containing the range of the current-field relations. The inverse 
decay length of the propagator is proportional to the off-diagonal part of the self- 
energy A. The factor 1 — g is a measure of the superfluid density, vanishing in the 
normal state Qq = 1. The three summands in ( |3.22| ) stem from the quasi-particle 
trajectories connecting x and x', either directly or by one reflection at the normal- 
metal boundary. 

For illustration we reproduce the current response of a half-infinite superconduc- 
tor. Setting d = the solution of the Eilenberger equation (|3.18| ) takes the simple 
for m g = u n /Q n , f = /J = A/Q n , where Q n = (A 2 + w 2 ) 1/2 , see (pT7|). Inserting 
in ( 3.22| ) we obtain the linear- response kernel 

Vp 

a 2 r , i-u 2 /v 2 F 



K s (x,x') 



e 2 v F N nT ^ — / <fu- 

UJn>0 n q 



(3.23) 



-(2n„+l/r) 



\x— x'\ 



+ e 



which describes the current response of a superconductor with arbitrary impurity 
concentration [[Abrikosov 1975|| , which here additionally includes the effect of the 
boundary. For fields varying rapidly spatially we arrive at a non-local current-field 
relation of the Pippard-type ||Pippard 1953|| , while for slowly varying fields the kernel 
can be integrated out in Eq. ( |3.16|) , producing the local result ( |3.8|) . We recall here 
certain generic features of this kernel, which are of importance below. In a clean 
superconductor (1/r A), the range is given by the superconducting coherence 
length £ ) while in a dirty superconductor (A ^ 1/r) the range is given by the 
mean free path / = v f t, and is thus in both cases (nearly) temperature independent. 
We discuss below how, in the proximity effect, the range of the kernel varies from 
infinity to / and £at(T), exhibiting a strong temperature dependence, which leads to 
non-trivial screening properties. 
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3.3.2 Clean limit 



The linear response in the clean NS slab illustrates the nonlocality of the current- 
field relation in the proximity effect. Inserting the clean limit solution ( [3.10|) (with 
A = 0) into (|3.22| ), we obtain straightforwardly fZaikin 1982|| , 



h 



1 



All\ 2 (T)d 



A(x)dx . 



(3.24) 



The temperature-dependent penetration depth that is given explicitly by 
47re 2 n 1 



1 



A1,(T) 



\ 2 ' 







(3.25) 



-2T/T A T~>T A 



Here n = kp/37f 2 is the normal electron density producing the London length typ- 
ically Xn <C d. The propagator trivially becomes m(v x ,x,x') = 1 implying an 
infinite kernel range. Interestingly, the thermal length £,n{T) does not enter the ker- 
nel, as one might expect from the analogy to a clean superconductor with range £o- 
The spatial extent of the current-field dependence is thus cut off by the geometry 
< x < d. Eq. ( p.24| ) produces a constant current over the normal-metal layer, de- 
pending on the average vector potential. Solving Maxwell's equation ( p,15| ) with the 
boundary conditions A y (x = 0) = (perfect superconductor) and d x A y (x = d) = H 
(applied magnetic field), we easily find for Aat(T) <C d, 

3H 



Jy\ x ) 



BJx) w —x - — 



4/0*0 



87Tcf 

3H H 

x — 
2d 2 

H 3H 2 

Y x + ld x - 



(3.26) 



The current j ~ H/ d screens the applied field on the geometric scale d. The magnetic 
induction B changes sign inside the normal metal approaching B(x = 0) = —if/2 at 
the NS interface: the field is overscreened. While the screening currents minimize 
the average vector potential to 



A(x')dx' ~ Hd 



d 2 ' 



the average magnetization 



4vrM 



d 4 ' 



(3.27) 



(3.2* 
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makes up only 3/4 of the value for an ideal diamagnet. Thus the nonlocality results 
in a less efficient expulsion of the applied field, with overscreening as its signature. 
Overscreening is observed at low temperature, as long as \n{T) <C d, the diamag- 
netic susceptibility leveling off at 4nx = M/H « —3/4. At high temperatures the 
susceptibility is suppressed exponentially, 

4tty = — e" 2T/TA (3 29) 

see Fig. |3]2] (clean limit). 



3.4 Sensitivity to impurities 

For an arbitrary impurity concentration the screening problem has to be solved 
numerically. First the Green's function in the absence of fields go is determined, 



according to (|3.18|) . For the discussion of a stable algorithm therefore we refer 



the reader to | Bclzig 1999j| . Second, the linear response kernel (|3.22j ) is evaluated 



and solved self-consistently for A(x) with Maxwell's equation. Depending on the 
relative size of the thermal length £n{T), the mean free path I, and the thickness 
d we distinguish various regimes of either ballistic or diffusive electron propagation 
which are shown in Fig. |3l| Note that the dividing line between ballistic and 
diffusive behavior is not given by the simple relation £at(T) = I (dotted line). Rather, 
the range of the current-field relation £jy(T), I, or oo is needed for an accurate 
characterization of the magnetic response. The numerical results for the magnetic 



susceptibility 4ttx as a function of temperature are shown in Fig. |3.2| . The different 
curves show the clean limit and the mean free paths l/d = 10 4 , 10, 1,0.1. Clearly, 
a finite impurity concentration has a strong influence on the susceptibility, even if 
I > d, and can either increase or decrease the diamagnetic screening, depending on 



temperature. In the following, we explain Figs. 3.2 and 3.3 in detail. 



3.4.1 Ballistic regime 

The ballistic regime is limited by I ^> min{<i, £at(T)}, which ensures a ballistic 
propagation of the electrons over the thickness or the thermal length of the normal 
layer, respectively. As a limiting case, we have already considered the clean limit 
(Z->oo). ' 

An estimation using the zeroth order Eilenberger equation Eq. ( [OH ) easily shows 
that the clean-limit solution ( |3.10| ) remains valid in the absence of fields for 



I >dexp(2d/£v(T)), if £v(T)<d, (3.30) 
I » d, if 6v(T) > d, (3.31) 
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Figure 3.2: Numerical results for diamagnetic susceptibility 4irx of the normal 
metal layer (Xn = 0.003c?) . The clean limit is indicated by a thin line saturating 
at 4ttx = —0.75 for T < 5T^. The low temperature screening is enhanced by a 
finite mean free path / > d by reducing the kernel range from oo to I. Similarly, the 
high temperature signal is enhanced with decreasing mean free path. As I < d, the 
superfluid (screening) density is reduced, suppressing the diamagnetic susceptibility 
at all temperatures. The dirty limit is shown to overestimate the screening at large 
temperature. 



see Fig. |3~3| . Note that this includes the region d <ti I £ N (T) } the finite thickness 
preventing the small mean free path / -C £,n{T) from becoming effective. In the eval- 
uation of the kernel, the off-diagonal potential A = (/o)/2r enters the propagator 
( |3.21|) , yielding the kernel (|3.22| ), which approximately takes the form 



K(x, x') 



8TT\ 2 (T)d 



+ e 



(3.32) 



Since / 3> d, the exponentials may be expanded to first order. As a result, we obtain 
two contributions to the current j = j c i ea n + jimp 



-1 



Jclean 



47rA 2 (T)rf 



A(x)dx, 



(3.33) 
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Figure 3.3: Dependence of the magnetic response on thermal length £at(T) = 
vfIIkT and mean free path I. In the ballistic regime I 3> mm{d,^(T)} we 
distinguish three regions: (a) the clean limit with infinite range of the kernel ex- 
hibiting a reduced diamagnetism (overscreening) , (b) the quasi-ballistic limit with 
finite range £tv(T) increasing screening at large temperatures, (c) the ballistic limit 
where the finite range I enhances the screening although I 3> d. In the diffusive 
regime I <C l;N(T),d the range of the kernel is given by I. The dirty limit with 
nearly isotropic Green's functions is restricted to I <C £ ,!; N (T),d. Note that the 
current-field relations can still be local or non-local depending on the relative size 
of penetration depth and mean free path. 



Jimp \Xj 



8%\ 2 (T)d 



\x — x'\ + 2d — x — x 1 



I 



A(x')dx' . 



(3.34) 



According to fl3.27| ) solving the screening problem in the clean limit, we find j c i oa n ~ 
H/d. Taking ( |3.34|) as a small perturbation, we obtain j imp ~ Hd 2 / \%(T)l, which 
becomes of the same order of magnitude as j c iean if 



l~d 



d 2 



(3.35) 



As long as / 3> d 3 / X 2 N (T), the range of the kernel is effectively oo and the clean limit 
applies, see Fig. |3~3"| . As / < d 3 /\ 2 v (T) — which is the case for I ^> d in general - 
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the mean free path I enters the kernel range, as shown in Fig. |3.3j . The diamagnetic 
susceptibility is enhanced at low temperature (£at(T) ^> d), as Fig. |3]2] (/ = 10 4 <i) 
illustrates. The screening currents are even sensitive to a mean free path / <C d. 
The applied field is found to be screened on the effective penetration depth 

A eff = f/WTji < d, (3.36) 

in analogy to the effective penetration depth of a Pippard superconductor \ ef[ = 
f/\%(T)£ , see e.g. [ |lmkham 1996). 

In the remaining part of the ballistic regime, see Fig. [3.3| (range £jv(T)), the full 
solution is not known analytically, but we may produce an approximate solution, 
which characterizes well the numerical results at high temperature (£at(T) <C d). 
Limiting ourselves to high temperature with £,n{T) <C d allows us to consider the 
Green's function for the first Matsubara frequency ujq = ttT only. We restrict 
ourselves to the forward direction v x = +v F . From Eq. fl3.18|) we find that fo(v F ) ~ 
2 exp(— x/£n(T)) remains unchanged as in ( |3.10| ), and fo(v F ) <C 1 and 1 — go(v F ) <C 
1 obey the approximate equations, 

' ^fl{x,v F ) = -\{fo(x)) (3.37) 



dx £ N (T)J J{JK ' r> I 

— (1 - g (x,v F )) = —r(fo{x))f (x,v F ) 
with the solutions 

/.W) = ^e-l^) (3.38) 



(3.39) 



Here we have used (fo) ~ fo(v F )/2. fo(v F ) and 1 — go(v F ) are of order £n(T)/1 
rather than exponentially suppressed /q ~ exp — g?/£jv(T), 1 — g ~ exp —2d/^N{T) 
as in the clean limit ( [3.10|) . The correction to go given by Eq. (|3.39|) increases the 
superfluid density in the vicinity of the superconductor via the factor 1 — g in the 
kernel. The range of the propagator is modified by the correction (|3.38|) to /', 
leading to 

m{v F ,x,x') ^exp (^J^fjj . (3.40) 

The range of the kernel is now given by £n(T), which is strongly temperature depen- 
dent, and the current flows in a layer of thickness Cn{T) close to the NS interface, 

~ / dx'e~ 2X tNW ] A(x'). (3.41) 



Thus in the high temperature limit of Fig. |3.2| , as the mean free path I becomes 
smaller the diamagnetic susceptibility is enhanced. 
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3.4.2 Diffusive regime 



If impurity scattering dominates, as described by (go)/ r ^ ^ an d (/o)/ r 3> A, 
Eq. (|3.18| ) can be reduced to the Usadel equation fUsadel 1970|| for the isotropic 
part (fo(x)). Assuming T< A the solution in the normal metal takes the form 



(/o(a0) = cosh(< 



'-^(d-x))/ cosh(W^d) 



(3.42) 



where D = Vpr/3 is the diffusion constant. Equation ( |3.42|) shows that the im- 



portant energy scale is the Thouless energy E c = D/2ftd 2 [Edwards 1972j1 . The 



/-function ( p^2|) features the decay length f D (T) = {Djl-nTfl 2 reflecting the dif- 
fusive nature of the electron motion. 

In the normal metal I -C £at(T), d are necessary conditions for the Usadel theory 
to be valid. However, the numerical solution of (|3.18j ) shows that the Usadel theory 
in the normal metal, which requires the Green's functions to be nearly isotropic, is 
only a good approximation for I < £o : the dirty limit in Fig. |3.3| . 

Using the fact that the zeroth-order Green's function is nearly isotropic and 
varies on a scale C,d{T) ^> I, we find for the kernel ( |3.22| ) 



K(x,x>) = (/o(x)) 2 - 

N w„>0 



VF 



du 



4 



— u 



VpU 



+ e 



(3.43) 



The prefactor contains the full temperature dependence here and defines a space 
dependent penetration depth X(x,T). Qualitatively, the range of the kernel retains 
the same length scale I as in the adjacent ballistic regime, while the screening density 
is suppressed by the impurities. In Fig. |3.2| we thus find that a small mean free path 
/ <ti d produces an overall reduction of the susceptibility. 

For X(x, T) ^> / the vector potential may be taken out of the integral in Eq. (|3.16 ) 
and after integration over the kernel we recover the well-known local current-vector 
potential relation used in the Usadel theory HUsadcl 1970 ], with penetration depth 
X(x, T). We note that as long as X(x, T) < I, although the Green's functions are 
nearly isotropic, and in absence of the field are given by the Usadel theory, the 



current response is nevertheless nonlocal. In Fig. [3.2| the magnetic susceptibility for 
/ = O.ld is compared to the local Usadel result (dirty limit), showing that the latter 
overestimates the screening at high temperature. 

The high temperature (T ^> E c ) behavior is characterized by screening limited 
to a layer of thickness ^d(^) close to the superconductor, and consequently the 
diamagnetic susceptibility is given by 



Attx oc 



1 



d 



T 



(3.44) 
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This estimate is in agreement with the Ginzburg-Landau theory ||de Gennes 1966b ] 



as well as the numerical Usadel results [[Narikiyo 1989 



3.5 Discussion 

There are two main differences in the observable properties of the induced screening 
between the clean and the dirty limit. First, the zero temperature value of the 
susceptibility — Aitx saturates at 3/4 in the clean limit, while it reaches unity in 
the dirty limit. Second, the asymptotic behavior at high temperature differs: in 
the clean limit \ decays exponentially x ^ exp(— 2T/Ta), while in the dirty limit it 
follows an algebraic law x T -1 / 2 . The intermediate behavior illustrated in Fig. |3^ 
is determined by the competition of the nonlocality range and the screening density. 
With increasing disorder, a reduction of the range from oo to £tv(T) or / enhances 
the screening, while a reduction of the superfluid density weakens the diamagnetic 
currents. We have found several ballistic and diffusive regimes shown in Fig. |3l| 
with different characteristics. Remarkably, due to the absence of the small scale £o 
in the proximity effect, the nonlocal kernel (|3.22| ) is crucial to the understanding of 
all the regimes. 

Our results imply a non-monotonic dependence of the diamagnetic susceptibil- 
ity on the mean free path. Starting from a clean sample, with decreasing purity 
first the range of the linear response kernel is augmented and Meissner currents are 
enhanced, then the superfluid density is reduced and the screening currents are sup- 
pressed. Assuming a temperature-dependent scattering mechanism with decreasing 
mean free path as a function of temperature, such as electron-electron or electron- 
phonon interaction, we might speculate to observe a non-monotonic (i.e. re-entrant) 
behavior of the susceptibility (here the smallness of the scattering rate is compen- 
sated by the high sensitivity of the non-local current-field relation). However, as is 
evident from Eq. ( |3.21|) , the largest off-diagonal self-energy (A) which includes e.g. 



impurity scattering will provide a (low-temperature) cutoff for this behavior. 

Our results allows for a quantitative agreement with the experiments. The fit 
of the experimental data for the diamagnetic susceptibility gives an independent 
determination of the mean free path in these samples, which has been compared 
to the results of the resistivity measurements UMiiller-Allinger 1999|| . The magnetic 



susceptibility emerges as a sensitive indicator of the impurity concentration due to 
the nonlocality of the constitutive relation. 
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Chapter 4 

Magnetic breakdown in a 
normal-metal — superconductor 
proximity structure 



4.1 Introduction 



In this Chapter we extend the study of magnetism from the linear response to the 
finite field regime. Recent experiments have demonstrated the non-trivial screening 
properties of hybrid normal-metal-superconductor structures, exhibiting a magnetic 
breakdown at finite fields |0da 198Q , |Mota 1982| , |Bergmann 19871 , IVisani 19905 1. 
The investigated samples have typical dimensions comparable to the thermal length 
£jv(T) in the normal metal, attributing a key role to the quantum coherence of the 
electrons coupled to the macroscopic phase of the superconductor. 

We consider a clean normal-metal slab of thickness d in proximity with a su- 
perconductor, as shown in Fig. jjTl[ , The self-consistent study of the finite field 
effect in a dirty NS sandwich within the framework of the Ginsburg-Landau (GL) 
equation was carried out a long time ago by the Orsay group ||de Gennes 1966b| . 
Their work has provided the first understanding of the magnetic breakdown, which 
denotes the breakdown of the full magnetic flux expulsion in the normal layer, in 
analogy to the critical field of the superconductor. A quantitative agreement of the 
GL results with recent experiments is not achieved, as is discussed in [ [Mota 19"89 |. 
This is not surprising considering the role of nonlocality in the proximity effect, 
the local current-field relation of the GL theory applying only to the dirty limit. 
The quasi-classical Green's function technique offers the possibility to describe the 
opposite clean limit using the Eilenberger equations fEilenberger 1968| , which we 
pursue here. In ||Zaikin 1982|| the nonlinear response in the clean limit has been ob- 
tained along these lines. Using numerical methods, in fBelzig 1996a|| the non-linear 
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Figure 4.1: H — T phase diagram of the normal metal slab of thickness d. The 
breakdown field Hj,(T) ~ & /^n (T)d marks the first order transition between 
the diamagnetic and the field penetration phase (Ajv (T) denotes the penetration 
depth, $ the superconducting flux unit). The critical point at the intersection 
of the spinodals H SC (T) ~ <& /d 2 and H sh {T) ~ $ Q /A^(T) separates the first 
order transition for (T) <C d from the continuous cross-over at large temperature 
(Ajv (T) >d). 



field regime of the screening problem has been investigated, and the two (meta)- 
stable solutions in both the clean and the dirty limit have been determined (the 
dirty limit extends the GL results to low temperature using the Usadel equations 
|Usadel 1970|1 ). In this chapter we determine the H-T 'phase' diagram shown in 
Fig. |4.1| of the normal metal layer in the clean limit, where the bistable regime is 
particularly extended. In thermodynamic equilibrium, we find a magnetic break- 
down at Hi, (T), which is a first order transition separating the phase of diamagnetic 
screening from the phase of magnetic field penetration. The results presented here 
have been published in fFauchere 1997|| . 

Recent experiments [[Mota 1989| , |Mota 1994| have investigated the magnetic re- 
sponse of metallic cylinders with a superconducting core, finding data, which was 
claimed to be characteristic for the ballistic limit |[Mota 1994|| . This has motivated 
us to derive the analytic dependence of the clean limit expression for the breakdown 
field Hj, on temperature T and thickness d of the normal layer and to compare it to 
the experiment. From the free energy of the normal layer, which allows us to iden- 
tify the two (meta)-stable states, we determine the spinodals, the thermodynamic 
breakdown field Hi, (T) and find the jumps in magnetization and entropy at the tran- 
sition. Furthermore, we obtain a critical temperature which marks the upper limit 
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of the bistable regime (see Fig. [4. 1| ) . Finally, we relate our results for the breakdown 
field to the signatures of the nonlocality in the ballistic regime as they show up in 
the magnetic susceptibility \ and compare them with the experimental data. The 
following discussion is divided into four sections: the analysis of the constitutive 
relations (Sec. |4.2p , the solution of the magnetostatic problem (Sec. (4.3| ), the deter- 
mination of the breakdown field from the free energy (Sec. thermodynamics), and 
finally the comparison with the experiment (Sec. |4.5|) . 



4.2 Constitutive relations 

The quasi-classical Green's function technique provides an appropriate description 
of a metal with nearly spherical Fermi surface. In a finite magnetic field, the vector 
potential A (x) can be included as a phase factor along unperturbed trajectories, 
provided the dimensions of the normal metal are smaller than the Larmor radius 
(tl — hkpc/eH). In the ballistic limit, the quasi-classical 2x2 matrix Green's 
function g^ n (x, vp) satisfies the Eilenberger equation (e = \e\, h = ks = c = 1), 

- (vf • V) g Un (x, v F ) = [{u n + iev F • A(x))f 3 + A (x) fi, g Wn (x, v F )] , (4.1) 

discussed in Chapter [3|. We exclude impurity scattering here by assuming r — > oo. 
We consider a normal metal slab of thickness d on top of a bulk superconductor 



as shown in the inset of Fig. %% The vector potential A = (0, A (x) , 0) describes 
a magnetic field B = (0,0, B (x)) applied parallel to the surface, which induces 
screening currents j = (0, j (x) ,0). We make the usual idealizations in the de- 
scription of the NS sandwich: The superconducting order parameter follows a step 
function A (a;) = A6 (—x) (A real), no attractive interactions being present in the 
normal layer. We assume a perfect NS interface as well as specular reflection at the 
normal-metal boundary. 

In the subsequent analysis we restrict our attention to the magnetic response 
of the normal layer. In the proximity effect, the macroscopic coherence of the su- 
perconducting condensate induces correlated electron-hole pairs in the normal layer 
through the process of Andreev reflection. The basic process consists of an electron 
traveling forward and a hole traveling backward along a quasi-classical trajectory as 
shown in Fig. [4.2| (at discrete energies, bound Andreev states are found along these 
trajectories). In the presence of a magnetic field, the area enclosed by the trajectory 
(see Fig. [4.2Q is threaded by the flux 



$ (a, &,(p) = j) A(x) ■ dx = 2tan$cos^y A (x) dx, (4.2) 

which can be expressed through the integral a = A (x) dx times a geometric factor 
due to the inclination of the trajectory (the spherical angles $ and ip parameterize the 
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direction of the trajectory with respect to the x-axis). The current carried along a 
trajectory depends on the phase factor <E> (a, <p) / $ acquired by the propagation of 
both the electron and the hole along the Andreev loop, and we arrive at an intrinsic 
non-local current-field dependence j (a). The total current is determined by the 
sum over the currents along the quasi-classical trajectories, see Q3.6| ), from which we 
obtain the current expression for finite fields [[Zaikin 1982|| , 

j (a) = dd dtp j (0, <p, <£> (a, 0, <p)) , (4.3) 
Jo Jo 

where (a n = 2u n d/vp cosfi) 

j <p, $ (a, ip) ) = -2ev F N T ^ sin 2 & cos ip (4.4) 

Wn>0 

A 2 sin27r$/$ 

X 2 ■ 

(uj n cosh a n + + A 2 sinh a n j + A 2 cos 2 7r$/$ 

The induced currents for each trajectory depend only on the flux $ modulo the su- 
perconducting flux quantum $ = irhc/ e, reflecting gauge invariance. At small fields 
(a/$ *C 1), the current response is diamagnetic for all trajectories and the proxim- 
ity effect produces screening currents in the normal metal. As the field increases to 
a/$ ~ 1, some of the more extended trajectories produce paramagnetic currents, 
since the reduced flux $ G [— $ /2, $ /2] they enclose becomes negative, and the 
net diamagnetic current response is reduced. As we reach large fields (a/$ 3> 1), 
the Andreev loops become mutually dephased due to a uniform distribution of the 
reduced flux. The associated currents are randomly dia- or paramagnetic and the 
net current vanishes. Note that the proximity effect, i.e., the existence of the An- 
dreev levels is not destroyed in this limit, leading to a finite kinetic energy of the 
currents induced by the magnetic field. 



4.3 Magnetostatics 

Owing to the independence of j on x, the Maxwell equation —d 2 x A (x) = Anj and the 



constitutive equation (4.3) combined with the boundary conditions A (x = 0) = and 
d x A (x = d) = H can be given a formal solution. We arrive at a parabolic dependence 
for 

A (x) = Hx + 4vrj (a) x (d - -J , (4.5) 

parameterized by a = A (x) dx, which in turn is determined through the self- 
consistency condition 

Hd 2 47T , s ,o , s 

a = — + —j(a) d\ (4.6) 
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The total magnetization Ai (per unit surface) is defined by 



d 



AnM = / dx {d x A (x) -H) = 2nj (a) d 2 . (4.7) 



Eq. (O) contains the essential physics of the problem: For small fields (a — > 0), the 
current j ~ —H/d linearly suppresses the magnetic induction to the B (0) — > —H/2 
at the NS boundary (overscreening) . The current is given by the linear response 
expression in this limit, 

Ha/ *- <<1) ~-^i(Trd a - (48) 

which depends on penetration depth A at (T) <C d, see (Eq. (|4.15| )). When inserted 
back into ( |4.6| ), the vector potential is found to be strongly suppressed to a ~ 
H\ 2 N (T), and we obtain a consistent diamagnetic solution (i.e., a/$ D <C 1) for 
fields up to H < $ /A^. At large fields, the current vanishes (j — > 0) and the 
magnetic field penetrates the normal layer. From Eq. ( |4.6|) we find a ~ Hd 2 /2, 
consequently this metallic behavior is expected down to magnetic fields H > $> /d 2 , 
as follows from the condition a/$ 3> 1 for the Andreev levels to be dephased. With 
$ /(i 2 -C $o/A^ the diamagnetic and field penetration solution coexist in the regime 
<& /d 2 < H < $ /A^. These simple estimates for the limits of the bistable regime 



elucidate the numerical data given in [Belzig 1996a 



4.4 Thermodynamics 

In the phase diagram of Fig. |4.1| the upper and lower bounds of the bistable regime 



found from the above mean-field analysis are identified with the spinodals of the 
transition, the super-cooled field H sc ~ & /d 2 and the super-heated field H s h ~ 
$ /A^(T). In the thermodynamic equilibrium, a magnetic breakdown occurs at 
an intermediate field, connecting the diamagnetic regime to the field penetration 
regime by a first order transition. In the following, we determine this breakdown 
field and the associated entropy and magnetization jump from the free energy. 

The energy (per unit surface) of the currents j (x) = —SF/ 8 A (x) is obtained via 
an integration over the non-linear current expression, 

F (a) = — j j (a') da' 
Jo 

= hv F N Tj2 / dpsintfcostf (4.9) 

(u n cosh ol u + a/co>2 + A 2 sinh a n j + A 2 



log 



(u n cosh a n + \Ju)n + A 2 sinh a n j + A 2 cos 2 7r$/$ 
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F (a) describes the difference in free energy between the metal layer under proximity 
and in the normal state. F (a) is a monotonous and strictly positive function, 
reflecting the absence of condensation energy in the normal layer, and expresses 
the cost of the induced proximity effect lying in the kinetic energy of the currents 
induced by the vector potential. The free energy T (T, M) is constructed by adding 
the electro-magnetic field energy and subtracting the vacuum field contribution, 

F(T,M) = F(a)+£dx (^^L - *\ . (4.10) 

We do not include the condensation energy and the kinetic energy of the screening 
currents in the superconductor. The field dependent term of the condensation energy 
might in fact be of the order of the free energy in the normal layer and would be 
expected to produce numerical corrections in the results, which are not accounted 
for by our idealized choice of the order parameter A (x) = A9(—x). The kinetic 
energy of the screening currents ~ if 2 A may be neglected. 

After a Legendre transformation, we obtain the Gibb's free energy 

G(T,H) = F(T,M)-MH 

= F(a)+ /'d,Mt£!l. (4.U, 



o 



8tt 



The field term in Eq. (|4.11 ) describes the work necessary to expel the magnetic field. 



The extrema of the free energy Q with respect to a reproduce the equation of state 



( f4.6| ). Fig. shows the free energy Q (H) as obtained from the parameterization of 
Q and H through a. The breakdown field Hb (T) is determined by the intersection of 
the free energies Q of the two (meta)-stable solutions. We note that this procedure is 
equivalent to the Maxwell construction in the magnetization curve M. = —dQ/dH 
of Fig. 

In the following, we consider the free energy ( |4.9| ) in the two temperature limits 
T = and Ta <C T < A and obtain (Ta = hvp/^nd), 

F T=0 (a) = -v F N Q T A i dd I dip sin$ cos 2 ■& {arctan [tan7r$/$ ]} 2 ,(4.12) 
2 Jo Jo 

rir/2 pn/2 

F T ^ TA (a) = Av F N T-f 2 (T, A) / d<d / dy?sin$cos$ 

Jo Jo 

exp f- 2 ) sin 2 7r$/$ . (4.13) 



Ta cos $ 

The finite value of the superconducting gap A is accounted for by the dimensionless 
parameter 7 (T, A) = A/ ( J A 2 + (ttT) 2 + ttT ) < 1. 
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Figure 4.2: Magnetization M. (H) and free energy Q (T, H) at a temperature T = 
1.5 Ta (T4 = vp/2iTd). The representation is universal in the thickness d. The small 
and large field branches represent (meta)-stable solutions describing the diamagnetic 
and field penetration phases, which overlap in the field interval H sc < H < H s h. 
The Maxwell construction determines the first order transition between the phases 
at the breakdown field (dashed line). Inset: Cross-section of the normal-metal 
slab in contact with the bulk superconductor. The flux enclosed by the quasi- 
classical electron-hole trajectory (Andreev loop) enters as a phase factor in the 
Green's function. 



The free energies of the two (meta)-stable states can be approximated by their 
asymptotic forms in the limits a — > and a — > 00. In the diamagnetic regime, the 
expansion in a/$ D up to quadratic order of (|4.12| ) and ( |4.13|) provide the result 



a 2 



(T)d' 



(4.14) 



Eq. ( f4.14j ) is valid both in the low and high temperature limits using the penetration 
depth 1/X 2 N (0) = 1/\ 2 N = (Anne 2 /m) at T = and 

^7 2 (T,A)^e- 2 ^ (4.15) 



\%(T) \ 2 N ' T 

for T ^> T4. Note that the derivative j = —OF /da applied to Eq. ( |4.14j ) produces 
the linear response constitutive relation of Eq. ( |4.8| ). The Gibb's free energy follows 
from Eqs. (|4.10|) and ([4.14|) , using the solution of the Maxwell equations, 

G (a < $ ) w J-H 2 d. (4.16) 

OZ7T 
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Eq. Q4.16D is dominated by the magnetization work necessary to expel the field, which 



is parametrically larger (by {d/X^ (T)) 2 ) than the kinetic energy of the currents. 



We note that the scaling in the free energy (|4.16|), Q ~ H d, is parametrically 



independent of the nonlocality range of the constitutive relation ( |4.4|) . 

In the field penetration regime we approximate the free energy by its asymptotic 
value at a — ► oo. In this limit we replace the strongly oscillating functions of $ in 
( f4.12p and ( }4.13| ) by their average value ((arctantan $) 2 ) = n 2 /12 and (sin 2 <E>) = 1/2 



and obtain 



■N 

3 „,„, $ 2 



Wtc 6 X N d 

The magnetization energy vanishes in this limit. The corrections to the free energy 
( [4.17D are of relative order ($ D /a) 2 . 



The magnetic breakdown field (T) is determined by the intersection of the 
two asymptotics of the free energy Q given by Eqs. ( (4.16[ ) and ( |4.17j ), 



H i [T»T A ) « — 7 (T,A) A- e -«»< T >. (4.19) 

We note three important features of this result: The temperature dependence is 
a simple exponential with the exponent d/£ (T) = T/Ta, where £jv (T) = vf/2tiT 
denotes the thermal length. The amplitude of the breakdown field scales inversely 
proportional to the thickness of the normal layer, ~ 1/d. In the limit T — » 
the magnetic breakdown field saturates to a value which is suppressed by the factor 
7r/6\/2 ~ 0.37 as compared to the extrapolation of the high temperature result. 
The formula for the breakdown field if& ~ <& / Xx{T)d is similar to the critical 
field H c ~ <3> /A£o of a Type I superconductor ||Blatter 1994|| , the geometric scale d 
replacing the superconducting coherence length £ - 

We arrive at the H—T phase diagram shown in Fig. Q[ The first order transition 
between the diamagnetic and the field penetration regime takes place between the 
spinodals H sc ~ & /d 2 < H b {T) < H sh ~ $ d /Aat(T) 2 which delimit the (meta)- 
stable regime. Their intersection marks the critical temperature 

T crit ^T A log{d/X N ), (4.20) 



where (T crit ) w d. Below T crit the penetration depth is small, Xn (T crit ) < d, and 
we observe a first order transition. Above the critical point T crit , where A at (T cr j t ) > 
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d, a continuous and reversible cross-over between the diamagnetic and field pene- 
tration regime is expected. 

The latent heat (at T 3> T4) of the transition follows from Eqs. Q4.16 ), ( |4.17[ ), 
and ( fclgp using S = -dQ/dT, 

TAS*^-^-H 2 b (T)d, (4.21) 

and is related to the magnetization jump 

4nAM « (T) d (4.22) 

via the Clausius-Clapeyron equation. 

In the derivation of the breakdown field we have used the asymptotic expansions 
of the free energies in a/$ and $ /a, respectively. Their quality at the transition 
point is determined by the range of overlap between the diamagnetic and the field 



penetration regimes in Fig. 42, which is governed by the parameter A^v (T) / d. In the 
diamagnetic phase, the corrections are of the order of (a/$ ) 2 ~ (H b X% (T) /$ ) 2 ~ 
(\%(T)/d) 2 , and similarly in the field penetration regime. The expansion thus 
breaks down at Xn {T) ~ d, which is the critical point of the transition line. We 
note that the total magnetization changes from its diamagnetic value Ai ~ H b d to 
the strongly suppressed value M. ~ H b d (X^ (T) / d) 2 at the transition, reflecting its 
strong first order character. 



4.5 Experiment 

The breakdown field has been measured in fairly clean Ag-Nb cylinders recently 



Mota 198"9"| , [Mota 199H|] . Let us compare our results from the clean limit theory 
with the experimental data. We neglect the difference in geometry, cylindrical for 
the sample and planar in the theoretical model. 

In Fig. 43| we show the two data sets for the breakdown field data obtained 
on heating and cooling a sample of thickness d = 5.5/xm exhibiting hysteresis. We 
note that the theoretical values of the super-cooled and super-heated fields H sc and 
H s h are not reached in the experiments, as the phase transition is smoothed by the 



breakdown of local (Andreev) current loops, see Fig. |47^. The data saturates at 
low temperatures, in qualitative agreement with our theoretical analysis. Given the 
electron density in Ag, n = 5.8 ■ 10 22 cm~ 3 (=>• A at = 2.2 • 10 -2 /im) and d = 5.5/im, 



the breakdown field is determined by Eq. ( |4.18| ) and ( |4.19|) . Due to the idealization 



of our model, which assumes a step function for the order parameter, as well as 
the difference between the planar and cylindrical geometry, we expect a numerical 
factor correcting the amplitude of H b . Making use of the scaling factor « 0.56 in Eqs. 
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Figure 4.3: Breakdown field Hj, (T) and linear susceptibility 4ttx from theory and 
experiment (the analysis applies to a Ag-Nb sample of thickness d = 5.5/im). The- 
ory: results of Eqs. (|4.19| ) and ( [4.24|) shown as solid lines; H b (T) is rescaled to fit 
the zero temperature value Hj, (0) (horizontal line) to the experiment. Experiment: 
data shown as solid dots, the dotted line is a guide to the eye ||Mota 1994|| . Note that 
the logarithmic slope of the breakdown field is reproduced precisely (dashed line), 
while the one of the susceptibility is much smaller than expected. 



( [4.18|) and ( |4.19| ), we calibrate the theoretical result to fit the zero temperature value 
H b (0), as shown in Fig. [4.3| . The theoretical prediction for the high temperature 
behavior then follows from Eq. ( |4.19|) and is shown as a solid line in Fig. |4.3| . Eq. 
( |4.19|) accurately reproduces the logarithmic slope —1/Ta of the experimental data, 
thus correctly tracing the signature of the Andreev levels. The amplitude of Hb (T) 
deviates from the data by the constant ratio m 0.64, which can be attributed to the 
presence of a barrier at the NS interface, see below. 

A further agreement between theory and experiment is found in the scaling of 
the breakdown field with sample thickness d, which was reported to be oc 1/d, in 
accordance with Eq. ( [4.19[ ) (the experimental study involved 10 samples | [Mota 19"9~4* [ 
with thicknesses ranging from d = 2.9 /mi to d = 28/mi []). Similarly, the critical 
temperature determined in the experiment [[Mota 1989|| exhibits the same scaling 
oc 1/d, in agreement with Eq. ( |4.20| ) (T4 oc 1/d). 



1 We acknowledge unpublished data provided by B. Miillcr-Allinger and A. Mota, ETH Zurich. 
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For comparison, we cite the dirty limit result, which follows from the GL equa- 



tions de Gennes 1966b 



H D (T) « 1.9-%-exp(-d/&), (4.23) 

where ^ = ^vfI/QtcT and Xr> is a fitting parameter for the penetration depth 
at the NS interface. Both the exponential dependence on c£/£d and the 

independence of the overall amplitude on d clearly deviate from the experimental 
data. 

The good agreement between the clean limit theory and experiment for the 
breakdown field does not trivially generalize to other physical quantities, however. 
In particular, the linear susceptibility x — -M-/H sensitively depends on the non- 
locality of the constitutive relation j (a) . From Eqs. ( |4.7|) and ( |4.8|) we obtain the 
susceptibility 

_ AtiM _3 1 
X ~ ~1T ~ 4 l + 3A^ (T) /rf 2 ' { } 

which exhibits a temperature dependence with the following characteristics: 4irx de- 
cays exponentially oc 1/ A^ (T) at large temperatures, twice as fast as the breakdown 
field. The susceptibility takes half its maximal value at 3X 2 N (^1/2) /d 2 ~ 1, which 
roughly coincides with the critical temperature T cr #. The logarithmic derivative at 
T — T1/2 is predicted to be x' (^1/2) lx (^1/2) — V-^a- Below the critical point, the 
susceptibility saturates as the penetration depth decreases below the sample thick- 
ness (Aat (T) < d). Due to the non-locality, the penetration depth drops out of the 
expression for 47r% « —3/4 and we are in the regime of over-screening. In Fig. |4.3| we 
show the linear susceptibility according to the clean limit predictions Q4.24J) (there 
is no fitting parameter). The experimental data fails to show the typical saturation 
of the susceptibility expected below the critical temperature. At low temperature 
the experimental value clearly exceeds the maximal diamagnetic value —3/4 found 
in the clean limit. The decay at large temperature is slower than the decay of the 
breakdown field, while Eq. ( [4.24] ) predicts a decay with twice the logarithmic slope, 
see Fig. [4 . 3| . This discrepancy finds a natural explanation in the different sensitivity 
of Hb(T) and x(T) to the degree of nonlocality in the constitutive relation, as we 
discuss in the following section. 

Let us consider the influence of an insulating barrier at the NS interface. The 
consequences of a finite reflectivity at the NS interface on the linear current response 
has been analyzed in fHigashitani 1995|j. Their results allow for the reflection co- 



efficient R to be included in the penetration depth A jv (T) by redefining the factor 
1r (Ta < T < A) = (1- R) / (1 + R), in Eq. ((LTJ); X N (0) = X N remains un- 
changed. Inserting the modified penetration depth into Eq. ( f4.24|) we obtain the lin- 
ear susceptibility. The additional factor 7 does not change the characteristic shape 
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of the susceptibility (saturation, logarithmic slope at Ty 2 , exponential decay), but 
only lowers the position of the half- value of x t° ^1/2 ~ l°g [d(l — R) /Ajv (1 + R))- 
The finite reflection does not remedy the qualitative discrepancy between the sus- 
ceptibility in theory and experiment, in consistency with the above considerations. 
Considering the structure of the equations we may expect the dependence on the 
reflection R to enter in a similar fashion into the breakdown field Hb (T). Eq. ( |4 . 24| ) 
inserted in Eq. ([4.1 9|) gives the high temperature behavior, while the zero temper- 
ature result of Eq. ( f4.18| ) remains unchanged. We fit the breakdown field data by 
using first an overall scaling factor needed to adjust Hb (0) and second, a finite 
reflectivity, which only enters at high temperatures. The fit of the high tempera- 
ture behavior provides us with an estimate of the reflectivity R ~ 0.21 of the NS 
interface, see Fig. [O (dashed line). 



4.6 Discussion 



We have calculated the clean limit expression for the breakdown field separating the 
diamagnetic phase and the field penetration phase by a first order transition. We 
have determined the spinodals, the critical temperature as well as the latent heat of 
the transition. 

In Chapter [5] we found that the solution of the screening problem crucially de- 
pends on the nonlocality range. The high sensitivity to impurities is a consequence 
of the self-consistency problem with the Maxwell equations. In the ballistic regime 
I > d the linear screening was found to be enhanced as compared to the clean 
limit, at low temperature due to the kernel range / and at high temperature due 
to the range £n(T). Here the clean limit result for the breakdown field has been 
found to be in good agreement with the experimental data on fairly clean samples 
(/ y d). We resolve this apparent discrepancy by showing the the self-consistent 
screening problem drops out of the derivation of the breakdown field: The break- 
down field Hb ~ &o/ \]y{T)d is obtained from matching the magnetization energy 
Q ~ H^d in the diamagnetic phase with the kinetic energy in the penetrating field 
Q ~ $q/A^(T). The magnetization energy overshadows the kinetic energy of the 
currents in the diamagnetic phase, which contain the solution of the screening prob- 
lem. Thus the breakdown field is remarkably stable towards a finite impurity con- 
centration. A correction to the breakdown field of order d/l is expected from the 
kinetic term in the field penetration phase at temperatures T ~ T&. We note that 
the high temperature corrections, which are relevant beyond T1/2 hardly affect the 
breakdown field, as T c and Ti/ 2 coincide. The inclusion of a finite reflection at the 
NS interface permits an accurate fit of the breakdown field and gives an estimate 
for the quality of the NS interface. 



Chapter 5 



Paramagnetic instability of 
Andreev electrons 



5.1 Introduction 

A normal metal in contact with a superconductor exhibits the phenomenon of prox- 
imity — the superconductor exports its coherent state across the interface into the 
normal metal. On a microscopic level, this phenomenon is described through the 
Andreev reflection of the normal-metal quasi-particles at the NS interface, con- 
verting normal- to supercurrent. Proximity superconductivity exhibits a rich phe- 
nomenology and has attracted considerable interest recently [[Kouwenhoven 1997|| . 
A particularly puzzling finding is the ultra-low-temperature reentrance observed in 
normal-metal coated superconducting cylinders | |Visani 1990a| |, where, contrary to 



expectation, the fully diamagnetic cylinder develops a paramagnetic response at 
low temperatures. Recently, it has been speculated that some novel kind of persis- 
tent current states circling the cylinder might be responsible for this phenomenon 
Bruder 199&| ], but closer inspection of the experimentally measurable quantities re- 



veals that the predicted effect is by orders of magnitude too small fFauchere 1999a 



In this chapter, we demonstrate that the presence of a repulsive electron-electron 
interaction in the normal metal naturally leads to the appearance of a paramagnetic 
instability at very low temperature, offering a possible explanation of the reentrance 
effect in the NS cylinders. 

For simplicity, we consider a clean normal-metal slab of thickness d (0 < x < d), 
in perfect contact with a bulk, conventional superconductor. The proximity effect is 
mediated by the Andreev reflection at the interface with the superconductor, which 
binds the quasi-particles states to the normal layer for E < A5. In the usual free 
electron gas description of the normal metal, the Andreev bound states are found 
at E n = hv x {2n + l)n /Ad {n = 0, 1, ..; v x = Vp cos-#) producing a linear suppression 
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of the DOS ||de Gennes 1963| , [Saint- James 1964" | 



N(E) ~ N Ed/hv F , (5.1) 

close to the Fermi level E = (A/q = rnkp/h^ir 2 ). In the following we assume that 
the electron-electron interaction in the normal layer, which follows from the delicate 
balance between the phonon-mediated- and the Coulomb-interaction, is repulsive. 
As a consequence, a finite order parameter A(x) is induced in the metal, opposite in 
sign as compared to A s in the superconductor, see Ref. ||de Gennes 1964|| . The NS 



junctions then behaves like a Josephson junction with a phase difference n, trapping 
quasi-particle states at the Fermi energy close to the NS interface. The local density 
of states N(E, x) exhibits a peak at zero energy on top of the Andreev density of 
states, as shown in Fig. |5.1|. This peak involves a macroscopic number of states with 



density n p ~ k F / d, which in the following we call the 7r-states. 

The change in the DOS crucially affects the response of the proximity metal. The 
linear current response j[A] can be divided into two contributions j = jdia+ jpara, the 
diamagnetic current jdia = — (e 2 n/mc) A giving the rigid response of the bulk density 
n = k F /37r 2 and the paramagnetic current j pa ra following from the deformation of 
the wavefunction at the Fermi surface ||Schriefier 1988|| , 



i^=*i A OtSl (5.2) 

Jpara mc J V dE J N K J 

for slowly varying fields A (/ is the Fermi occupation number). While in a bulk 
superconductor the paramagnetic current is quenched by the energy gap at low tem- 
peratures producing a net diamagnetic response, the paramagnetic current of a bulk 
normal metal cancels the diamagnetic current exactly. In the non-interacting metal 
under proximity, the linear density of states suppression N(E) oc E is still sufficient 
to suppress the paramagnetic current at zero temperature [ JZaikin 1982| |. Including 
a repulsive interaction places the system in the opposite limit: The sharp DOS peak 
at the Fermi level produces a paramagnetic signal which over-compensates the dia- 
magnetic response. Such a paramagnetic response naturally leads to an instability: 
The free energy SF = —cj5A < can be lowered via a non-zero magnetic induc- 
tion induced by spontaneous currents along the NS interface. The interface currents 
are associated with an orbital magnetization M(T) producing a low-temperature 
reentrance in the magnetic susceptibility. 



5.2 Zero energy bound states 

In the following we present a quantitative analysis of the paramagnetic instability 
induced by the 7r-states. The magnetic induction B z (x) parallel to the surface 
is described by the vector potential Ay{x) which drives the currents jy{x ). The 
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Figure 5.1: Local DOS N(E, x) at the NS interface x = and at the metal boundary 
x = d (Ea = hvp/d), as it follows from the self-consistent solution of the Eilenberger 
equation, Eqs. ( |5.3D and (|5.4J) , for a thickness d = lOhvp/As and the coupling con- 
stants Vs = —0.3 and Vn = 0.1. Inset: Spatial dependence of the order parameter 
A(x) and local DOS N(E = 0, x) at the peak energy (S: x < 0, N: x > 0). 



electron-electron interaction in the superconductor is accounted for by an effective 
coupling constant Vs < and similarly Vn > in the normal metal, see also 
|Zhou 1995 , [Nazarov 1996|| . Two self-consistency problems have to be solved: First, 



we evaluate the order parameter A(x) accounting for the different coupling constants 
in the superconductor and the normal metal, and obtain the local DOS N(E,x). 
Second, we determine the current functional j[A] which we solve together with 
Maxwell's equation to find the spontaneous interface currents. 

We use the quasi-classical description following from the Eilenberger equation, 
see O, 

- v x d x g = \{uj n + iev y A y (x)}f 3 + A(x)f 1 ,g] , (5.3) 



where the 2x2 matrix g have been defined in Chp. |3.1| . Eq. (|5.3| ) is completed by 
the self-consistency relation for the pair potential ((..) is the angular average), 

A(x) = ~VN nTj2(f^(x,v x )). (5.4) 

^n>0 
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The self-consistent numerical solution of Eqs. ( |5.3| ) and ( |5.4| ) is shown in the inset 
of Fig. |5.1] . The course of the order parameter in the normal layer is asymptoti- 
cally given by A(x) ~ —V^Nohvp/x, as expected from the /-function in the non- 



interacting case Vn = UFalk 1963fl . A(x) decays from a value ~ — \Vn/Vs\A$ at the 



NS interface, to ~ — V^N^hvp / d at the outer boundary. Close to the NS interface, 
the local DOS 

N{E,x) = N Re[(g_ iE+5 (x,v x ))], (5.5) 



exhibits a pronounced peak at zero energy, as shown in Fig. |5.1| . At the outer 
metal boundary x = d the DOS is suppressed linearly, as in the noninteracting case 
dc Gcnncs 1963||. Such a pseudogap is typical in proximity induced superconduc- 



tivity |[Hara 1993| , |Frahm 1996|| . We note that a peak structure as a consequence of 



interactions could be related to one found for a Luttinger liquid in proximity with 
a superconductor [[Fazio 1996| . 



A (a;) 



In order to proceed with analytical results, we approximate the order parameter 
by a step function, 

A s , x < 0, 

— Atv, < x < d, 

where A^v oc Vn enters as a parameter. The Green's function in the normal layer 
x > can be determined exactly and takes the form 

_ u n sinh [x(d) - 7] + Ajy cosh [x(d - x)} 
9 ^ Vx) ~ fi n cosh[ X (d)- T ] ' [b - b) 

where x{ x ) = ^n,x/v x , = A 2 N + u 2 , and tanh7 = An/Qti (we consider the limit 



As ^> An, T). The second term in (|5.6| ) describes the ^-states at the NS interface. 
The poles of the Green's function at u n — > —%E + yield the bound state energies. 
For E > An the bound states given by 

hv t ( E 



E 2 — A 2 N = I nir + arccos —J , (n = 0, 1, ...), (5.7) 

down-shifted by 6E n m — 2Ajv/(2n + 1)tt with respect to the Andreev states of the 
free electron gas. Below the gap E < An we find the 7r-states at 



E = A N / cosh 2 ^ A " N EH ~ A N e- 2A » d ^ , (5.8) 

v x 

exponentially close to Fermi energy. All trajectories with v x = vpcosi) A^rci/ 
possess a bound state at E w 0, thus producing the macroscopic weight of the zero 
energy DOS peak: For A N > Vp/d the number of ^-states per unit surface iV sur f 
is equal to the number of transverse levels N sur{ ~ k F , while for A N < Vp/d it 
is reduced to iV sur f ~ k 2 F (A N d/vp) 2 via the reduction of the available solid angle 
cos?? < ANd/vF. 
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5.3 Spontaneous currents and magnetization 

We derive the current-field relations at low temperatures, assuming T <C v F /d, A N . 
This implies a thermal length £tv(T) = Hvf/^ttT larger than the thickness d and no 
thermal smearing on the scale An. Only the trajectories with cosi? < A^d/vF 
contribute to the current at low temperatures. We describe them in the limit 

v x j A N d -»• by 

g Wn ( X:Vx ) = ^ + UtiAn c~ x(x) . (5.9) 
Q n fl n (Q n — An) 

The current in the presence of a slowly varying vector potential A, follows from Eq. 
( |5.9|) after replacing uj n by u> n + iev y A and inserting it into the quasi-classical current 
expression, see (|3.6| ). In addition to the diamagnetic current 



Jdia = -(c/4vrA^)A (5.10) 

[Ajv = (mc 2 /47rne 2 ) _1//2 denotes the London length], we obtain the paramagnetic 
current 

c ^~/ n to 3$ ev F A , 

? D ara ~ — r^-e Kn a jr- arctan — — . 5.11) 

Jp 4vrA 2 v 2ir$f ttT ' v ; 

in the limit evpA/c,T An. Here, $o = vr/ic/e denotes the flux quantum, and the 
coherence length in the normal metal 

e N = hv F /2A N , (5.12) 

gives the extent of the 7r-states. Under the assumption An > vp/d we set a = 1. 
At temperatures T ^> eupA/c the paramagnetic current j para ~ (1/ \ N )(A N /T)A is 
linear in A and oc 1/T, a signature of the thermally smeared zero energy DOS peak, 
and competes with the diamagnetic current on the scale £ N . At T — > 0, Eq. ( |5.11| ) 
is nonlinear in the field and generates the spontaneous paramagnetic current. This 
paramagnetic interface current results from the energy splitting of the 7r-states in the 
field, E fa ±evpA, allowing the system to gain energy by shifting the DOS below the 
Fermi surface For An < hvp/d, the paramagnetic current is reduced by the factor 
a = (A N d/hvF) < 1 in Eq. ( |5.11| ). The surface current I = f jdx ~ a 2 c§o/\ N is in 
agreement with the current estimate I n ~ N sm fevF based on the number of 7r-states 
at zero energy. Eq. ( |5.11|) thus always produces a net paramagnetic response at low 



temperature and fields. 

The evaluation of the induced magnetization requires the self-consistent solution 
of Maxwell's equation — d^A(x) = 4irj(x)/c together with the current functional 

1 A similar mechanism producing spontaneous currents is discussed by Honerkamp et al. 



[Honerkamp 199S 
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Figure 5.2: Magnetization M(T, H) curve at various temperatures (in units of Ajy), 
for A a? = 0.3<i and £jy = d. The two meta-stable branches P± exhibit a spontaneous 
magnetization in zero field, the diamagnetic branch D is unstable. Inset: Zero field 
magnetization M (T). 



j[v4(x)]. The solution of the screening problem requires the full dispersive relation 
between j(q) and A(q). According to the study of Chp. |3], the coherence length in 
the normal metal £o gives the nonlocality range, as A^v is the off-diagonal self-energy 
in ( |3.2| ). Eq. ( |5.11[ ) represents the long wavelength limit q — > 0. For simplicity, here 
we use Eq. ( |5.11| ) under the assumption of local response. Including the nonlocal- 
ity in the constitutive relations does not alter the spontaneous interface currents 
qualitatively. 

The magnetization curves M(T,H), which follow from Eq. (|5.11|) are shown in 
Fig. |5.2| . Approaching from large fields, Fig. |5.2j shows two paramagnetic branches 
P± with a linear diamagnetic slope exhibiting a spontaneous magnetization in zero 
field. They result from the superposition of the paramagnetic magnetization M (T) 
and the Meissner response to the applied field H. As the field is decreased (in- 
creased) past H = 0, the branch P + (P_) becomes meta-stable. The spontaneous 
magnetization Mq(T) appears below a second order transition point T C M and sat- 
urates at low temperatures, as shown in the inset of Fig. |5.2|. The magnetization 
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curve includes a diamagnetic branch D, which arises from the competition between 
the paramagnetic instability and the thermal smearing and is thermodynamically 
unstable. 

In the following we give a semi-quantitative analysis of the magnetization M = 
J dx M(x)/d, first at zero temperature and field [M ], proceeding to finite temper- 
atures [M (T)], and finally including an applied magnetic field H [M(T, H)]. The 
boundary conditions are given by A(0) = and d x A(d) = H. We concentrate on 
the most relevant limit where ^,d A at- At T = 0, according to Eq. ( |5.11Q , 
the paramagnetic interface current j ~ a$o/ ^"n^n remains unscreened until being 
matched by jdia ~ — A/\ 2 N , producing a vector potential A ~ a§ /!;% on the scale 
Aat. The vector potential A saturates beyond A at, as the para- and diamagnetic cur- 
rents cancel each other. Assuming that the 7r-states extend up to the outer metal 
surface (a < 1), the induced magnetization M = A(d)/4:7rd is given by 

We note that although the spontaneous currents increase as A N > hvp/d (a = 1) 
they are screened exponentially beyond the extent of the ^-states in this limit, 
giving a magnetization M ~ ($ /£^<i) exp[— (d — £%)/\n]- We assume a < 1 in the 
following. 

At finite temperature, the spontaneous magnetization is suppressed by the factor 
arctan(et>irv4/7rT), which itself depends on the magnetization via A ~ Md, implying 
the implicit equation 

M (T) M (T)aA N 

— — — ~ arctan — — . 5.14 

M M T v ; 

The spontaneous magnetization appears below a second order transition at T C M ~ 
aA N , saturating at low temperatures, as shown in the inset of Fig. |57^. The tran- 
sition temperature is equal in magnitude to the energy splitting of the DOS peak 
E ~ evpA/c ~ aA N . 

Under an applied magnetic field H, the Meissner current jdia screens both the 
spontaneous interface current and the applied field. At zero temperature we deal 
with a linear problem and the magnetization is given by the superposition M(H) = 
Mo + xH of the spontaneous magnetic moment Mo and the Meissner response xH . 
As the temperature increases, Mo(T) decreases and the meta-stable regime shrinks. 
At T > T C M the spontaneous magnetization in zero field has disappeared, the sig- 
nature of the paramagnetic currents remains, however, reducing the diamagnetic 
susceptibility x & t small fields. At large temperature T 3> T C M we recover the pure 
Meissner response. 

Note that the two meta-stable branches P + and P_ in the magnetization curve, 
see Fig. p.'2[ imply a first order transition with changing field at H = 0. The 
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first order transition is similar to the magnetic breakdown occurring in the same 
system at large fields between the fully diamagnetic phase and a field penetration 
phase fFauchere 1997| |. The rotation of the magnetic moments to the energetically 
more favorable polarization will show the hysteretic behavior typical for a first order 
transition. The transition from P + to P_ implies a paramagnetic slope in the ther- 
modynamic cfc-magnetization (M(T, H)), which will link the meta-stable solutions 
P± in Fig. |5.2| and cross the origin at M(H = 0) = 0. In summary, we find that 
on approaching T C A/ from above, the diamagnetic susceptibility Xdc = (M(T, H))/H 
is reduced, exhibiting a low-temperature reentrance. Below T C M , the spontaneous 
interface currents produce a net paramagnetic susceptibility \dc- 



5.4 Discussion 

In the following we discuss our results in the context of the experiments by Mota 
and co-workers, who have measured the magnetic response of normal-metal coated 
superconducting cylinders at low temperatures jVisani 1990a| , [Mota 1989fl . In the 



previous chapters ^| and f| we have established a quantitative understanding of the 
magnetic response of these samples at higher temperatures, studying the screening 
at small impurity concentration and the magnetic breakdown at finite field. The 
samples are typically characterized by a mean free path I ~ d and an interface 
transparency of order unity. The Nb-Ag and Nb-Cu cylinders show an anomalous 
paramagnetic signal in the magnetic response in the low-temperature - low-field 
corner of the H — T phase diagram ||Visani 1990a] , |Mota 1991 . A direct comparison 



with our theory requires the magnetization curve (M(T,H)} which has not yet 
been measured. The observed dc-susceptibility XdciT) as a function of temperature 
shows an increase at low temperature fMota 1997fl . The measured ac-susceptibilities 



Xac(T) and Xac(H) exhibit a reentrance both as a function of temperature and 
field [[Visani 1990a| . The reentrance is accompanied by an out-of-phase response 



signaling dissipation and by hysteresis in the field dependence. These features are 
in qualitative agreement with our results for the magnetization curve. We find 
that theory and experiment agree in order of magnitude for a 2 ~ 0.1, implying a 
transition temperature T C M ~ lOOmK and a spontaneous magnetization Mq ~ 1G. A 
more quantitative comparison with experiment requires a self-consistent treatment 
of the spontaneous currents with the pair potential, accounting for the nonlocality 
of the current-field relation and its sensitivity to disorder. 

In conclusion, we have demonstrated that the inclusion of a finite electron- 
electron repulsion in a proximity coupled normal-metal layer naturally produces 
spontaneous interface currents leading to a paramagnetic reentrance in the mag- 
netic response: The sign change in the coupling across the NS interface leads to 
the trapping of 7r-states at the Fermi energy. The frustrated NS junction relaxes 
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through the generation of spontaneous interface currents, inducing a paramagnetic 
moment. The spontaneous magnetization implies a first order transition in zero 
field in the low-temperature sector of the H — T phase diagram. A non-trivial is- 
sue remains the requirement that the electron-electron interaction be repulsive at 
the low energy scales involved. Interesting consequences of this assumption have 
been discussed in the context of the proximity effect ||cle Gennes 1964j| and most 
recently in relation to the low temperature transport in mesoscopic NS structures 
[Pctrashov 1995| , [Nazarov 1996|| . In fact, the noble metal coatings used in the exper- 
iments of Mota and co-workers [[Visani 1990a|| appear to be the most plausible can- 
didates for a repulsive electron-electron interaction. Turning the argument around, 
in the light of our findings the experimental observation of a paramagnetic reen- 
trance can be taken as an indication of the presence of a repulsive interaction in 
these materials. 
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Chapter 6 

Nonlocality in Josephson 
junctions: Anomalous current — 
flux periodicity 



6.1 Introduction 



The equilibrium transport in Josephson junctions is caused by the quantum inter- 
ference of the two overlapping reservoir wavef unctions ||Josephson 1962|| . The su- 
percurrent is driven by the nonlocal phase difference between the superconducting 
leads. This is of particular interest in superconductor-normal-metal-superconductor 
(SNS) junctions, where the thickness of the interlayer can be much larger than the 
superconducting coherence length £ , without suppressing the Josephson current. 
The transport relies on the coherence of the quasi-particle population between the 
superconductors, on the mesoscopic scale of the thermal and phase coherence lengths 
of the normal metal IILikharev 1979 . 



In early studies on SNS junctions [ [Kulik 19T5|| the supercurrents were already 
expressed in terms of the quasi-particles which are bound to the interlayer due to the 
Andreev reflection. These bound states consist of a forward propagating electron 
and a back propagating hole and carry the current of the double electron charge 2e, 
which is converted from an to supercurrent at the NS interfaces. Recently, there has 
been a renewed interest in those states as they can be brought out of equilibrium 
in a controlled fashion ||Octavio 1983| , |Averin 1995| , |Gorelik 1998| , [Lehnert 1999|| . 
Due to the advances in nanofabrication technology, which have achieved a good 
coupling of a two-dimensional electron gas (2DEG) to superconducting reservoirs 
|Takayanagi 1995a| , |Takayanagi 1995T5| , ballistic properties of these quasi-particles 



like the quantization of supercurrent ||Beenakkcr 1991| have been observed 



In the present chapter we determine the critical-current-flux relation in a ballistic 
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SNS junction. The current- carrying quasi-particles traversing the weak link are 
sensitive to the Aharonov-Bohm phase, producing a nonlocal dependence of the 
current density on the magnetic induction in the junction and the superconducting 
phase difference. As a consequence, the nonlocality and the finite size junction 
produce an anomalous doubling of the critical-current-flux periodicity I c ($). 

The motivation for our work is drawn from an experiment on S-2DEG-S junc- 
tions |Heida 1998|1 of width w comparable to length d, where the a 2<3>o periodicity 
was found instead of the usual $o- A first attempt to explain this finding is due 
to fBarzykin 1998|| ; considering the point-contact geometry of Fig. 1(a) with open 
boundary condition in the metal, they indeed recover a 2$o periodicity for a geo- 
metric ratio w/d — > 0. However, the experiment in |[Heida 1998|| is carried out in the 
strip geometry of Fig. 1(b) and involves dimensions w ~ d of the same order. Here 
we determine the critical current I c as a function of flux, taking proper account of 
the reflecting boundaries in the normal-metal characteristic for the strip geometry 
of Fig. |Tl](b). We find that the periodicity of the critical current changes from <3> 
to 2$ as the flux through the junction increases, i.e., as a function of field. At low 
temperatures the crossover to the 2$ periodic current appears at a flux ~ Q w/d, 
thus explaining the result of Heida et al. |[Heida 1998|| , who found a 2$ periodic 
pattern for all fields in devices with w/d ~ 1. The results of this chapter have been 
accepted for publication [[Ledermann 1999j| . 



6.2 Current density pattern 

In our study, we neglect the screening of the applied fields by the supercurrent in 
the normal-metal layer, which is justified for a small critical current j c . The scale 
over which the fields can be screened is given by the Josephson length 

Aj ~ Vc^o/jc, (6.1) 

where d = d+2\$ is the effective length of the junction, extended by the penetration 
depth \s in the superconductors. We assume that the Josephson vortex distance 
ITinkham 1991 , 

a = $oM (6.2) 

is the smaller length scale, a <C Aj. Expressing the critical current by the screen- 
ing density n oc 1/\^{T) 2 introduced in Chp. [| j c ~ c§ / \ 2 N (T)d, we find that 
the field regime of interest lies above the 'critical field' of the normal interlayer, 
H 3> $o/AAr(T)(i. We consider the SNS-junction sketched in Fig. |6.1| . In the quasi- 
classical Green's function technique, the current density ( |3.6| ) in a point P results 
from the contributions over all quasi-particle trajectories connecting the two NS 
interfaces through P. In a junction of infinite width, the trajectories involve no 
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Figure 6.1: (a) Junction with a point-contact (open) geometry as discussed in [10], 
where w is the width of the two superconductors, (b) The junction studied here 
has a strip geometry with w the width of the normal conductor, (c) The magnetic 
field H is applied in the ^-direction and the coordinate system is chosen symmetric 
with respect to the junction center. The current density in the point P involves 
contributions from all trajectories T parameterized by the angle ip. (d) The phase 
difference 7 along the trajectory T can be expressed through the enclosed flux 0; 
using the trajectory T as our reference, the flux through the areas above (below) 
T contributes with a positive (negative) sign. 



reflection at the boundaries. In the case of a finite junction, boundary conditions at 
the normal- metal- vacuum boundary have to be applied, which we idealize through 
the assumption of specular reflections. Furthermore, we adopt the usual approxima- 
tions: perfect Andreev reflections at the SN-interfaces and a coherence length £ m 
the two superconductors with £ d, allowing for a step-like approximation of the 
order parameter A [[Kulik 19701 . The quasi-classical Green function is calculated by 
matching the partial solutions in N and S at the interfaces. For the current density 
j, we arrive at a generalization of the results given by Svidzinskii and co-workers 
Antsygina 1975||. For finite temperatures with d 3> £jv, the current j takes the form, 



3(x,y) 

Jc,0 



w/2 d 
d V p sin( 7 )— ===== 

-tt/2 V&vHV 9 ) 



cxp 



K<p) 



(6.3) 
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while in the low temperature limit, d <C £ N , 



Jcfi 



k=l 



k 



dip p sin(A;7) 



vr/2 



i{tpy 



(6.4) 



where p = (cos(y), sin(^), 0), Z(y?) = d/ cos(<£>) is the length of a trajectory with 
slope if, and the zero temperature critical current density is 



Jc,0 



ne 2 $ 
mc 2d 



(6.5) 



In ( |6.5| ), n denotes the electron density in the normal metal, and we have assumed 
T <C A as usual. While in the low temperature limit all harmonics sin(/c7) (k = 
1,2,...) contribute to the current density | )lshii 1970[ |, at finite temperatures, the 
thermal smearing of the Andreev levels leads to a suppression of the higher harmonics 
oc exp(— kd/^ N ) and only the first term oc sin (7) survives. An individual trajectory 
contributes with a weight oc exp(— l/£ N ) at finite- and oc d/l in the low temperature 
limit. In a wide junction, 7 takes the formF] [|Antsygina 1975|| , 



2-7T 

y; (p) = 70 - ir-Hd [y-x tan(^)] . 
<±>o 



(6.6) 



In the more general result derived here, 7 is given by the gauge invariant phase 
difference 



Aip / A • ds, 

% Jr 



(6.7) 



where Aip denotes the phase difference between the two superconductors and T is the 
path which goes through the point (x, y) with slope (p. Combining the current ex- 
pressions ( |6.3p or ( p\4p and fl6.7| ) with the Maxwell equation V 2 A = — 47rj[A, A<p\/c, 
we obtain the transverse vector-potential A allowing to solve the full screening prob- 
lem. The solution of the screening problem has only been carried out for a tunnel 
junction so far |pwen 1967 . 



In the following, we neglect screening and concentrate on junctions with the strip 
geometry of Fig. |6.1| (b), including the (reflecting) trajectories F in (|6.3| ) and ( |6.4|) . 
We express the gauge invariant phase difference ( p.7| ) in terms of the flux <fr enclosed 
by T and the reference path Tq and obtain, 



t(z> y;<p) = 70 



2vr0(x,?/;v?) 







(6i 



1 In the discussion of the point-contact geometry, Barzykin and Zagoskin make use of the 
above results, but limit the integration in (|6.3|) and (6.4) to those trajectories connecting the 



superconducting contacts [Barzykin 1998 
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Figure 6.2: The 'bow ties' of width r display the ensemble of trajectories contributing 
to the current density in a given point. The arrows indicate the current flow, (a) For 
weak nonlocality, r < a , the current flows straight through the junction, (b) The 
strong nonlocality, r > a , leads to the formation of vortex-like domains of circular 
flow. 



where for negligible screening <p(x, y\ <p) = HS(x, y; cp) and S is the properly weighted 
enclosed area, see Fig. p.l| (d). Note the analogy of this result to the Aharonov-Bohm 
phase picked up by an Andreev loop in a SN structure, see Fig. [|72|. The surface S 
is calculated as a function of the number of reflections the trajectory T undergoes 
(in the following called the 'order' of the trajectory). The point-contact geometry 
of Fig. |6.1| (a) then is described by the order-zero trajectories alone [Barzykin 1998|| , 



while in the strip geometry of Fig. |6.1| (b), higher orders have to be included. 

The geometrical pattern in the current density j depends strongly on the sample 
dimensions d and w, the normal metal coherence length £ N , and the applied field H. 
At finite temperature, the current density in P draws its weight from trajectories 
with ip < a/^jv/c/, allowing us to introduce the transverse nonlocality range r = y/£ N d 
(in the low temperature limit, tp ~ 1 and we define r = d). This range of nonlocality 
has to be compared to the scale a = &o/Hd of transverse variations in j (see 
Fig. |6.2|) : For weak nonlocality, r < a , the flow is uniform along x with amplitude 
j c and changes direction on a distance a /2 along the y-axis. This contrasts with the 
strongly nonlocal case r > a which is found with increasing field, where the current 
density forms domains of left-and right-going circular flow. While the local case is 
similar to that in a tunnel junction, the pattern in the nonlocal situation reminds of 
the usual vortex structure in a superconductor, see Fig. |Q| . For finite temperatures 
with d 3> £jv, the current density of the order-zero trajectories is given by 



3x(x,y) 

jc,0 



12 

sm ( 7q 

7T 

12 

7T 



2ny 



exp 



-a\x) 



cos 70 



a 

2ny \ £ N 2nx 
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d a 



exp [—at (a;) 



(6.9) 
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Figure 6.3: The current density for d/£ N = 3, r/a w 3.5, and $/$o = 3. Apart from 
the vortex-like domains in the middle produced by order-zero trajectories, additional 
circular flow is set up by the order-one trajectories. 

where 

For weak nonlocality, a < r, the exponent remains approximately constant in the 
normal part, a(x) ~ «(0), leading to a uniform current flow, while for strong non- 
locality, a > r, a(x) grows as x approaches the interfaces, a(±d/2) ^> a(0), such 
that the current concentrates in the middle of the junction. For a > r, the higher 



order trajectories lead to a refinement of the current pattern, see Fig. |6.3| . Similar 
results are obtained in the low temperature limit, see [[Ledermann 1998"|. 



6.3 Critical current — flux dependence 

The ratio r/a and its associated characteristic current pattern manifest themselves 
in the (pseudo)-periodicity of the critical current, 



/w/2 
dy jx(0,y;7o,*), (6.11) 
-w/2 



versus flux $ = Hdw in the junction. In the case of weak nonlocality, r < a , the 
relevant contribution to the critical current comes from the order-zero trajectories 
resulting in a $o periodicity. For strong nonlocality, r > a , higher orders are 
relevant and lift the order- zero result as shown in Fig. |6.4| — the periodicity of the 
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Figure 6.4: The critical current for d/£ N — 5 and w/d = 1/3. The solid curve shows 
the full critical current and the dashed curves are the contributions from the orders 
0, 1, and 2. The orders and 1 oscillate with periodicity 2<£> , while the second order 
decreases monotonically, remaining always positive. The current pattern produced 
by the orders and 1 is lifted by the order 2 contributions, and the critical current 
attains the periodicity 2<£> . 



critical current changes to 2$ - F° r a critical current j c < §o/cd the cross-over 
ao ~ r lies within the negligible screening regime. To be specific, we discuss in 
detail the orders 0, 1, and 2 for the case of finite temperatures with d >> £ N (the 
qualitative arguments for d <C £ N are similar). 

For w > r > a , the critical current due to the order- zero trajectories takes the 
form 

T«»(<b) = v / 24t w cos(7t$/$o) 

c 1 ' v^F r (7r$/<l>o) 2 ' 1 1 

where I c> t = w (12 / ir)j c ,o exp[—d/ £n(T)]. For the first-order trajectories, we numer- 
ically find a 2$ (pseudo-) periodic contribution as well. Both components vanish 
with field oc l/$ 2 . The second- and all following even-order trajectories exhibit a 
large current amplitude of order j c on a scale a oc 1/$ in the junction center (0, 0), 
a consequence of the (^-independence of the gauge invariant phase difference 7 along 
trajectories through (0, 0). Their contribution scales oc 1/$ and therefore dominates 
over the zeroth- and first-order terms at large enough fields — as the strongly non- 
local limit with a < r is reached, the periodicity changes to 2$ . Samples with a 
small width w < r are always in the strongly nonlocal limit and their current pattern 
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Figure 6.5: The critical current in the low temperature limit, d <C £ N , for w/d = 
0.8, 0.9, 1.1, 1.5 (from top to bottom). Successive plots have been offset by 0.35. 
Measured data ([9], diamonds) is shown for the case w/d = 0.9. 



is 2$o periodic throughout the entire field axis. At low temperatures, the condition 
w < r transforms into the geometric requirement w < d. The same arguments apply 
for the r > w > do, where the zero-order current is given by [[Heida 1998| 



/(■»($) 



I C>T w ( sin(7r$/2$o) 



2v/?r V 7r$/2$ c 



(6.13) 



While Eq. ( |6.13|) trivially implies a 2$o periodicity when taking the absolute value, 
more generally Eq. ( 6.12j ) only produces a 2$ periodicity when accounting for the 
lifting by the higher order contributions. We give a complete classification of the 
regimes of simple or double flux periodicity, depending on the nonlocality range r, 
and the Josephson vortex distance qq in Table |6.1|. 



6.4 Discussion 

We have demonstrated that the current density and the critical current in a clean 
SNS-junction with strip geometry depend on the ratio r/a between the nonlocality 
range r and the vortex distance a . The period of the critical current changes with 
increasing field from a $ periodicity for a > r to a 2$ pseudo-periodicity for 
strong nonlocality clq < r. 

Recently, Heida et al. [ Heida 1998|| observed such a 2$o periodicity in strip 



like (w ~ d) S-2DEG-S junctions made from Nb electrodes in contact with InAs 
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operating at low temperatures T = 0.1 K. As the total flux through the junction is 
difficult to determine in the experiment, Heida et al. had to infer their 2<3>o periodic 
structure from a fit on four samples with different ratios w/d ranging from 0.9 to 
2.2. In Fig. |6.5| , we present the results of our numerical calculations for the strip 
geometry, where we have properly taken into account the finite penetration depth 
of the flux into the superconducting banks. While geometries with w/d < 1 clearly 
exhibit a 2$o periodicity throughout the entire field region, a <3>o-component starts 
to develop at low fields in wide junctions. The comparison with the data of Heida 
et al. (w/d = 0.9) gives a satisfactory description of the pseudo-periodic structure. 

For further experimental studies on wider junctions, we predict a crossover from 
a $o- to a 2$o-periodicity with increasing fields. The cross-over is accompanied by 
the detachment of the current pattern form the superconducting leads, which form 
a vortex structure shown in Fig. 



Table 6.1: The periodicity of the critical current is controlled by the parameters 
w/r, and r/a . The table has to be read as a flow chart, starting at the top row and 
selecting the proper condition proceeding down the rows. The nonlocality range is 
given by r = ^^ N {T)d for £,n(T) <C d and r = d for ^at(T) 3> d. 
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Chapter 7 

Quiet SDS Josephson junctions for 
quantum computing 



7.1 Introduction 



Quantum computers take advantage of the inherent parallelism of the quantum 
state propagation, allowing them to outperform classical computers in a qualita- 
tive manner. Although the concept of quantum computation has been introduced 
quite a while ago ||Feynman 1996|| , wide spread interest has developed only recently 
when specific algorithms exploiting the character of coherent state propagation have 
been proposed |[Ekert 1996|| . Here we deal with the device aspect of quantum com- 
puters, which is flourishing in the wake of the recent successes achieved on the 
algorithmic side. Two conflicting difficulties have to be faced by all hardware 
implementations of quantum computation: while the computer must be scalable 
and controllable, the device should be almost completely detached from the en- 
vironment during operation in order to maximize phase coherence. The most ad- 
vanced propositions are based on trapped ions | |Cirac 19"9"o] , [Monroe 1995| |, photons in 
cavities [[lurchette 1995|| , NMR spectroscopy of molecules ||Gershenfeld 199?|j , and 
various solid state implementations based on electrons trapped in quantum dots 
[ Loss 1998 , Burkard 1999|| , the Coulomb blockade in superconducting junction ar- 
rays fShnirman 19971 , IShnirman 1998| , |Averin 1998|| , or the flux dynamics in Super- 
conducting Quantum Interference Devices (SQUIDs) [ [Bocko 1997|| . Nanostructured 
solid state quantum gates offer the attractive feature of large scale integrability, once 
the limitations due to decoherence can be overcome [[Haroche 1996|| . 

In this chapter we propose a new device concept for a (quantum) logic gate ex- 
ploiting the unusual symmetry properties of unconventional superconductors. The 
basic idea is sketched in Fig. 1: connecting the positive (100) and negative (010) 
lobes of a <i-wave superconductor with a s-wave material produces the famous 
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S' 



(010) ^ (100) 




Figure 7.1: Geometrical arrangements between s-wave and <i-wave superconductors 
producing a 7r-loop (used in the phase sensitive experiment by Wollman et al. [1]) 
and a qubit, the basic building block of a quantum computer. 



7r-loop with a current carrying ground state characteristic of <i-wave symmetry 
|Wollman 1993| , |Sigrist 1995|| . We make use of an alternative geometry and match 
the s-wave superconductors (S) to the (110) boundaries of the d-wave (D) mate- 
rial. As a consequence, the usual Josephson coupling oc (1 — cos0) vanishes due to 
symmetry reasons and we arrive at a bistable device, where the leading term in the 
coupling takes the form E^ cos 20 with minima at <fi = ±7r/2 (0 denotes the gauge 
invariant phase drop across the junction). In our design we need the minima at the 
positions <fi = 0, it — the necessary shift is achieved by going over to an asymmetric 
SDS' junction with a large DS' coupling, see Fig. 1. The static DS' junction shifts 
the minima of the active SD junction by the desired amount <p = ±7r/2. A similar 
double-periodic junction has recently been realized by combining two <i-wave super- 
conductors oriented at a 45° angle jlPichev 1998j| . The ground states of our SDS' 
junction are degenerate and carry no current, while still being distinguishable from 
one another: e.g., connecting the junction to a large inductance loop, the 7r state is 
easily identified through the induced current. It is this double-periodicity and the 
associated degeneracy in the ground state of the SDS' junction which we want to 
exploit here for quantum computation: combining the SDS' junction, a capacitor, 
and a conventional s-wave junction into a SDS' SQUID loop, we construct a bistable 
element which satisfies all the requirements for a qubit, the basic building block of 
a quantum computer. This work has been submitted for publication flotfe 1998 . 
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7.2 Device functionality 

In the following we give a detailed account of the operational features of our device. 
Consider a small- inductance (L) SQUID loop with IjL <C $ , where /, denotes 
the (Josephson) critical current of the loop and <£>o = hc/2e is the quantum unit 
of flux. Such a loop cannot trap magnetic flux ($ = 0) and the gauge invariant 
phase differences 0i and 02 across the two junctions are slaved to each other, as the 
uniqueness of the wave function requires that 0i — 02 = 27r<l>/$o- Combining a SDS' 
junction with a coupling energy E d and a conventional s-wave junction (coupling 
E s ) into a SDS' SQUID loop, we obtain a potential energy 

V{4>) = E d (l - cos 20) + E s (l - cos0), (7.1) 

exhibiting two minima at = 0, 7r, see Fig. 2. The switch s allows us to manipulate 
their energy separation, choosing between minima which are either degenerate or 
separated by 2E S . 

In the quantum case, the phase fluctuates as a consequence of the particle- 
phase duality [[Tinkham 1996|| . The phase fluctuations are driven by the electrostatic 
energy required to move a Cooper pair across the junction and are described by the 
kinetic energy 

T(0) = (£/2e) 2 C0 2 /2, (7.2) 

where C denotes the loop capacitance. The dynamics of is manipulated by in- 
serting a large switchable (switch c) capacitance C ext into the loop acting in parallel 
with the capacitances C d and C s of the d- and s-wave junctions. Note that the 
Lagrangian L = T — V of our loop is formally equivalent to that of a particle with 
'mass' m oc C moving in the potential U(0). 

With the switch settings c on and s off, see Fig. 2(a), the loop capacitance is large 
and the junction exhibits a doubly degenerate ground state which we characterize 
via the phase coordinate 0, |0) and \tt). Closing the switch s, see Fig. 2(b), the 
degeneracy is lifted and while |0) becomes the new ground state, the |7r)-state is 
shifted upwards by the energy 2E S of the s-wave junction, the latter being frustrated 
when = 7r. On the other hand, opening the switch c, see Fig. 2(c), completely 
isolates the ci-wave junction and leads to the new ground and excited states |±) = 
[|0) ± |7r)]/v2 separated by the tunneling gap 2A^. The latter relates to the barrier 
2Ed and the capacitance Cd of the c?-wave junction via ||Tinkham 1996 1 



A d KE d exp(-2^C d E d /e 2 ). (7.3) 

Closing the switch c, the capacitance is increased by C cxt and the tunneling gap is 
exponentially suppressed. Using the above three settings, we can perform all the 
necessary single qubit operations: 
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Idle-state: The switch settings c-on and s-off define the qubit's idle-state. While 
the large capacitance C ext inhibits tunneling, the degeneracy of |0) and \ir) guaran- 
tees a parallel time evolution of the two states. This idle-state is superior to other 
designs, where the two states of the qubit have different energies and one has to 
keep track of the relative phase accumulated between the basis states. 

Phase shifter: Closing the switch s separates the energies of the basis states 
|0) and \tt) by an amount 2E S . Using a spinor notation for the two-level system, 
the relative time evolution of the two states is described by the Hamiltonian 7i s = 
—E s a z , with <j z a Pauli matrix. Keeping the switch s on during the time t, the time 
evolution of the two states is given by the unitary rotation u z ((p) = exp(—ia z ip/2) 
with ip = —2E s t/h. 

Amplitude shifter: Assume we have prepared the loop in the ground state |0) and 
wish to produce a superposition by shifting some weight to the \ir) state. Opening the 
switch c in the loop, see Fig. 2(c), the time evolution generated by the Hamiltonian 
7i d = A d a x of the open loop induces the rotation u^i?) = exp(— io- x $/2) with 
$ = 2A d t/h. The system then oscillates back and forth between |0) and \ir) with 
frequency uj = A d /h and keeping the switch c open for an appropriate time interval 
t we obtain the desired shift in amplitude (note that the qubit remains isolated from 
the environment during these Rabi oscillations). 

Imposing the condition E d 3> E s , A d on the coupling energies, we make sure 
that the two states |0) and \tt) are well defined while simultaneously involving only 
the low energy states |0) and \n) of the system. Furthermore, all times involved 
should be smaller than the decoherence time Td ec , requiring E s , A d 3> h/r dcc . 

The present setup differs significantly from the conventional (large inductance) 
SQUID loop design, where the low-lying states are distinguished via the different 
amount of trapped flux and their manipulation involves external magnetic fields H or 
biasing currents / citeblatter:qtunneling. SQUID loops of this type are being used in 
the design of classical Josephson junction computers fLikharev 1991|| and have been 
proposed for the realization of quantum computers, too, see ||Bocko 1997|| . However, 
this setup suffers from the generic problem that the flux moving between the loops 
leads to a magnetic field mediated long-ranged interaction between the individual 
loops and further produces an unwanted coupling to the environment. By contrast, 
our device remains decoupled from the environment, the operating states do not 
involve currents, and switching between states can be triggered with a minimal 
contact to the external world — we therefore call our qubit implementation a ' quiet 
one. 

Next, we discuss how to perform two-qubit operations within an array of SDS' 
SQUID loops. A two-qubit state is a coherent superposition of single qubit states 
and can be expressed in the basis {\xy)}, where x, y 6 {0,7r} denote the phases 
on the rf-wave junctions of the first (x) and second (y) qubit, respectively. Unitary 
operations acting on these states are represented as 4 x 4 unitary matrices. Single- 
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(d) 

Figure 7.2: Energy-phase diagrams for the SDS' SQUID loop, (a) Idle-state: The 
switches are set to c-on and s-off — the relative dynamics is quenched, leaving 
the state unchanged, (b) Phase-shifter: With the switch settings c-on and s-on 
the relative phase between |0) and \n) increases linearly with time, (c) Amplitude- 
shifter: The switch setting c-off isolates the rf-wave junction. An initial state |0) 
oscillates back and forth between |0) and \ir), allowing for a shift of amplitude, (d) 
A SDS' junction, a n junction, and a s-wave junction combined into a SQUID loop 
and serving as a switch. 



qubit operations u acting on the second qubit take the block-matrix form 

and a similar block form selecting odd and even rows and columns defines the single- 
qubit operations on the first qubit. As all logic operations on two qubits can be 
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constructed from combinations of single-qubit operations and the Controlled-NOT 
gate [Ekert 1996|| it is sufficient to define the operational realization of the latter. 



The Controlled NOT gate performs the following action on two qubits: with the 
first (controller) qubit in state \x) and the second (target qubit) in state \y) the 
operation shall leave the target qubit unchanged if x — 0, while flipping it between 
and 7r when x = it, in matrix notation 

Wcnot = ( J ^ J ■ (7.5) 

The above Controlled NOT operation can easily be constructed from the two-qubit 
'phase shifter': Connecting two individual qubits in their idle-state over a s-wave 
junction into a SQUID loop, the states 1 00) and \-kti) become separated from the 
states |0tt) and 1 7r0) by the energy 2E Sb of the s-wave junction. Keeping the two 
qubits connected during the time t introduces a phase shift x — —2E Sb t/h between 
the two pairs of states, 



The Controlled NOT gate (|7.5| ) then can be constructed from the phase-shifter ( |7.6| ) 
via the following sequence of single- and two-qubit operations (see [ Loss 199~8fl for a 
similar realization of the CNOT gate), 

W CNOT = exp(-z7r/4)W 2 j / (7r/2)W 1 ,(-7r/2)W 2z (-7r/2)W ps (7r/2)W 2 ,(-7r/2), (7.7) 

where the single qubit operations Ui^{6) rotate the qubit i by an angle 6 around the 
axis fi {u^{6) = exp(—ia fl 9/2) acting on i) while leaving the other qubit unchanged. 

A key element in our design are the switches and a valid suggestion is the single 
electron transistor discussed in the literature [|Joyez 1994|| . Here we propose a quiet 



switch design optimally adapted to our SDS' qubits. The underlying idea is to insert 
a SQUID loop which allows to either stiffen or relax the phase, thus producing 
a phase switch. Combining a SDS' junction with energy E^, a 7r-junction with 
E n <C Ed-, and a s-wave junction with E s = E^ into a (small inductance) SQUID 
loop, see Fig. 2(d), we obtain the following switching behavior: The phase = 
on the SDS' junction frustrates the remaining junctions, the loop's energy-phase 
relation is a constant, E sw (cj) n = 4> s — ir) = 0, and the switch is open. A voltage 
pulse coming down the signal lines and switching the SDS' junction into the \ir) 
state changes the phase relation between the 7T- and the s-wave junctions and closes 
the switch: the energy E sw ((f) n = 4> s ) = 2E 7T (1 — cost/v) implies the current-phase 
relation / = (2e/h)d^E sw that stiffens the phase. The appropriate voltage pulses 
can be generated by driving an external SDS' SQUID loop unstable. 
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7.3 Discussion 

The quiet device concept proposed above heavily relies on the double periodicity of 
the SD junction. As the second harmonic is strongly suppressed in a SID tunnel 
junction, a more feasible suggestion for the realization of a cos 20 junction is the 
SND sandwich, where the superconductors are separated by a thin metallic layer N. 
For a clean metallic layer, the coupling energies for the n-th harmonic are large and 
of order E 3 ~ k 2 Ahv F /d, producing the well known saw-tooth shape in the current- 
phase relation [jlshii 1970|| (v F denotes the Fermi velocity in the N layer while d and 



A are its width and area). In reality, it seems difficult to deposit a clean metallic 
film on top of a <i-wave superconductor and we have to account for the reduction in 
the coupling E 3 due the finite scattering length / in the metal layer. Using quasi- 
classical techniques to describe a dirty SN D D junction, we obtain a second harmonic 
coupling energy E d ~ klA(hv F /d)(l/d) 3 ~ (R Q /R)(l/d)E T , where I denotes the 
scattering length in the normal metal, R Q = h/e 2 is the quantum resistance, and 
E T ~ (hv T / d)(l / d) is the Thouless energy. 

The second important device parameter is the tunneling gap Ad, which depends 
quite sensitively on the coupling to the environment. The usual reduction in the 
tunneling probability produced by the environment [|Caldcira 198311 is modified if 



the system is effectively gapped at low energies l|Ambcgaokar 1982|. This is the 



case for our SNDN'S' junction where the low-energy quasi-particle excitations in the 
metal are gapped over the Thouless energy E T ||Golubov 1988|| . The dynamics of the 



junction is only affected by the presence of virtual processes involving energies larger 
than E T , leading to a renormalized capacitance C ren ~ h/RE T (cf. [fSchon 1990| ) and 
resulting in the reduced tunneling gap 



A d (xE d exp[-is(R Q /R)^l/d}, (7.8) 

with v of order unity. Consistency requires that the tunneling process is 'massive' 
and hence slow, h/r < E T . With a tunneling time r ~ S/E d (S ~ H(R Q / R) sJYfd = 
tunneling action) we find that the constraint H/tE t ~ yljd < 1 is satisfied. 
The condition A^ <C E d requires the tunneling gap A^ to be small, but large 
enough in order to allow for reasonable switching times, requiring (R Q /R) \jljd 
to be of order 10. With typical device dimensions d ~ 1000 A, I ~ 10 A, and 
R/Rq ~ (d I 7) (1 / Ak 2 ,) ~ 1/100, this condition can be realized. Finally, the oper- 
ating temperature T is limited by the constraint S/h > E d /T, guaranteeing that 
our device operates in the quantum regime, and the requirement T < E T that 
thermal quasi-particle excitations be absent. The first condition takes the form 
T <ti h/r ~ y^l/dE T and is the more stringent one. Using the above parameters 
and a typical value v F ~ 10 8 cm/s, we obtain a Thouless energy E T ~ 1 K and hence 
T < 0.1 K. 
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In conclusion, we have discussed a novel device concept for logic gates in su- 
perconducting computers. The SDS' SQUID loop realizes a number of attractive 
features which are potentially relevant both in classical Josephson computers based 
on RSFQ logics as well as in superconducting quantum computers. The most obvi- 
ous advantage over previous designs is the quietness of the device: The SDS' SQUID 
loop is a naturally bistable device and does not involve external bias currents or mag- 
netic fields. Second, the basic states of the loop do not involve currents or trapped 
flux, hence long-range interactions between various elements of the computer are 
eliminated. Third, the qubits do not accumulate phase differences during idle time. 
And fourth, all operations can be carried out via simple switching processes. 



Chapter 8 

Conclusions and Outlook 



In this thesis, we have discovered the variety of facets of mesoscopic supercon- 
ductivity induced by proximity, ranging from the non-equilibrium transport to the 
thermodynamic screening properties. A tremendous progress in nanofabrication 
technology and the perfection of low-temperature cryostats in the last two decades 
have made the small scales of the Andreev or Thouless energies accessible to ex- 
periment. Proximity induced superconductivity combines the macroscopic phase 
of the superconductor to the quantum coherence of the microscopic electron wave- 
functions, allowing for the observation of single-particle properties on a macroscopic 
scale. The term 'mesoscopic' coins the importance of the phase coherence of the 
involved quasi-particles preserving quantum interference effects. 

The study of mesoscopic superconductivity has repeatedly lead us to the micro- 
scopic process of the Andreev reflection (AR) governing the quasi-particle exchange 
across the normal-metal-superconductor contact. The AR pins the single-particle 
physics to the Fermi level, combining electrons above and holes below symmetri- 
cally. The non-trivial AR phase-shift 7r/2 translates into the suppression of the 
quasi-particle spectrum at the Fermi energy, the resonance structure in transport 
and the diamagnetic response in magnetism. These phenomena are understood in 
the idealized free electron gas description of the proximity metal, in good quantita- 
tive agreement with experiments. 

In this concluding chapter, we give a perspective on open questions and future 
work on Andreev physics. On one hand we discuss some specific problems that can be 
addressed based on the framework provided in this thesis and are of direct relevance 
to experiments. On the other hand, we focus on some questions which appear of 
long term interest, in connection to our work and the present developments in the 
field of mesoscopic superconductivity. 

(i) The scattering matrix approach used in chapter ^| gives a direct theoretical 
access to the finite voltage shot noise in NS junctions. The numerical realization 
of disorder average based on the shot noise expression would allow for a quantita- 
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tive comparison with shot noise experiments on diffusive NS junctions presently in 
progress ||Jehl 199"9| , [Kozhevnikov 1999(1 - 



(ii) In chapters [5] and |6| we have found several signatures of non-locality in the 
proximity effect, which have just started to be traced experimentally ||Heida 1998|| . 
We predicted a finite field cross-over from a $o to a 2$ current-flux periodicity 
in ballistic SNS junctions, and the detachment of the current pattern as we reach 
the non-local regime. Future experiments on more extended ballistic junctions will 
require the self-consistent solution of the screening problem, where the over-screening 
and the sensitivity to impurity or boundary scattering crystallizes in the critical- 
current-flux relation. Non-locality is also an issue in High-T c superconductors, where 
the quasi- particle gap vanishes due to the d-wave symmetry of the order parameter, 
see ||Kogan 199§. 

(iii) The experimental feasibility of ballistic SNS weak links provides an access 
to non-equilibrium phenomena: The superconducting reservoirs act as a confin- 
ing potential on the Andreev quasi-particle population which can be driven out of 
equilibrium in a controlled fashion fScheer 1998|| . This allows for the observation 



of supercurrents at large temperatures [[Lehnert 1999|| or the manipulation of the 



Josephson relation by a repopulation of the Andreev levels Paselmans 1999| . The 
study of phase relaxation processes in this system are within experimental reach. 

(iv) The rapidly developing area of quantum computing relies on the existence 
of phase coherent devices as it exploits the time evolution of a quantum two- 
level system (qubit). The Josephson junctions offer themselves as a scalable, solid 
state implementation of a qubit, by making use of the quantization of the charge 
|Shnirman 1997|| or of the phase in mesoscopic weak links, see chapter |7|. Before 



speculating on the efficient use of quantum algorithms, the near future challenge 
lies in the experimental realization of single- and two-qubit operations and in the 
theoretical investigation of the phase coherence problem, as it arises from the cou- 
pling with the environment. 

(v) In chapters ^ and [| we have addressed the nonlinear magnetic response in 
proximity structures, revealing a rich behavior in the low-temperature low-field cor- 
ner of the H — T phase diagram. The breakdown of the diamagnetic state studied in 
Clip. [|was found to agree quantitatively with experimental data [ Mota 1989| . The 



polarization of the spontaneous moments discussed in Chp. |5] (see also below) imply 
a similar phase boundary at zero field. We predict a latent heat as a signature of 
both first order transitions. 

Much of the underlying motivation for the detailed study of the orbital mag- 
netism in this thesis was provided by the paramagnetic reentrance observed in 
the susceptibility of normal-metal coated superconducting cylinders ||Visani 1990a 



While the studies within the free electron gas approximation of Chps. ^| and [| could 
not account for the low-temperature anomaly, they allowed a good characteriza- 
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tion of the typical mean free path and the quality of the NS interface, see also 
IMiiller-Allinger 1999|| . The effect of the cylindrical topology had already been pre- 



viously considered in fBelzig 1995| | using the quasi-classical technique, as well as 
more recently in Pruder 1998|| , failing to produce either a non-monotonic tempera- 
ture dependence or the order of magnitude observed in experiment. This prompted 
us to consider the possibility of a repulsive electron-electron interaction in the nor- 
mal metal in Chp. y. We found that a repulsive coupling constant in the normal 
metal layer frustrates the NS contact by a phase difference ir across it. A density of 
states peak is accumulated at the Fermi energy due to the 7r-states trapped at the 
NS interface. This peak induces an paramagnetic instability towards spontaneous 
interface currents which naturally give rise to magnetic moments. The experimen- 
tally observed signatures being the paramagnetic reentrance in the susceptibility, a 
hysteretic behavior, dissipative response and creep qualitatively support our inter- 
pretation in terms of spontaneous magnetic moments and indicate the presence of 
repulsive electron-electron interactions in these material (Ag, Cu). In order to test 
our ideas, the following points merit consideration. 



(vi) The local spectroscopy of the NS contact by means of a (scanning) tunneling 
microscope should reveal a (split) density of states peak at the chemical potential. 

(vii) The experimental measurement of the <ic-magnetization curve would provide 
a direct evidence for our findings. For a quantitative comparison, the numerical 
solution of the Eilenberger equation is required, accounting for the self-consistency 
of the pair potential, the non-locality of the current-field relations that enter the 
Maxwell equation, and the presence of impurity scattering. 

(ix) Recent theoretical work also addresses the low-temperature corrections to the 
proximity effect due to interactions [|Zhou 1995|, [Nazarov 1996||. First experiments 



on ferromagnet-superconductor junctions show non-trivial behavior ||Giroud 199(S|| 
challenging our understanding within the free electron gas approximation. A den- 
sity of states peak at the Fermi level may affect the transport properties beyond the 
perturbative level that has been considered so far. 



As experiments progress to ultra-low temperatures we unravel new facets of An- 
dreev physics. The correlations induced by the proximity to a bulk superconductor 
could allow for the observation of metallic electron-electron interaction effects on 
a macroscopic scale. The low-temperature anomalies provide the opportunity to 
determine both the size and the sign of the coupling in these materials. 
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Appendix A 

Bogoliubov-de Gennes equations 



We diagonalize the BCS-Hamiltonian for superconductivity by a Bogoliubov trans- 
formation, considering a spatially inhomogeneous system. The transformation is car- 
ried out by determining the quasiparticle wavefunctions that fulfill the Bogoliubov- 
de-Gennes (BdG) equations. We show that the symmetry in the solutions of the 
(BdG) equations with respect to reversing the energies e Q — > — e a is a consequence 
of the spin degeneracy. Making use of the spin-reversal symmetry, we express the 
Hamiltonian, the density, and the current operators in terms of the quasiparticle 
operators j a , and arrive at the current expression used in chapter 0. 

The mean-field, spin singlet BCS Hamiltonian for superconductivity can be ex- 
pressed in the form (h = (— zV + eA) 2 /2m) 

H = /<fx:*t(x)( V _ A f + „)*(x): (A.1) 



using the Nambu spin-up annihilation operator 



*<"> = ! 1$). (A2) 

The pair potential is given by A(x) = A(^ (x) ^-f (x)) (coupling constant A), the 

colon (:) denoting normal ordering with respect to ^ and The Hamiltonian 
(|A.1|) can be diagonalized by a basis transformation, 

*(x) = ^$ q (x) 7q , (A.3) 

a 

7« = / d 3 x& a (x)*(x), (A.4) 
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with the eigenfunctions (x) of the BdG equations (which follow from the insertion 
of the Hamiltonian ( |A. 1| ) into [H, j a ] = —e a lc 



la) i 



h — fi A 
A* -h* + /j 



$ a (x) = e Q $ Q (x) , (A.5) 
<*>« (x) = 



Note that 7 a are annihilation operators of spin-up states and thus the excitations are 
all described in terms of spin-up quasi-particles. In order to preserve the (fermionic) 
commutation relations we need a complete orthonormal set of wavefunctions of the 
hermitian operator in ( |A.5|) , 



E*«( X )*1M = W(x-x'), (A.6) 

a 

d 3 x$i(x)$ Q ,(x) = 5 a , a , (A.7) 
(involving both positive and negative energy eigenstates) ||INiho 195^ . 



Consider the spin-reversal transformation S, 
SfcS" 1 ^^ (A.8) 

S linear (Sip^S -1 = tpi, SipiS^ 1 = — Noting that the order parameter A(x) 
is invariant under the transformation ( |A.8| ), it is easily seen that the Hamiltonian 
( [A. 1|) is spin-reversal symmetric, [H, S] = (this symmetry extends also to finite 
magnetic field, if the Zeeman splitting is neglected). By means of spin-reversal we 
may attribute to each quasi-particle operator j a a linearly independent operator 7^ 
through 

74 = S la S- 1 = J Sx *t (x) f J "J W (x) . (A.9) 
7„ describes an excitation with opposite energy 

£a ^a 

(according to [H, 7s] = 

£a7«)- From Eqs. ( |A.9j ) and ( |A.4|) we infer the effect of spin-reversal on the electron- 
hole wavef unction, 

*.<*>=(; )*••«=( ^ 

We arrive at a complete set of quasi-particle excitations in spin-up space, which 
are grouped into pairs 7^,7,5 with energies ±e a , as a direct consequence of spin 
degeneracy. Note that within the Nambu picture, all quasi-particles carry spin up 
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[spin | e and spin j /i (instead of spin j e)], explaining the opposite energy of the 
related wavefunction ([A.lOj) . The Hamiltonian (|A.1|) takes the form, 

U = Yl €a ( 7 - 7 « + ^ -2jd 3 x \v a \ 2 ^j . (A.ll) 

The ground state is realized by filling all the negative energy quasi-particle states. 

The spin-reversal symmetry allows us to express all equations using only half 
of the eigenstates, i.e., one representative $ n from each pair of states $ a }. 
In the following we choose the positive energy eigenstates, expressing the negative 
energy eigenstates through ( [A.lUp . We keep the positive energy states for the spin- 
up quasi-particles ( r y n = y n -\), and reinterpret the related quasi-particle states of 
opposite (negative) energy as spin-down excitations (7^ = 7L). The Bogoliubov 
transformation ( |A.3| ) then takes the well-known form 

V" ( U n \ „, . , ( ~K \ „/f 

V 
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*w=U 7 *t+ „• K> (A.12) 



and the Hamilton operator is expressed by TC = J2e n >o e ™(7nT7 n T+7fij7»U — ^ / d 3 x |t> n | 2 
displaying the spin degeneracy in the usual fashion. In the same way, we give the den- 
sity and current operators in both representations, using all (indexed by a) or only 

the positive energy eigenstates (indexed by n), respectively (u\7v = uVv — (Vu)v), 
p (x) = -e : (x) ( I ° ) * (x) : 



x -1 , 

e Yl { U *a( X ) u a'(^)lha' + <( X K' {*h*'ll) (A- 13) 



- e K( x K(x) - <(x)u n (x)] (^7,L T 7nT + 7,k i 7n| / 

£ra,£»>0 

-2e^K(x)| 2 (A.14) 

£n>0 

_e 5^ {k(x)« n (x) + w m (x)M n (x)]7 mT 7 fii + h.c. }, 



em,£n>0 



j(x) = --^:#t (x) v*(x): 

2m« 



e 



2mi 



2mi 



^<(x)VMa'( X )7i7a' -<(x)V^a'( X )7a'7i (A.15) 
J] ( M m( X )VWn(x) + <(x)V^ n (x)J (7^ T 7nT + 7^|7n|) 



,£m,en>0 
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-2 5>;(x)V«„(x) (A.16) 



€n>0 



+ [(-"m(x)Vfn(x) +V m (x)V« n (x)j 7mt7n| ~ t.C. 



>o 



The current operator is easily generalized to finite magnetic field. We note that 
current and density operators obey the continuity equation 

9P(X) + V • j(x) = 2.e*t ( X ) ( A(x) W (x) _ (A _ 1?) 



at v ' v y v ~ A *( x ) o 

Clearly, in the presence of the pair potential A the single quasi-particle currents are 
not conserved. However, as the right hand side of (|A.17 ) vanishes when taking the 
self-consistent expectation value, the total current is conserved, the quasi-particle 
currents being balanced by the current of the condensate. The Andreev reflection 
represents a good example, where the quasi-particle current of 2e entering the su- 
perconductor decays and is converted to supercurrent. 



Appendix B 



Quasi-classical Green's function 
technique 



Here we give give a short, but self-contained derivation of the Eilenberger equation 
fl3.2| ) and the corresponding quasi-classical current expression ( |3.6| ). We show how 
the renormalization procedure takes care of the cancellation of bulk diamagnetic 
and Fermi surface paramagnetic currents in the normal metal. Eq. ( |3.6| ) depends on 
the quasi-classical Green's function at the Fermi surface and its deviation from the 
normal state. The treatment given here follows the spirit of |Kopnin 1997 . 



B.l Eilenberger equations 

We define the 2x2 matrix Green's function in imaginary time [x = (x, r)] 

G M = -(T T * (ll )**(x 2 )> = ( ) , (B.l! 



by the Nambu-space operator ( A.2| ) in the Heisenberg picture with respect to the 



spin-symmetric Hamiltonian ( |A.1|) . G fulfills the equation of motion 



G- 1 *G(x 1 ,x 2 ) = 5( y x 1 -x 2 ). (B.2) 
where G' 1 is the matrix operator 

G _1 (xi,x 2 ) = - ( & " ^ + ^ ) 6( Xl - X2 ), (B.3) 

h = — g— (Vi + ieA(xi)) 2 , and the contraction is defined by A * B(xi,x 2 ) = 
J dx'A(xi, x')B(x', x 2 ). Similarly, the conjugate equation of motion takes the form 

G * G~ 1 (x 1 ,x 2 ) = S(x 1 - x 2 ). (B.4) 
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In the following we assume homogeneity in time d Tl = —d T2 . We also want to neglect 
the A 2 term in ( |B.3|) . While the Green's function G varies on the Fermi wavelength 
1/kp, the vector potential varies on a much larger scale, which we denote by (. With 

(V + ieA(x)) 2 = V 2 + ie (V • A(x)) + 2ieA(x) ■ V - e 2 A(x) 2 , (B.5) 

and typically V ~ kp, we find that the term e 2 A 2 ~ e 2 H 2 ( 2 is negligible in com- 
parison to 2eA • V ~ eHkpC under the condition 

C < ^ = r L , (B.6) 

where tl denotes the Larmor radius, which is the cyclotron radius of an electron 
traveling at Fermi velocity vf- 

The quasi-classical approximation is best done in the Fourier representation 

G(pi, p 2 ; r)=J d 3 Xl d 3 x 2 G(x ls x 2 ; T ) e -^+^^ ( B .7) 

where the equation of motion for the first Gorkov Green's function in ( |B.1| ) takes 
the form 

~ || + A*) G( Pl , p 2 ; r) - J d 3 k' A(k') P\ Pl - k', p 2 ; r) 

J d 3 k' k! • A(k') G( Pl - k', p 2 ; r) (B.8) 

/" rf 3 fc / A(k / } . (pi _ ^ _ k / r) = ( 27r )3 ( J(p _ p 2 ) 5 ( r ). 

m y 



Introducing the center-of-mass and relative coordinates 

x=-(xi + x 2 ), r = xi-x 2 , (B.9) 
and their conjugate momenta 

k = Pi-p 2 , P = ^(pi + p 2 ), (B.10) 

we find that r varies on scale 1/kp, while x varies on the much larger scale £■ Eq. 
QB.8|) translates to [G(k, p; r) = G(pi, p 2 ; r)] 



d (P + |) 2 
<9t 2m 



+ ^ G(k, p; r) - y rf 3 fc' A(k') F+(k - k', p - ^; r)(B.ll) 
"2^ / ^ / (2p + k-k')-A(k')G(k-k',p-|;r) = (2n) 3 6(k)6(r). 
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In the quasi-classical approximation, we retain only contributions for k' -C kp, 
assuming a smooth variation of A and A on Fermi wavelength (kp( <C 1), and 
obtain 

(~^~ P2 2m' k + /i ) Gfrnr) - J M*fflF*Qt-V,nT) (B.12) 
J d 3 k' p • A(k') G(k - k', p; r) = (27r) 3 <5(k)5(r), 



e 

m 



Transforming back to the real space center-of-mass coordinate and carrying out a 
Fourier expansion in time, we arrive at 

p 2 - ip ■ V e . . A . 

za;„ h fM p ■ A(x) G(x, p; u n ) 

2m m J 

-A(x)Ft(x,p;u; n ) = l. (B.13) 

The fast oscillations of G on Fermi wavelength are described by the p 2 term. The 
above procedure is carried out similarly for the Green's functions F and F\ obtain- 
ing the quasi-classical equation of motion for the matrix Green's function G and its 
conjugate 

G" 1 (x,p;cj n )G'(x,p;c; n ) = 1, (B.14) 
G(x,p;uj n ) G _1 (x,p;u; n ) = 1, (B.15) 

where (£ p = p 2 /2m — jj) 

v-i/ . \ • e _ * , ( ^ , ,P-V 



G (x, p; cuA = ?u;„ p • A + — £ p + z— — f 3 - ReAfi + ImAf 2 , (B.16) 

mc \ 2m J 

and the gradient V is to be taken with a negative sign when acting to the left in 
( |B . 1 5|) . We have thus arrived at first-order differential equations with respect to the 
center-of-mass coordinate x, the dependence of which accounts for the inhomogene- 
ity of the slowly varying potential and field. The momentum dependence of G lies 
at the Fermi surface |p| = kp. 

We are now ready to define the quasi-classical Green's function by the integral 
over the momentum energy = £ p , 

<L„(x,v F ) = -f 3 / ^G(x,p;w n ), (B.17) 



7T 



"id 



keeping the dependence on the momentum direction as denoted by the Fermi ve- 
locity vp. To provide a cutoff for the high energy contributions at E c , we use the 
integration path 70 shown in Fig. |B.1| (T, A Ec Ep). From the subtraction 
of the equation of motion ( |B.14p and its conjugate ( |B.15[ ), appropriately multiplied 
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with T3, we eliminate the £ p term in ( |B.16| ) and after integration over the momentum 
energy we find the Eilenberger equation 



-(vF-V)i„(x,V F ) = 

{u n + iew F ■ A(x)}f 3 + ReA(x)f 2 + ImA(x)fi , £ w „(x, v F ) 



(B.18) 



The quasi-classical equations need to be complemented by a normalization condition, 
as the source term on the right hand side of (|B.14j ) and (|B.15|) has been canceled 
by the subtraction. The normalization is provided by 

g 2 = 1, (B.19) 

as follows from the homogeneous solution and the fact that g 2 is a constant of motion 
of (|B.19|) . The following symmetry relations for the quasi-classical Green's functions 



and 



g{u n )* = -g(-Un), f(u n )* = P(-co n ), 



g(v F y = g(-v F ), f(v F y = f\-y F ) 



(B.20) 



(B.21) 



simplify the solution of (p.!9| ). Eq. ( |B.20| ) follows from the definition ( |B . 1 7|) and 
Eq. ( B.21|) directly from the equation of motion ( B.19|) . The quasi-classical Green's 
function takes the form 



9 



9 f 
P -9 



(B.22) 



For convenience, we redefine the pair potential A —>■ iA and obtain the form of the 
equations used in this thesis, see 



- (vf • V)g = {u n + iew F ■ A(x)}f 3 + ReA(x)fi - ImA(x)f 2 , g 



^B.23) 



B.2 Quasi-classical current expression 

The current density in terms of the usual Green's function (with respect to center- 
of-mass and relative momentum) is given by 

j( k ) = / 7^3 P G ^ P5 °") " — A ( k )- (B.24) 

The first term in ( |B.24| ) gives the paramagnetic current which is induced by the 
perturbation of the Green's function G in the presence of an applied field and lives 
at the Fermi surface. The second term in ( |B.24| ) gives the diamagnetic current giving 
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Figure B.l: Integration paths: 70 denotes the two inner half circles, the two 
outer high energy contributions. Together 70 and 7^ reduce to the integration over 
the real line. 

the rigid response of the bulk density n to the vector potential. In the normal state 
G = G^ n \ the paramagnetic current is known to cancel the diamagnetic current up 
to small correction of the order of Landau diamagnetism, see e.g. [|Tinkham 1996 1 - 
Using this fact, we easily find that the current 

j(k) = -— y d£ P N(Q J (G(k,p;0-)-GW(k,p;0-)), (B.25) 

can be expressed by the deviation of the Green's function from the normal state 
value (N(£ p ) is the density of states). Splitting the integration of the momentum 
energy £ p along the paths 70 and 700, we obtain 

j(k) = -eNoj^^j d£ p (G(k,p;0-)-G^(k,p;0-)) (B.26) 

-± f d£ p N(Q f^p(G-G<»)). (B.27) 

The high energy contribution in ( [B.27| ) vanishes as the Green's function G only 
deviates from at energies of the order £ p ~ T, A <C Eq- Inserting the Fourier 
expansion in time into the first we find 

j(k) = -eiV / J ^T^(G(k,p;a; n )-GW(k,p;a; n )). (B.28) 

Having introduced the high-energy cutoff by 70 we can exchange the Matsubara 
sum and the energy integration, which otherwise would have lead to a divergence, 
see [[Abrikosov 1975||. Using the definition of the quasi-classical Green's functions 



( |B.17| ), the current takes the form 

j(x) = ieNQixT J2 J Mx, y F ) - £(x, v F )) . (B.29) 



104 APPENDIX B. QUASI-CLASSICAL GREEN'S FUNCTION TECHNIQUE 



The quasi-classical Green's function for the normal state g = sgn(a; ri ) that is found 
from G(p) = \/{iu n — £ p ) conveniently drops out of Eq. (|B.29|) and we obtain the 
quasi-classical current expression (]3.6p. 



Appendix C 
Response kernel 



Starting from Eqs. ( |3.19| ) and ( p.2(J| ) we derive the linear response kernel K(x,x r ). 
The transformation of g, f, and /• introduced in ||Schopohl 1995| , [Schopohl 1998|| , 



a(x,v x ) 



f(x,v x 



a\x,v x ) 



(C.l) 



-v x d x a (x,v x ) = 2Q(x)a (x,v x ) + A(x)(al(x,v x ) 
v x d x al(x,v x ) = 2Q(x)al{x,v x ) + A(x)(al 2 (x,v x ) 
After linearization of ( |C.1|) we obtain 

a x (x, v x , v y ) - al(x, v x )a\(x, v x , v y ) 



l + g{x,v x y ' ' l + g(x,v x )' 

decouples the full equations of motions (|372|) to two Riccati differential equations. 
For the zeroth order equations we obtain 



(C.2) 

(C.3) 
(C.4) 

(C.5) 
(C.6) 



f 1 (X,V Xj Vy) 
fl(x,V x ,Vy) 



(1 + a (x,v x )al(x, v x )) 2 
a\(x,v x ,v y ) - a^(x,v x )ai(x,v x ,v y ) 



(1 + a (x,v x )al(x,v x )) 2 
and Eqs. (|3.19| ) and (|3.20|) are decoupled to 



- —a x a 1 {x,v x ,v y 



d x a\(x,v x ,v y 



u(x) + A(x)a (x,v x ) a 1 (x,v x ,v y ) 
+ev y A y (x)a (x,v x ) 

u(x) + A(x)al(x,v x ) a\(x,v x ,v y ) 
+ev y A y (x)al(x,v x ). 



As a consequence of Eq. (|C.5|) we find ai(—v y ) = —a\{v y ) and the same for a\, 
implying (fx) = (fl) = 0. Furthermore, since a\(v x ) = ai(—v x ), we only have 
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to consider one of the two equations (e.g. the first one). Equation ( |C.5[ ) is an 
inhomogeneous first-order differential equation, which can be integrated analytically. 
Assuming that / and f' do not change sign as a function of x, with the help of ( |3.18| ) 
the solution can be written as 

/ \ / .m(v x ,x,x ) 
a 1 (x,v x ,v y ) = c{v x ,v y )— r — (C.7) 



2ev,, 



v x fo(x,v x ) 



[1 - g (x', v x )]m(v x , x, x)A y (x')dx' 



.I'D 



where 



/ 2 f x> K(x") \ 
mfe,x,x') = exp^-y i J^'lz'') (&8) 

In this equation xq is an arbitrary reference point and the constant c has to be 
determined by the appropriate boundary conditions, m satisfies the relations of a 
propagator, m(v x ,x, x') = rn(u, x r , and m(v x ,x,x")m(u,x",x') = m(v x ,x,x'). 
Now we determine the constant c for a system of size [—d s , d]. We assume specular 
reflection at two boundaries at x — —d s , d and ideal interfaces between different ma- 
terials inside the system. The appropriate boundary conditions are the continuity of 
the Green's function along the classical trajectories, i.e., g(x, v x , v y ) = g(x, —v x ,v y ) 
at x = —d s ,d for the specular reflection and g(x = 0—,v x ,v y ) = g(x = 0+, — v x , v y ) 
at the NS interface. The same conditions are valid for a\ and a\. We thus obtain 



c(v x ,v y ) = 2e 



v y f m(v x ,d,x') + m(—v x ,d,x') 



v x J-d a m(v x , d, -d a ) - m(-v x , d, -d s ) 
x [1 - g Q (x', v x )]A(x')dx' . (C.9) 



The current is determined by Eq. ( |3.6| ), expressed by the solution ( |C.7|) . We obtain 
the following general result for the linear current functional, 



j y (x) = - J K(x,x')A(x')dx' , (CIO) 
where the kernel K(x, x') is given by 

V F 

du — [1 + go(x, u)][l — go(x', u)\ 







6(x — x')m(u, x, x') + Q(x' — x)m(—u, x, x') 
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m(—u, x, d)m(u, d, x') 
1 — m(u, d, —d s )m(—u, —d s , d) 

m(u, x, —d s )m(—u, —d s , x') 
1 — m(u, d, —d s )m(—u, —d s , d) 
m(—u, x, d)m(u, d, —d s )m(—u, —d s , x') 

1 — m(u, d, —d s )m(—u, —d s , d) 
m(u, x, —d s )m(—u, —d s , d)m(u, d, x') 
1 — m(u, d, —d s )m(—u, —d s , d) 



(C.ll) 



Equation ( |C.11| ) gives the exact linear-response kernel of any quasi-one- dimensional 
system, consisting of a combination of normal and superconducting layers extending 
from x = —d s to x = d. The kernel is expressed in terms of the quasi-classical 
Green's functions in absence of the fields, which may be specified for the particular 
problem of interest. We note two characteristic features of Eq. ( C.ll| ): The factor 
1 — go measures the deviation from the normal state Green's function go = 1, which 
produces no screening current. The propagator m(u, x, x') shows up in six summands 
which represent all the ballistic paths from i to i', accounting for multiple reflection 
at the walls at — d s and d. Thus the first two summands connecting x and x' 
directly constitute the bulk contribution, while the additional four summands are 
specific to a finite system (assuming specular reflection at the boundary). We note 
that a form similar to ( |U.11| ) may be derived for non-ideal interfaces between the 
normal and superconducting layers, if the appropriate boundary conditions following 
Zaitsev 19841 are taken into account. 



For the NS system considered in Chapter |3|, the linear response kernel (|C.11| ) 
may be simplified using m(u,x,— oo) — > and m(— u, — oo,x) — » as d s —* oo 
(u > 0), obtaining 



K(x, x') = e 2 v F N nT / du 

w n >0 n 



Vh - u 



[1 + g (x,u)] [1 - g (x',u)] 



Q(x — x')m(u, x, x) + Q(x' — x)m(—u, x, x') 



+m(—u, x, d)m(u, d, x') 
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